MatLab在工科数学中的应用

学校代码： 10128

《MatLab在工科数学中的应用》结课作业

学生姓名：梁文隆

学 号：201020509049 学 院：化工学院

专 业：无机非金属材料工程 班 级：2010级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。

?01021?A???34684??

??97322??

上机实验步骤 输入

>>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

运行结果如下：

输入>> B\\C 运行结果如下：

输入

>> det(B) 运行结果如下：

输入

>> det(C)

运行结果如下：

输入 >> B/C

运行结果如下：

二、绘图方面的题，包括二维和三维图形（参照第二章练习题）

（1）、二维图形的绘制

用ezplot绘制函数??x?2(t?sint),t?[0,6?]的图形。

?y?2(1?cost)（2）、三维图形的绘制

在同一平面中的两个窗口分别用polar绘制心形线??3(1?cos?)和马鞍面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 运行结果如下：