1、矩阵运算方面的题（参照第一章练习题3）
生成矩阵a,通过修改矩阵a的第1行第3列元素为2,把矩阵扩充为4行5列，其中第4行第3列元素为7，其余列为0,删除第2列元素的操作生成矩阵b, 按行的逆顺序取a的1,2,4,5列，且第四行为自然数1到4，生成新的矩阵c;计算b*c，b.* c, b+c,2*b, b和c的行列式，b\c, b/c,并找出a中大于3的元素，且将其替换为1。
上机实验步骤。
输入。>a=[0,1,0,2,1; 3,4,6,8,4; 9,7,3,2,2];b=a;
>b(1,3)=2;b(4,3)=7;b(:,2)=[
运行结果如下：
输入。>c=a(3:-1:1,[1,2,4,5])
运行结果如下：
输入。>c(4,:)1:4
运行结果如下：
输入。>b*c
运行结果如下：
输入。> b.*c
运行结果如下：
输入。> b+c
运行结果如下：
输入。> 2*b
运行结果如下：
输入。> det(c)
运行结果如下：
输入。> b/c
运行结果如下：
输入》 b\c
运行结果如下：
输入。> det(b)
运行结果如下：
2、绘图方面的题，包括二维和三维图形（参照第二章练习题）
1）、二维图形的绘制。
用ezplot绘制函数的图形。
2）、三维图形的绘制。
在同一平面中的两个窗口分别用polar绘制心形线和马鞍面z=x*x-y*y 的图形，并以不同的角度观察马鞍面。
上机实验步骤。
1）、在matlab命令窗口输入。
> ezplot('2*(t-sin(t))'2*(1-cos(t))'0,6*pi])
运行结果如下：
1）、在matlab命令窗口输入。
>clf
subplot(1,2,1)
theta=linspace(0,2*pi,30);
rho=3*(1-cos(theta));
polar(theta,rho,'b')
subplot(1,2,2)
x=-3:0.1:3;
y=x; x,y]=meshgrid(x,y);
z=x.^2-y.^2;
surf(x,y,z)
shading flat
运行结果如下：
三、m文件程序设计（参照第三章练习题1）
已知函数，计算f(-0.5),f(0),f(1),并作出该函数的曲线图，说明x=1.5,x=0,x=1处的连续性和可导性。
上机实验步骤。
1）在matlab命令窗口中，点击菜单file－>new->m-file，打开m文件编辑窗口输入程序如下：
function y=fun1(x)
if (x>=-1)&(x<0)
y=x+1;
elseif x<1
y=sin(x)+1;
elseif x<=2
y=x.*x;
end2）在m文件编辑窗口中，点击file－>s**e，在默认路径下存盘，文件名为。
3）在matlab命令窗口中，输入。
fun1(-0.5)，fun1(0)，fun1(1)
得f(-0.5)＝0.5000，f(0)＝1，f(1)＝1
4) 在matlab命令窗口中，输入。
plot(0,fun1(0),'1,fun1(1),'1.5,fun1(1.5),'
hold on
fplot('fun1',[1,2])
运行结果如下：
由图像观察得各点的连续性和可导性：
在x=1.5处，连续可导；
在x=0处 ，连续可导；
在x=1处，不连续不可导。
MATLAB结课作业
学号 20091110 姓名 李 班级 机械09 成绩 时间 2012 6 12 一任务。用迭代法求方程的一个根。二流程。考察方程 x x 是隐式不能直接求出根，但如果给出某个猜测值x0，将它带入上式的右端，即可求得x1 x0 然后又可取x1作为猜测值，进一步得到x2 x1 如此反复迭代得。xk 1...
matlab结课作业
北京林业大学全校公共选修课结课作业。matlab程序应用。2013年12月。1 指选题1 分别用矩阵求逆 矩阵除法求下面线性方程组的解 思路和方法 1 应用矩阵求逆的方法，解线性方程组。源程序 a 2 1 1 1 1 1 1 1 7 2 2 4 7 1 1 5 b 1 2 5 8 x a b 运行结...
matlab程序应用结课作业
北京林业大学全校公共选修课课程作业。matlab程序应用。2013年12月。1 分别用矩阵求逆 矩阵除法求下面线性方程组的解 问题描述 用矩阵除法解线性方程组。建模方法 x a b 程序 a 2 1 1 1 1 1 1 1 7 2 2 4 7 1 1 5 b 1 2 5 8 x a b运行结果 x ...