MATLAB刘卫国版本Word文档格式.docx

1、0.34780.42060.4841Column49through510.53790.58150.6145Column52through540.63660.64740.6470Column55through570.63510.61190.5777Column58through600.53270.47740.4126Column610.3388t2z0t1t2411t2，其中t0:0.5:2.5t22t1,2t3程序为t=（0:2.5）;z4=（t=0&t=1&2）.某（t.2-1）+（2A=1234-4;34787;3657;B=13-1;203;3-27;I=100;010;001;A+6某

2、Ban=1852-10467105215349A-B+I1231-3328840671A3372262338244860424737014918860076678688454142118820A.3172839304-643930434365850327274625343A/B516.4000-13.60007.600035.8000-76.200050.200067.0000-134.000068.0000BA109.4000-131.2000322.8000-53.000085.0000-171.0000-61.600089.8000-186.2000A,B1234-413-1347872

3、0336573-27A（1,3,:B21234-436574511101920-54012345306789101763.设有矩阵A和B：A1112131415,B0231617181920972122232425413求它们的乘积C；将矩阵C的右下角32子矩阵赋给D；1694011查看MATLAB工作空间的使用情况。A=12345;678910;1112131415;1617181920;2122232425;B=3016;17-69;023-4;970;41311;C=A某BC=93150776258335237423520397588705557753890717D=C（3,4,5,2,

4、3）D=520397705557890717whoNameSizeByteClaAttributeA5某5200doubleB5某3120doubleC5某3120doubleD3某248double4.完成下列操作：求100,999之间能被21整除的数的个数；建立一个字符串向量，删除其中的大写字母。A=100:1:999;K=find（rem（A,21）=0）;length（K）an=43ch=aHDJLK143663CFHI;H=find（ch=A&chch（H）=ch=a143663实验二MATLAB矩阵分析与处理二（1）E=eye（3）;R=rand（3,2）;O=zero（2,3）;

5、S=diag（2,3）;A=E,R;O,S;A2B=E,（R+R某S）;O,S2an=1.0000002.85041.9439701.000000.69343.5652001.00001.82053.04840004.0000000009.0000B=1.0000002.85041.943901.000000.69343.5652001.00001.82053.04840004.0000000009.0000（2）H=hilb（5）P=pacal（5）Hh=det（H）Hp=det（P）Th=cond（H）Tp=cond（P）H=1.00000.50000.33330.25000.20000.

6、50000.33330.25000.20000.16670.33330.25000.20000.16670.14290.25000.20000.16670.14290.12500.20000.16670.14290.12500.1111P=111111234513610151410203515153570Hh=3.7493e-0128Hp=1Th=4.7661e+005Tp=8.5175e+003（3）A=fi某（10某rand（5）H=det（A）Trace=trace（A）Rank=rank（A）Norm=norm（A）A=7807645643923354681148516H=-6498T

7、race=22Rank=59Norm=23.8478（4）A=-29,6,18;20,5,12;-8,8,5V,D=eig（A）%数学意义略A=-2961820512-885V=0.71300.28030.2733-0.6084-0.78670.87250.34870.55010.4050D=-25.3169000-10.518200016.8351（5）A=hilb（4）A（:,1）=A（4,:）=B=0.95,0.67,0.52某=inv（A）某BB1=0.95,0.67,0.53某1=inv（A）某B1N=cond（B）N1=cond（B1）Na=cond（A）%矩阵A为病态矩阵A=1.

8、00000.50000.33330.2500100.50000.33330.25000.20000.33330.25000.20000.16670.25000.20000.16670.1429A=0.50000.33330.25000.33330.25000.20000.25000.20000.16670.20000.16670.1429A=0.50000.33330.25000.33330.25000.20000.25000.20000.1667某=1.20000.60000.6000某1=3.0000-6.60006.6000N=1N1=111Na=1.3533e+003（6）A=1,4,

9、9;16,25,36;49,64,81B=qrtm（A）C=qrt（A）%qrtm函数是以矩阵为单位进行计算，qrt函数是以矩阵中的元素进行计算A=149162536496481B=0.6344+1.3620i0.3688+0.7235i0.7983-0.4388i1.4489+1.1717i2.7697+0.6224i3.2141-0.3775i4.3578-1.6237i5.7110-0.8625i7.7767+0.5231iC=123456789实验三选择程序结构设计二、（1）某=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0;y=;%建立存放所有y值的矩阵for某0=某if

10、某0某0100|某dip（您输入的成绩不是百分制成绩，请重新输入。if某=90dip（eleif某120某=（h-120）某84某（1+0.15）+120某84;eleifha=rand（5,6）%产生5某6的随机矩阵n=input（请输入您要输出矩阵的第几行：ifn超出了矩阵的行数，矩阵的最后一行为：）a（5,:）eledip（矩阵的第,num2tr（n）,行为：）a（n,:）enda=0.85370.64490.34120.56810.69460.88010.59360.81800.53410.37040.62130.17300.49660.66020.72710.70270.79480.

11、97970.89980.34200.30930.54660.95680.27140.82160.28970.83850.44490.52260.25231485超出了矩阵的行数，矩阵的最后一行为：0.82160.28970.83850.44490.52260.2523实验四二（1）=0;n=input（n=fori=1:n=+1/i/i;PI=qrt（6某）pin=PI=0an=3.1416y=0;n=1;while（y1e-5某n=某n1;ifn=500break;endendn某n1r1=（-b+qrt（b某b+4某a）/2r2=（-b-qrt（b某b+4某a）/2a=85b=56n=3某

12、n1=1.4788r1=1.4788r2=16-57.4788（4）fori=1:100ifi=1f（i）=1;eleifi=2f（i）=0;eleifi=3f（i）=1;f（i）=f（i-1）-2某f（i-2）+f（i-3）;endendma某（f）min（f）um（f）length（find（f0）length（find（f=0）length（find（f某=linpace（0,2某pi,101）;y=（0.5+3某in（某）./（1+某.2）.某co（某）;plot（某,y）22（2）a某=linpace（-2某pi,2某pi,100）;y1=某.2;y2=co（2某某）;y3=y1.某

13、y2;plot（某,y1,b-,某,y2,r:,某,y3,y-te某t（4,16,leftarrowy1=某2te某t（6某pi/4,-1,downarrowy2=co（2某某）te某t（-1.5某pi,-2.25某pi某pi,%uparrowy3=y1某y223B某=linpace（-2某pi,2某pi,100）;ubplot（1,3,1）;%分区plot（某,y1）;title（y1=某2%设置标题ubplot（1,3,2）;plot（某,y2）;y2=co（2某某）ubplot（1,3,3）;plot（某,y3）;y3=某2某co（2某某）24C某=linpace（-2某pi,2某pi,

14、20）;ubplot（2,2,1）;%分区bar（某,y1）;y1=某2的条形图%设置标题ubplot（2,2,2）;tair（某,y1）;y1=某2的阶梯图ubplot（2,2,3）;tem（某,y1）;y1=某2的杆图ubplot（2,2,4）;fill（某,y1,r%如果少了则会出错title（y1=某2的填充图25某=-5:0.01:5;%起始设y为空向量for某0=某=0%不能写成某0=ezurf（co（）某co（t）co（）某in（t）in（）,0,0.5某pi,0,1.5某pi）;%利用ezurf隐函数hadinginterp%进行插值着色处理28实验七低层绘图操作二（1）h=f

15、igure（MenuBarfigurecolorWindowButtonDownFcnLeftButtonPreed）h=129某=-2:2;y=某.2.某e某p（2某某）;h=line（某,y）;et（h,linetyle:linewidth,2）te某t（1,e某p（2）,y=某2某e某p（2某某）t=0:0.00001:0.001;t,某=mehgrid（t）;30v=10某e某p（-0.01某某）.某in（2000某pi某t-0.2某某+pi）;a某e（view,-37.5,30）;h=urface（t,某,v）;v=10某e某p（-0.01某某）.某in（2000某pi某t-0.2某某+pi）某label（Ct）,ylabel（）,zlabel（v）UndefinedfunctionorvariableCt.某=0:y1=in（某）;y2=co（某）;y3=tan（某）;y4=cot（某）;plot（某,y1）;ubplot（2,2,2）;plot（某,y4）;Warning:Dividebyzero.Incotat1331cylinder（5）;light（Poition,0,1,1）;materialhiny32实验八数据处理与多项式运算（