matlab语言及应用上机实习报告

matlab语言及应用上机实习报告

ID:11650905

大小:740.68 KB

页数:50页

时间:2018-07-13

matlab语言及应用上机实习报告_第1页
matlab语言及应用上机实习报告_第2页
matlab语言及应用上机实习报告_第3页
matlab语言及应用上机实习报告_第4页
matlab语言及应用上机实习报告_第5页
资源描述:

《matlab语言及应用上机实习报告》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

1、MATLAB软件与基础数学实验指导老师:报告人:学号:班级:报告所有代码均采用MATLAB2016a编写,和之前版本不同,它首次支持实时编辑,以下内容在可能在命令行窗口、编辑器或实时编辑器中编写。实验1MATLAB基本特性与基本运算例1-1求[12+2×(7-4)]÷32的算术运算结果。>>(12+2*(7-4))/3^2ans=2例1-2计算5!,并把运算结果赋给变量y>>y=factorial(5)y=120例1-2计算2开平方2^(1/2)ans=1.4142例1-4计算2开平方并赋值给变量x(不显示)查看x的赋值情况>>x=2^(1/2);>>xx=1.4142例

2、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+Bans=043569312>>A*Bans=779131921131>>norm(A)ans=9.5758>>A^(-1

3、)ans=-0.83330.33330.5000-0.33330.3333-1.00000.8333-0.33330.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=ab

4、x=-30.000022.5000-13.0000例1-11试求矩阵方程的解。>>a=[1,2,1;4,2,-6;-1,0,2];>>b=[1,2,3;1,1,1];>>x=b/ax=3.0000-2.0000-6.00002.0000-1.5000-5.0000例1-12建立同时计算,的函数。即任给a,b,n三个数,返回y1,y2.functionmianfunction()disp(num2str(y1(1,2,3)))disp(num2str(y2(1,2,3)))functiony=y1(a,b,n)y=(a+b)^n;endfunctiony=y2(a,b,n)

5、y=(a-b)^n;endend运行显示:>>mainfuction27-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);gridon;例如:对于例题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.131618018099606466370490592781621.29954968258482178923

6、173272637340.61603416775739233856937893320409-0.22190336322056809485939175107611i0.61603416775739233856937893320409+0.22190336322056809485939175107611i例如:求一元函数最小值(fminbnd命令)由于fminbnd只能求一个最小值,通过观察函数可知极小值点为-∞和+∞>>x=fminbnd(y,0,inf)x=Inf>>y=@(x)1/((x-0.3)^2+0.01)+1/((x-0.9)^2+0.04)-6;>>[x,y]

7、=fminbnd(y,0,inf)x=Infy=-6>>[x,y]=fminbnd(y,-inf,0)%运行该命令时会报错,因为这里的x为NAN,这可能是因为fmindnd或optimget函数的代码有bug导致例如:求例题1-13中所定义f(x)在[0,1]上的定积分.>>symsx;>>int(1/((x-0.3)^2+0.01)+1/((x-0.9)^2+0.04)-6,0,2)ans=13*pi+atan(35641981496012/67550630868559)-12例1-14求二重积分及三重积分。>>symsxy

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。