MATLAB语言及应用上机实习报告.docxVIP

MATLAB软件与基础数学实验 指导老师： 报告人： 学号： 班级： 报告所有代码均采用MATLAB2016a编写，和之前版本不同，它首次支持实时编辑，以下内容在可能在命令行窗口、编辑器或实时编辑器中编写 。 实验1 MATLAB基本特性与基本运算 求[12+2×(7-4)]÷32的算术运算结果。 (12+2*(7-4))/3^2 ans = 2 例1-2 计算5!，并把运算结果赋给变量y y=factorial(5) y = 120 计算2开平方 2^(1/2) ans = 1.4142 例1-4 计算2开平方并赋值给变量x（不显示） 查看x的赋值情况 x=2^(1/2); x x = 1.4142 例1-5 设，计算的值。 a=-24/180;b=75/180; sin(abs(a)+abs(b))/(tan(abs(a+b))) ans = 1.7951 例1-6 设三角形三边长为，求此三角形的面积。 a=4;b=3;c=2;p=(a+b+c)/2; s=(p*(p-a)*(p-b)*(p-c))^(1/2) s = 2.9047 例1-7 设，，计算，。 A=[1,2,3;4,5,6;1,0,1]; B=[-1,2,0;1,1,3;2,1,1]; A+B ans = 0 4 3 5 6 9 3 1 2 A*B ans = 7 7 9 13 19 21 1 3 1 norm(A) ans = 9.5758 A^(-1) ans = -0.8333 0.3333 0.5000 -0.3333 0.3333 -1.0000 0.8333 -0.3333 0.5000 例1-8 显示上例中矩阵A的第2行第3列元素，并对其进行修改. A(2,3) ans = 6 A(2,3)=1; A(2,3) ans = 1 例1-9 分别画出函数和在区间[-6,6]上的图形。 x=(-6)*pi:0.01:6*pi; y=(x.^2).*cos(x); z=sin(x)./x; plot(x,y); figure,plot(x,z); 例1-10 试求方程组的解。 a=[1,2,1;4,2,-6;-1,0,2]; b=[2;3;4]; x=a\b x = -30.0000 22.5000 -13.0000 例1-11 试求矩阵方程的解。 a=[1,2,1;4,2,-6;-1,0,2]; b=[1,2,3;1,1,1]; x=b/a x = 3.0000 -2.0000 -6.0000 2.0000 -1.5000 -5.0000 例1-12 建立同时计算，的函数。即任给a,b,n三个数，返回y1,y2. function mianfunction() disp(num2str(y1(1,2,3))) disp(num2str(y2(1,2,3))) function y=y1(a,b,n) y=(a+b)^n; end function y=y2(a,b,n) y=(a-b)^n; end end 运行显示： mainfuction 27 -1 例1-13 设，试画出在[0,2]上的曲线段。 % 加坐标网格 x=linspace(0,2); y=1./((x-0.3).^2+0.01)+1./((x-0.9).^2+0.04)-6; plot(x,y); grid on; 例如：对于例题1-13中所定义的f(x)，求其零点c. x=solve(1/((x-0.3)^2+0.01)+1/((x-0.9)^2+0.04)-6=0,x) x = -0.13161801809960646637049059278162 1.2995496825848217892317327263734 0.61603416775739233856937893320409 - 0.221903363220568094