数值建模与仿真实验报告

《数值建模与仿真实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、1、启动和退出MATLABo在命令窗口认识help、demo命令，并査找cell、floor、fix、round、rem、sign等函数的用法。2A=[l,0,-l;2,4，l;-2,0,5],B=[0,-l,0;2,l,3;l,l,2】求2A+B、A2-3B、A*B、B*A、A.*B、A/B、AB、A./B、A.B.•••••»B=[0,-l,0;2,1,3;1,1,2]0-10213112»A=[l,0,-l;2,4,l:-2,0,5]A=10-1241-205»C=2*A+BC=2-1-2695-3112»D=A*A-

2、3*BD=33-6213-2•15-321»5=A*BE=-1_2-293145710?=B*A-2-2_1-444-11410G=A.*B000443-2010H=A/B•2.00003.0000■5.0000•5.00003.0000-4.00007.0000-9.000016.0000»I=ABI=0.3333-1.33330.66670.25001.00000.25000.3333■0.33330.6667»J=A./BWarning:Dividebyzero.J=Inf1.0000-2.000004.00000•In

3、f0.33332.5000»K=A.BWarning:Dividebyzero.K=01.0000-0.5000■Inf0.2500Inf03.00000.40003、利用函数产生3X4阶单位矩阵和全部元素都是4.5的4X4阶常数矩阵»L=eye(3,4)100001000010»^=4.5*ones(4)4.50004.50004.50004.50004.50004.50004.50004.50004.50004.50004.50004.50004.50004.50004.50004.50004、利用函数产生5X5阶随机分布

4、的矩阵和5X5阶正态分布的随机矩阵»N=rand(5,5)N=0.95010.76210.61540.40570.05790.23110.45650.79190.93550.35290.60680.01850.92180.91690.81320.48600.82140.73820.41030.00990.89130.44470.17630.89360.1389»0=randn(5j5)0=-0.43261.1909-0.18670.11390.2944•1.66561.18920.72581.0668-1.33620.1253-

5、0.0376-0.58830.05930.71430.28770.32732.1832-0.09561.6236-1.14650.1746-0.1364-0.8323-0.69185、练习文件读写操作1、gcd函数用于求两个整数的最大公约数。先用help命令查看该函数的用法，然后利用该函数求15和35的最大公约数。»:g,c,d】=ccd<15,35〉2、己知矩阵A=[8,9,5;36,7,11;2V8,5],B=[-1,3，-2;2,O,3;-3,1,9]求下列表达式:(1)A+5*B和A-B(2)A*B和A.*B(3)A/B

6、和BA(4)AA3和A.A3(5)[A,B]72-633.1341-1.2B01-1.分8781O.7I9G-1.26S39.47S61.OS544.8O9r-1.

7、;»fork=l:100»=(2*k-l)*2;s=s+m;end»s13333004、求下而分数序列前20项之和:2358T235••…»ss=0；forn=0:10-_S—s=s+l/ga»»a(n):end32.66031、1行100列的Fibonacc数组a，a⑴=a(2)=l,a(i)=a(卜l)+a(i-2),川for循环语句来寻求该数组屮第一个大于10000的元素，并指出其位置h要求编写M文件。n=.00;a=oncs(l,n);fori=3:na(i)=a(i-l)+a(i-2);ifa(i)>10000a(i)

8、,break;ans=10946end;212、创建一个4X3阶的服从0-1均匀分布的随机矩阵A,求岀A中各行的平均伉B(列句量)，将B补在A的右方构成4X4阶的矩阵C,并提取C的下三角矩阵。•Mrand(4>3);•B=A(:,1)+AG,2)+AG,3>•O