还剩17页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
第二章
3.设矩阵A为A=[24239216;6574241121;345987521;8424253121;4321456421];
(1)B=A2:51:2:5B=652421349821842121434521
(2)A7=[]A=246534843235422192498424521117553646212112121
(3)A+30
(4)sizeA;ndimsA
(5)题目有误
(6)reshapex34
(7)absx
(8)charx
4.L1=000010000L2=111110000L3=000111000L4=
4565.1B=A1:3:C=A:1:2D=A2:43:4E=B*CB=
23.
000010.0000-
0.
7780041.0000-
45.
000065.
00005.
000032.
00005.
0000032.0000C=
23.
000010.
000041.0000-
45.
000032.
00005.
00006.0000-
9.5400D=
65.
00005.
0000032.
000054.
00003.1400E=
1.0e+003*
0.9141-
0.
22391.
20802.
71231.1330-
0.21032EDans=010001EDans=110111E|Dans=111111~D|~Eans=001000findA=10A25ans=
156.allAans=0anyAans=1isnanAans=0100000isinfAans=0011000isfiniteAans=
10001117.A
1.x1=’学号’;A
1.x2=’姓名’;A
1.x3=’专业’;A
1.x
4.x41=’成绩1’;……….A
2.x1=’学号’;A
2.x2=’姓名’;A
2.x3=’专业’;A
2.x
4.x41=’成绩1’;……….A
3.x1=’学号’;A
3.x2=’姓名’;A
3.x3=’专业’;A
3.x
4.x41=’成绩1’;……….A
4.x1=’学号’;A
4.x2=’姓名’;A
4.x3=’专业’;A
4.x
4.x41=’成绩1’;……….A
5.x1=’学号’;A
5.x2=’姓名’;A
5.x3=’专业’;A
5.x
4.x41=’成绩1’;……….
8.1sizeBans=22ndimsBans=22B2ans=[3x3doubleB4ans={3x3cell}3B3=[]B=
[1][3x3double]{3x3cell}B{3}=[]B=
[1][3x3double][]第三章
1.1A=eye32A=100+100*rand563A=1+sqrt
0.2*randn10504B=onessizeA5A+30*eyesizeA6B=diagdiagA
2.B=rot90AC=rot90A-
13.B=invA;A的逆矩阵C=detA;A的行列式的值D=A*BE=B*AD=E因此A与A-1是互逆的
4.A=[42-1;3-12;1230];b=[2;10;8];x=invA*bx=-
6.
000026.
666727.
33335.1diagA;主对角线元素ans=1159triuA;上三角阵ans=1-12301-4200520009trilA;下三角阵ans=100051003050111509rankA;秩ans=4normA;范数ans=
21.3005condA;条件数ans=
11.1739traceA;迹ans=16
(2)略
6.A=[
110.5;
110.25;
0.
50.252]A=
1.
00001.
00000.
50001.
00001.
00000.
25000.
50000.
25002.0000[VD]=eigAV=
0.
72120.
44430.5315-
0.
68630.
56210.4615-
0.0937-
0.
69760.7103D=-
0.
01660001.
48010002.5365第四章
1.a=input请输入一个4位数:;whilea1000|a9999a=input输入错误,请重新输入一个4位数:;endb=fixa/1000;c=remfixa/10010;d=remfixa/1010;e=rema10;b=b+7;c=c+7;d=d+7;e=e+7;b=remb10;c=remc10;d=remc10;e=reme10;g=b;b=d;d=g;g=c;c=e;e=g;a=1000*d+100*e+10*b+c;disp[加密后num2stra]
2.a=input请输入a:;b=input请输入b:;c=input请输入c:;x=
0.5:1:
5.5;x1=x=
0.5x
1.5;x2=x=
1.5x
3.5;x3=x=
3.5x=
5.5;y1=a.*x.^2+b.*x+c;y2=a*sinb^c+x;y3=logabsb+c./x;y=y
1.*x1+y
1.*x2+y
3.*x3;dispy
3.x=fixrand120*89+10;x1=fixsumx/20;disp[平均数是num2strx1]m=remx2==0xx1;n=findm;disp[小于平均数的数是num2strxn];
4.A=input请输入20个数的一个行向量;m=A;a=m;b=m;form=Aifa=ma=m;elseifb=mb=m;endenddisp[最小数是num2stra]disp[最大数是num2strb]
5.s=0;a=0;forb=1:64c=2^a;a=a+1;s=s+c;enddisp[2的0次方到63次方的和是num2strs]6sum1=0;forn=1:100x=-1^n+1*1/n;sum1=sum1+x;enddisp[当n取100时sum=num2strsum1]sum2=0;forn=1:1000x=-1^n+1*1/n;sum2=sum2+x;enddisp[当n取1000时sum=num2strsum2]sum3=0;forn=1:10000x=-1^n+1*1/n;sum3=sum3+x;enddisp[当n取10000时sum=num2strsum3]
8.clearalla=input请输入一个矩阵;b=input请再输入一个矩阵注意两矩阵要可以相乘;[f1f2]=juzhenjiab;
9.clearally=qiuhe1001+qiuhe502+qiuhe10-1;dispy
10.1S=1082x=41220y=246第五章
1.1x=-10:
0.1:10;y=
100./1+x.^2;plotxy2x=-10:
0.1:10;y=1/2*pi*exp-x.^2/2;plotxy3ezplotx^2+y^2=14t=-10:
0.1:10;x=t.^2;y=5*t.^3;plotxy
2.1theta=0:
0.01:2*pi;rho=5*costheta+4;polarthetarho2theta=
0.001:
0.1:2*pi;rho=
12./sqrttheta;polarthetarho3theta=
0.001:
0.1:2*pi;rho=
5./costheta-7;polarthetarho4theta=
0.001:
0.1:2*pi;rho=pi/
3.*theta.^2;polarthetarho
3.1t=0:pi/100:2*pi;x=cost;y=sint;z=t;plot3xyz2u=0:pi/100:2*pi;v=0:pi/100:2*pi;x=1+cosu.*cosv;y=1+cosu.*sinv;z=sinu;plot3xyz
345.plot函数x=linspace-1010200;y=[];forx0=xifx00y=[yx
0.^2+1+x
0.^1/4+5];elseifx0==0y=[y0];elseifx00y=[yx
0.^3+sqrt1-x0-5];endendplotxyfplot函数fplotx
0.*x.^3+sqrt1-x-5+x==
0.*0+x
0.*x.^2+1+x.^1/4+5[-1010]第六章
1.A=randn105
(1)meanA;均值stdA;标准方差
(2)maxmaxA;最大元素minminA;最小元素
(3)B=sumA2;A每行元素的和sumB;A全部元素之和
(4)sortA;A的每列元素按升序排列sortA2’descend’;A的每行元素按将序排列
2.
(1)
(2)X=
[149162536496481100];Y=1:10;X1=1:100;Y1=interp1XYX1cubic
3.x=
[165123150123141];y=
[187126172125148];P=polyfitxy3P=
1.0e+003*-
0.
00000.0013-
0.
17798.4330所以它的线性拟合曲线为px=
1.3x2—
177.9x+
84334.1P1=
[032];P2=[5-12];P3=[10-
0.5];P=convconvP1P2P3P=
015.
00007.0000-
3.
50000.5000-
2.0000-
2.0000所以Px=15x5+7x4-
3.5x3+
0.5x2-2x-22rootsPans=
0.
70710.1000+
0.6245i
0.1000-
0.6245i-
0.7071-
0.66673i=0:10;xi=
0.2*i;polyvalPxians=-
2.0000-
2.3920-
2.6112-
1.
70242.
710415.
000042.
112094.
1408184.
9056332.
5264560.
00005.1建立函数文件functionf=fxyux=u1;y=u2;f=
3.*x.^2+2*x.*y+y.^2在命令窗口中输入以下命令[Ufmin]=fminsearchfxy
[11]结果U=
1.0e-004*-
0.
06750.1715fmin=
1.9920e-010
(2)f=inline-sinx-cosx.^2;fmax=fminbndf0pifmax=
0.
73106.
(1)x=[pi/6pi/4pi/3];f=inlinesinx.^2+cosx.^2;dx=difff[x5*pi/12]/pi/12可参见第157页例题
6.19dx=000x=pi/2时单独计算x=pi/2;f=inlinesinx.^2+cosx.^2;difff[xpi]/pi/2ans=0
(2)x=1:3;f=inlinesqrtx.^2+1;dx=difff[x4]结果dx=
0.
82190.
92620.
96087.
(1)f=inlinesinx.^
5.*sin5*x;quadf0pians=
0.0982
(2)f=inline1+x.^
2./1+x.^4;quadf-11ans=
2.2214
(3)f=inlinex.*sinx./1+cosx.^2;quadf0pians=
2.4674
(4)f=inlineabscosx+y;dblquadf0pi0pians=
6.
28328.N=64;%采样点数T=5;%采样时间终点t=linspace0TN;%给出N个采样时间ti(i=1:N)y=exp-t;%求各采样点样本值ydt=t2-t1;%采样周期f=1/dt;%采样频率Y=ffty;%计算y的快速傅里叶变换YF=Y1:N/2+1;%Fk=Ykf=f*0:N/2/N;%使频率轴f从0开始plotfabsF%绘制振幅-频率图
9.
(1)矩阵求逆法A=[235;374;1-71];b=[10;3;5];x=invA*bx=-
1.8060-
0.
53733.0448矩阵除法法A=[235;374;1-71];b=[10;3;5];x=A\bx=-
1.8060-
0.
53733.0448矩阵分解法A=[235;374;1-71];b=[10;3;5];[LU]=luA;x=U\L\bx=-
1.8060-
0.
53733.0448
(2)方法同
(1)
10.函数文件line_solutionAbfunction[xy]=line_solutionAb[mn]=sizeA;y=[];ifnormb0ifrankA==rank[Ab]ifrankA==ndisp原方程组有唯一解x;x=A\b;elsedisp原方程组有无穷个解,特解为x,齐次方程组的基础解系为y;x=A\b;y=nullAr;enddisp方程组无解;x=[];endelsedisp原方程组有零解x;x=zerosn1;ifrankAndisp方程组有无穷个解,基础解系为y;y=nullAr;endend程序A=[21-11;42-21;21-1-1];b=[1;2;1];[xy]=line_solutionAb结果原方程组有无穷个解,特解为x,齐次方程组的基础解系为yWarning:Rankdeficientrank=2tol=
4.3512e-
015.Inline_solutionat11方程组无解x=[]y=-
0.
50000.
50001.
0000001.
00000011.
(1)f=inlinex-sinx./x;x=fzerof
0.5x=
0.87672f=inlinesinx.^
2.*exp-
0.
1.*x-
0.
5.*absx;x=fzerof
1.5x=
1.
673812.函数文件functionf=fxyux=u1y=u2f1=x-
0.6*sinx-
0.3*cosyf2=y-
0.6*cosx+
0.3*siny在命令窗口输入以下命令x=fsolvefxy[
0.
50.5]optimsetDisplayoff结果x=
0.
63540.
373415.A=[-120;-12-1;-12-1;-12-1;02-1];d=[-1;0;1];B=spdiagsAd55;b=
[10000];x=invB*bx=
0.
83330.
66670.
50000.
33330.1667。