《MATLAB语言程序设计》练习题

第1页《MATLAB语言程序设计》练习题1一、填空题1．设有程序A=[1,2,3,4;6,-5,6,3];B=[2,4,5,7];A1=A(:,[134]);A2=A(1,[2,3]).将下列命令的运行结果填在横线上size(A)=;A(1,:).*B=;length(B)=;A1=;A2=;A+3=;B.^2=2.A=[1,2,3,4;6,-5,6,3];B=[-2,4,-6,7,9];A1=(A3);A2=all(A(:,1)5);A3=any(A(2,:)10);A4=find(abs(B)3|abs(B)7);A5=find(abs(B)7&abs(B)3)A1=;A2=;A3=;A4=;A5=.3.n=1;a=1;whilesum(1:n)10n=n+1;a=a*n;endaa=.5.a=[90,70,65,54,76];n=length(a);b=zeros(1,n);fori=1:nifa(i)75b(i)=100-a(i);elseifa(i)=75b(i)=a(i)-60;endendbb=4.a=[];fori=1:3;forj=i:3;a(i,j)=i+j;end;endaa=.6.functiony=temf(a)iflength(a)2y=a;elsey=max(a)-min(a)endtemf([23456789])=;temf(78)=.第2页;).4(;2(3).);BAR(AR2).A(|;)(|).1(,214154101,321214153111VUAEBABAXABDBA：、TT与特征向量的特征值求解矩阵方程的秩设题试编写程序计算下列各二3512332232.3123322352:.2221432:.1321321321321321432143214321432143214321xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx、超定方程组恰定方程组欠定方程组程组的程序写出求解下列各线性方三.)(,cos3)sin()sin(cos2sin)(.3;)(,,1sin)ln(.2;02.1:,4321022223变换的计算并计算生成符号矩阵的根求方程下列各题运用符号运算功能计算四LaplacetAwtewttttettttettAdxxAdxAdxxesxxAxx、utttx的解函数求解微分方程用1)0(')0(,sin'2.4xxtxxxdsolve五、下列这组数据是美国1900—2000年人口的近似值（单位：百万）。时间t19001910192019301940195019601970198019902000人口y7692106123132151179203227250281（1）若.2cbtatyty的经验公式为与试编写程序计算出上式中的a、b、c;第3页（2）若.btaeyty的经验公式为与试编写程序计算出上式中的a、b;（3）画出数表中的散点图（红色点）,cbxaxy2中拟合曲线图(蓝色实心线),以及.btaey(绿色点划线_.).（4）图形标注要求:无网格线,横标注“时间t”，纵标注“人口数(百万)”，图形标题“美国1900—2000年的人口数据”。).10,5(]4,0[)2()),(:(;)1(,cos141)(12nnxfourieryfunctionntnxf、nk分别取内图形用子图画出该函数在程序的第一句为返回该函数并保存试编写函数式文件程序设六.2,,,%),,,(:)2(])()(4)([6)(:122122个数为区间的等分小区间的为积分的上下限的数值为函数名字符串程序的第一句为等分此公式中的区间被分为数式文件公式编写计算积分的函试用七nbafunnamenbafunnamesimpsonIfunctionn，xfxfxfnb-aISimpson、nkkkk.755)5(;755).4(;).3(;).2(].8,12,10,9,11[],,,,[],75,65,70,71,68[],,,,[.)5,4,3,2,1(90,）,（1000590).1(:54321543212分的人数于门课程平均成绩大于等统计出分的人数门课程成绩都大于等于统计出及格的人数分别统计出各门课程不均成绩成绩以及各名同学的平计算出各门课程的平均并假设的成绩课程课程课程课程课程名同学的五门课程矩阵理解为可将此分布之间的整数并服从正态阶矩阵使其元素在试产生编写程序解决下列问题八iiiN、第4页《MATLAB语言程序设计》练习题1参考答案参考程序：1.[2,4];[2,8,15,28];4;A1=[134;663];A2=[23];[4567;9-296;4162549];2.A1=[1100;0100];A2=0;A3=1;A4=[15];A5=[23];3.24;4.[234;045;006];4.[10105-624];6.66;78.二、参考程序：三、参考程序四、参考程序：1．solve('x^3-x^2+2=0')2．symsxsA=[x*sin(x),log(x);exp(-x),1/(s+x)]diff(A)int(A)A*A3．limit((cos(x)-cos(x)^(1/2))/x/(exp(sin(x))-1))4．symsxyz;F=[x^2+y^2+4*z^4;2*x*cos(y)+3*z^2;z^3*exp(x*y)]jacobian(F,[x,y,z])5．dsolve('D2x=-(2/t)*D1x-x+sin(t)','x(0)=1','Dx(0)=1')五、参考程序：t=[1900:10:2000];y=[7692106123132151179203227250281];p=polyfit(t,y,2)yp=p(1)*t.^2+p(2)*t+p(3)y1=log(y);A=[ones(11,1),t']B=A\y1';a=exp(B(1))b=B(2)ye=a*exp(b*t);plot(t,y,'r.',t,yp,'b-',t,ye,'g-.')第5页gridonxlabel('时间t')ylabel('人口数(百万)')title('美国1900—2000年的人口数据')六、参考程序：(1).functiony=fourier(x,n)y=0;fori=1:ny=y+cos(i*x)/(4*i^2-1);endy(2).x=0:pi/20:4*pi;y1=fourier(x,5);y2=fourier(x,10);subplot(211),plot(x,y1)subplot(212),plot(x,y2)七、参考程序functionI=simpson(fun,a,b,n)ifnargin4n=100;endh=(b-a)/2/n;s=0;fork=0:nx0=a+2*k*h;x1=a+(2*k+1)*h;x2=a+(2*k+2)*h;s=s+feval(fun,x0)+4*feval(fun,x1)+feval(fun,x2);endI=s*h/3;八、参考程序：mu=[68,71,70,65,75];sig=[11,9,10,12,9]cj=[];fori=1:5cj1=round(randn(90,1)*sig(i)+mu(i))cj=[cj,cj1];第6页endfori=1:90forj=1:5ifcj(i,j)100cj(i,j)=100;elseifcj(i,j)0cj(i,j)=0;endendendkn60=sum((cj60))%各门课程不及格的人数kmean=round(mean(cj));%各门课程的平均成绩std(cj);%各门课程的成绩的标准差studentmean=round(mean(cj,2));nn75=sum(sum((cj75),2)==5);k75=sum((mean(cj,2)75));disp('各门课程不及格的人数分别为:')disp('课程1课程2课程3课程4课程5')disp(sprintf('%4d%8d%8d%8d%8d',kn60(1),kn60(2),kn60(3),kn60(4),kn60(5)))disp('5门课程都大于75分的人数为:')disp(sprintf('%4d',nn75))disp('平均成绩大于75分的人数为:')disp(sprintf('%4d',k75))disp('课程1课程2课程3课程4课程5')disp(sprintf('%4d%8d%8d%8d%8d',kmean(1),kmean(2),kmean(3),kmean(4),kmean(5)))hist(studentmean)%全班同学平均成绩的分布图