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

《MATLAB语言程序设计》练习题1 填空题 设有程序 A=[1,2,3,4;6,-5,6,3];B=[2,4,5,7];A1=A(:,[1 3 4]);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= . 五、下列这组数据是美国1900—2000年人口的近似值（单位：百万）。 时间t 1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 2000 人口y 76 92 106 123 132 151 179 203 227 250 281 若试编写程序计算出上式中的a、b、c; 若试编写程序计算出上式中的a、b; 画出数表中的散点图（红色点）, 中拟合曲线图(蓝色实心线),以及 (绿色点划线_.). 图形标注要求:无网格线,横标注“时间t”，纵标注“人口数(百万)”，图形标题“美国1900—2000年的人口数据”。 《MATLAB语言程序设计》练习题1参考答案 参考程序： [2,4]; [2, 8, 15, 28]; 4;A1 =[1 3 4; 6 6 3];A2 =[ 2 3]; [4 5 6 7; 9 -2 9 6; 4 16 25 49]; 2. A1 =[ 1 1 0 0; 0 1 0 0];A2 = 0;A3 = 1;A4 =[1 5];A5 =[2 3]; 24; 4. [ 2 3 4;0 4 5;0 0 6]; 4. [10 10 5 -6 24]; 6. 66; 78. 二、参考程序： 三、参考程序 四、参考程序： solve(x^3-x^2+2=0) 2．syms x s A=[x*sin(x),log(x);exp(-x),1/(s+x)] diff(A) int(A) A*A 3．limit((cos(x)-cos(x)^(1/2))/x/(exp(sin(x))-1)) 4．syms x y z; 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=[76 92 106 123 132 151 179 203 227 250 281]; 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-.) grid on xlabel(时间t) ylabel(人口数(百万)) title(美国1900—2000年的人口数据) 六、参考程序： (1).function y=fourier(x,n) y=0; for i=1:n y=y+cos(i*x)/(4*i^2-1); end y (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) 七、参考程序 function I=simpson(fun,a,b,n) if nargin4 n=100;end h=(b-a)/2/n; s=0; for k=0:n x0=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); end I=s*h/3; 八、