川大《Matlab基础及应用1480》17春在线作业2

《Matlab 基础及应用 1480》15 春在线作业 1一、单选题（共 12 道试题，共 48 分）1. matlab 中求矩阵各列元素和的命令为( )A. diffB. sumC. CumprodD. cumtrapz正确答案：B2. 在命令窗口中输入“2*pi” ，则输出数值为（）A. 5.2832B. 6.2832C. 7.2832D. 8.2832正确答案：B3. 已知 A 为m×n 的矩阵，要提取其中 i 行 j 列（0C. 1D. 2正确答案：C7. MATLAB 中的一维插值命令是（）A. Interp1B. polyfitC. convD. deconv正确答案：A8. Matlab 图形用户界面开发环境中可以建立开关按钮的是（）A. Push ButtonB. Toggle ButtonC. Radio ButtonD. checkbox正确答案：A9. MATLAB 中的一维插值命令是（）A. Interp1B. polyfitC. convD. deconv正确答案：A10. 在 MATLAB 的若干通用操作指令中，清除命令窗口的所有显示内容的是（）A. clearB. claC. clfD. clc正确答案：D11. 当 a=2.4,使用取整函数得出 3，则该取整函数名为（）A. fixB. roundC. ceilD. floor正确答案：C12. 下面正确的说法是（）A. 函数文件与 M 文件一样，都可直接运行B. 函数文件中的变量是局部变量C. 命令与函数是同一个概念D. 矩阵运算与数组运算相同。

正确答案：B《Matlab 基础及应用 1480》15 春在线作业 1二、多选题（共 12 道试题，共48 分。

）1. 在 matlab 中可以完成绘图功能的命令是（）。

第二章 MATLAB 应用基础 练习题1、下列变量名中____A _________是合法变量。

A 、pi,exe_01B 、x*y,x1C 、1a,ifD 、abs, b.m2、已知a=0:4, b=1:5, 下面的运算表达式出错的为_______D ______。

A 、a+bB 、a ./bC 、a ’*bD 、a*b3、将矩阵A=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡987654321用________D_____命令可以变为A=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡741852963。

A 、A ’ B 、det(A) C 、inv(A) D 、rot90(A)4、已知x 为一个向量，计算ln(x)的运算为______B______。

A 、ln(x)）B 、log （x ）C 、Ln(x)D 、log10(x)5、____A____产生均匀分布的随机矩阵，元素取值范围0.0～1.0。

A 、 rand(m,n)B 、eye(m,n)C 、magic(N)D 、randn(m,n)6、____D____产生正态分布的随机矩阵。

A 、 rand(m,n)B 、eye(m,n)C 、magic(N)D 、randn(m,n)7、已知x 为一个向量，计算其余弦函数的运算为______C ______。

A 、COS （X ）B 、COS （x ）C 、cos(x)D 、cosx8、用“from ：step ：to ”方式得到从0到4π步长为0.4π的变量；使用linspace 函数生成向量从0到2*pi 等分成100个点的变量；使用logspace 函数生成向量从1到100对数等分成3个点的变量。

9、求矩阵⎥⎦⎤⎢⎣⎡4321的转置矩阵、逆矩阵、矩阵的秩、矩阵的行列式值、矩阵的三次幂、矩阵的特征值和特征向量。

10.我国人口按照2000年第五次全国人口普查的结果为12.9533亿，如果年增长率为1.07%，求公元2010年末的人口数。

计算人口的公式为：P1=P0(1+r)^n ，其中：P1为几年后的人口，P0为人口初值，r 为年增长率，n 为年数。

《Matlab基础与应用》课程网上考试题库第一章（MATLAB概述）一、单项选择题(每题2分)1、MATLAB的含义是( )A：矩阵B：实验室C：矩阵实验室D：矩阵分析室答案：C2、下列课程中哪个MATLAB不是其基本教学工具的是（）A：线性代数B：数学分析C：数字信号处理D：时间序列分析答案：B3、下列选项中，不是MATLAB数值运算功能的是( )A：多项式与有理分式运算B：优化处理C：解代数方程D：数值积分答案：C4、下列不是MATLAB中专业工具箱的是( )A：基本工具箱B：控制系统工具箱C：最优工具箱D：金融工具箱答案：A5、MATLAB的命令窗口无法显示的是( )A：向量B：图形C：输入命令D：矩阵答案：B6、MATLAB的帮助命令包括( )A：help和lookfor B：help和find C：find D：clear 答案：A7、启动MATLAB文本编辑器的方法不包括( )A：菜单操作B：命令操作C：命令按钮操作D：文本操作答案：D8、在MATLAB的使用中，如果用户对积分函数quad不了解，可以在命令窗口输入如下哪个命令得到帮助信息( )A：find quad B：? quad C：help quad D：know quad9、MATLAB中的lookfor命令只对M文件的哪一行进行关键字搜索( )A：第一行B：第二行C：第三行D：第四行答案：A10、在MATLAB命令窗口中输入以下哪项命令时，无法进入帮助窗口( )A：helpwin B：help C：helpdesk D：doc答案：B11、利用（）命令设置搜索路径。

A：lookfor B：load C：whos D：path 答案：D12、在（）中可键入各种MATLAB的命令、函数和表达式，并显示除图形外的所有运算结果。

A：当前目录浏览器窗口B：历史命令窗口C：命令窗口D：工作空间浏览器窗口答案：C13、（）用于显示所有MATLAB工作空间中的变量名、数据结构、类型、大小和字节数。

MATLAB 作业二参考答案1、试求下面线性代数方程的解析解与数值解，并检验解的正确性。

2932114010110503848246303356684953X -----⎡⎤⎡⎤⎢⎥⎢⎥----⎢⎥⎢⎥=⎢⎥⎢⎥---⎢⎥⎢⎥------⎣⎦⎣⎦【求解】求出A , [A;B ] 两个矩阵的秩，可见二者相同，所以方程不是矛盾方程，应该有无穷多解。

>> A=[2,-9,3,-2,-1; 10,-1,10,5,0; 8,-2,-4,-6,3; -5,-6,-6,-8,-4]; B=[-1,-4,0; -3,-8,-4; 0,3,3; 9,-5,3]; [rank(A), rank([A B])] ans = 4 4用下面的语句可以求出方程的解析解，并可以验证该解没有误差。

>> x0=null(sym(A));x_analytical=sym(A)\B; syms a; x=a*[x0 x0 x0]+x_analyticalx =[ a+967/1535, a-943/1535, a-159/1535][ -1535/1524*a, -1535/1524*a, -1535/1524*a][ -3659/1524*a-1807/1535,-3659/1524*a-257/1535,-3659/1524*a-141/1535] [ 1321/508*a+759/1535, 1321/508*a-56/1535, 1321/508*a-628/1535] [ -170/127*a-694/307, -170/127*a+719/307, -170/127*a+103/307] >> A*x-B ans =[ 0, 0, 0] [ 0, 0, 0] [ 0, 0, 0] [ 0, 0, 0]用数值解方法也可以求出方程的解，但会存在误差，且与任意常数a 的值有关。

>> x0=null(A); x_numerical=A\B; syms a; x=a*[x0 x0 x0]+x_numerical; vpa(x,10) ans =[ .53*a+.36, .53*, .53*] [*a+.89,*a+.88e-1,*a+.88e-1] [ *a, *a, *a][ .13*, .13*,.13*] [*, *a+, *a+.24] >> A*x-B ans =[ 1/1984*a, 1/1984*a, 1/1984*a] [ -5/*a, -5/*a, -5/*a][ -25/1984*a, -25/1984*a, -25/1984*a][ 13/1984*a, 13/1984*a, 13/1984*a]2、求解下面的联立方程，并检验得出的高精度数值解（准解析解）的精度。

实验一 MATLAB 工作环境熟悉及简单命令的执行一、实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。

二、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单（Menu ）、工具条 （Toolbar ）、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。

三、实验步骤：1、启动MATLAB ，熟悉MATLAB 的桌面。

2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。

（1）（365-52⨯2-70）÷3（2）>>area=pi*2.5^2（3）已知x=3，y=4，在MATLAB 中求z ：()232y x y x z -= （4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。

m1=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡11514412679810115133216 执行以下命令>>m1( 2 , 3 )>>m1( 11 )>>m1( : , 3 )>>m1( 2 : 3 , 1 : 3 )>>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1)（5）执行命令>>help abs查看函数abs 的用法及用途，计算abs( 3 + 4i )（6）执行命令>>x=0:0.1:6*pi;>>y=5*sin(x);>>plot(x,y)（6）运行MATLAB 的演示程序，>>demo ，以便对MATLAB 有一个总体了解。

四、思考题1、以下变量名是否合法？为什么？（1）x2（2）3col（3）_row（4）for2、求以下变量的值，并在MATLAB 中验证。

matlab基础与应用教程课后答案【篇一：matlab教程基本应用练习题及解答】txt>要求：将每题的答案（命令行和运行结果、图片或m文件的文件名及具体内容）直接拷贝插入到各题的下方：（1）若为命令行，要求将提示符“”一起拷入，并在右侧用“％”注明命令行的每条命令的作用；（2）若为多个运行结果，拷入后要求解释每个结果具体对应题目中的哪个要求，也在右侧用“％”注明；（3）如果为m文件，除了将文件名和此文件的具体内容全部拷入外，再将所有原始m文件和本练习题电子版放在一个以“专业班级＋本人姓名＋学号”命名的文件夹内一起上传。

ans =1 0 00 1 00 0 1ans =0 0 00 0 0ans =1 1 11 1 11 1 11 1 120+40*rand(1,10) %生成10个在区间[20,60]上均匀分布的随机数。

ans =44.617351.677556.872549.528327.050636.228257.418856.676236 .410855.7460二．1）计算向量（2，4，6，8）的最大值，最小值，平均值，中值，排序，总和值； max([2,4,6,8])%最大值ans =8min([2,4,6,8])% 最小值ans =2mean([2,4,6,8])% 平均值ans =5median([2,4,6,8])% 中值ans =5sort([2,4,6,8])% 排序ans =2 4 6 8sum([2,4,6,8])% 总和值ans =202）在行向量（2，4，6，8）和（1，3，5，7）之间实施加减乘除及幂运算； a=[2 4 6 8];b=[1 3 5 7];a+b%加ans =3 7 11 15a-b%减ans =1 1 1 1a.*b%乘ans =2 12 30 56a./b%除ans =2.0000 1.3333 1.2000 1.1429a.^b%幂ans =2 64 7776 20971523）在向量（1，3，5）与标量2之间实施加减乘除及幂运算。

第1章一、思考题4.( 1) B=A(2:5,1:2:5)(2)A(7)=[](3)A=A+30(4)size(A)( 5) t(find(t==0))=eps(6)t=reshape(x,3,4)(7)abs('matlab')(8)char(93)5.A=[97,67,34,10;-78,75,65,5;32,5,-23,-59;0,-12,54,7] (1)B=A(1:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*C(2)A(find(A>=50 & A<=100)) 二、实验题=-74/180*pi;y=-27/180*pi;sin(x*x+y*y)/sqrt(tan(abs(x+y)))+pi一・・,exp*a).*sin(a+=[2,4;,5];log(x+sqrt(1+x.*x))/24.A=[3,54,2;34,-45,7;87,90,15];B=[1,-2,67;2,8,74;9,3,0];(1)A*Bans =129 432 41977 -407 - 1052402 591 12489A.*Bans =3 -108 13468 -360 518783 270 0(2) A A3ans =-28917 240246 -4368137883 -259101 27669171333 252504 38673A.A3ans =27 157464 839304 -91125 343658503 729000 3375( 3) A/Bans =B/A ans =（4）[A,B]ans =3 54 2 1 -2 6734 -45 7 2 8 7487 90 15 9 3 0[&[1,3],:)砂2]ans =3 54 287 90 15600 183 -81684 282 72615 6 825=1+2i;b=3+4i; c=exp(pi*i/6) c =+ c+a*b/(a+b) ans =+ 第2章一、思考题=0; for n=0:63s=s+2A n;end disp(s)n=0:63;s=sum(2.An)二、实验题1.x=input( ' 输入一个四位整数：' );y=[fix(x/1000),mod(fix(x/100),10),mod(fix(x/10),10),mod(x,1 0)] z=mod((y+7),10)x=z(3)*1000+z(4)*100+z(1)*10+z(2)2. gh=input( ' 输入工号' ); h=input( ' 输入工时' );dj=84;if h>120gz=dj*120+*dj*(h-120);elseif h<60gz=dj*h-700;else gz=dj*h;end format bank; display([gh,gz])3. 循环结构n=input( 'input n:' );s=0;for k=1:n s=s+1/k A2;end display(sqrt(s*6)) 向量运算n=input( 'input n:' ); k=1:n;display(sqrt(sum(1./k.A2)*6))4.y=0;k=0;while y<3k=k+1; y=y+1/(2*k-1);end display([k-1,y-1/(2*k-1)])5.x0=0;x=1;k=0;a=input( 'a=' ); b=input( 'b=' );while abs(x-x0)>=1e-5 && k<500 x0=x;x=a/(b+x0);k=k+1;end display([k,x]); display([(-b+sqrt(bA2+4*a))/2,(-b-sqrt(bA2+4*a))/2]);6.y=fun(40)/(fun(30)+fun(20))(1)函数文件function f=fun(n) f=n+log(nA2+5);(2)函数文件function f=fun(n) a=1:n;f=sum(a.*(a+1));第3章一、思考题4.t=0::;y=sqrt(3)/2*exp(-4*t).*sin(4*sqrt(3)*t+pi/3);5.x=-10::10; y=linspace(-6,6,size(x,2))z=x.A3+3*x.*y.A2;plot3(x,y,z)6.x=100:100:400;y=100:100:400;z=[636,697,624,478;698,712,630,478;680,674,598,412;662,626,552,334]; [X,Y]=meshgrid(x,y);mesh(X,Y,z)二、实验题1.( 1 ) x=-10::10;plot(x,x-x.A3/3/2)( 2) plot(x,exp(-x.*x/2)/2/pi)(3)x=-8::8;plot(x,sqrt((64-x.*x)/2))( 4) t=0:pi/100:6*pi;plot(t.*sin(t),t.*cos(t))2.( 1) x1=linspace(0,1,100);y1=2*;t=linspace(0,pi,100);x=sin(3*t).*cos(t);y=sin(3*t).*sin(t);plot(x1,y1,'r-',x,y,'b:');text,1,'y='); text,,'x=sin(3t)cos(t)');text,,'y=sin(3t)sin(t)');(2)subplot(1,2,1);scatter(x1,y1,10)title('y=');subplot(1,2,2);scatter(x,y,10)3.subplot(1,2,1);x=1:1:100; y=sin(1./x);plot(x,y) subplot(1,2,2);fplot('sin(1/x)',[1,100])4.subplot(2,2,1);bar(t,y);subplot(2,2,2);stairs(t,y);subplot(2,2,3);stem(t,y) subplot(2,2,4); semilogy(t,y);3.theta=linspace(-pi,pi,100); ro=5.*cos(theta)+4; polar(theta,ro); fi=linspace(0,2*pi,100);a=1 r=a.*(1+cos(fi));polar(fi,r);4.plot3(x,y,z);t=0::1;x=t;y=t.A2;z=t.A3;plot3(x,y,z);7.x=-30::0;y=0::30; [x,y]=meshgrid(x,y);z=10.*sin(sqrt(x.A2+y.A2))./sqrt(1+x.A2+y.A2); meshc(x,y,z);8. x=linspace(-3,3,100);y=linspace(-3,3,100);[x y]=meshgrid(x,y);fxy=-5./(1+x.A2+y.A2);i=find(abs(x)<= & abs(y)<=; fxy(i)=NaN;surf(x,y,fxy)9.u=linspace(1,10,100);v=linspace(-pi,pi,100); [u,v]=meshgrid(u,v);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4*u.A2;x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.A2; surf(x,y,z);shading interp;light('position',[1,0,1]);10. t=0:pi/100:2*pi; y=sin(t);comet(t,y) 第4章一、思考题5.(1)A=eye(3);(2)C=100+(200-100)*rand(5,6);( 3) D=1+sqrt*randn(1,500);( 4) E=ones(size(A));( 5) A=A+30*eye(size(A));( 6) B=diag(diag(A))二、实验题1.P=pascal(5);H=hilb(5); Dp=det(P);Dh=det(H); Kp=cond(P);Kh=cond(H);P矩阵的性能更好，因为Kp较小2.A=[1,-1,2,3;0,9,3,3;7,-5,0,2;23,6,8,3] B=[3,pi/2,45;32,-76,sqrt(37);5,72,;exp(2),0,97] A1=diag(A);B1=diag(B);A2=triu(A);B2=triu(B); A3=tril(A);B3=tril(B); rA=rank(A);rB=rank(B);nA=norm(A);nb=norm(B); cA=cond(A);cB=cond(B);3.A=[31,1,0;-4,-1,0;4,-8,-2] ；[V,D]=eig(A);4.A=diag([-1,-1,-1,-1],-1)+diag([-1,-1,-1,-1],1)+diag([2,2,2,2,2])b=[1,0,0,0,0]';x1=inv(A)*b; x2=A\b; [L,U]=lu(A); x3=U\(L\b); [Q,R]=qr(a); [Q,R]=qr(A); x4=R\(Q\b) R=chol(A); x5=R\(R'\b) 5.B=sparse(A); x1=inv(B)*b; x2=B\b; [L,U]=lu(B); x3=U\(L\b); 第5章一、思考题3.A=randn(10,5); mean(A) std(A) max(max(A)) min(min(A)) sum(A,2) sum(sum(A)) sort(A,1) sort(A,2,'descend') 二、实验题1.A=rand(1,30000); mean(A) std(A) max(A) min(A) size(find(A>)/size(A)2.h=[466,715,950,1422,1635]; w=[,,,,];hh=[500,900,1500]; ww=interp1(h,w,hh,'spline')3.x=linspace(1,10,50); y=log(x);第6章一、思考题2.fx=i nlin e('1./(1+x.A 2)');[I,n]=quad(fx,-100000,100000,1e-3); [I,n]=quadl(fx,-100000,100000,1e-3); x=-100000::100000; y=1./(1+x.*x);f=polyfit(x,y,5); yy=polyval(f,x); plot(x,y,'r-',x,yy,'g.') 4. N=64; T=5;t=linspace(0,T,N); x=3*exp(-t); % dt=t(2)-t(1); f=1/dt; X=fft(x); F=X(1:N/2+1); f=f*(0:N/2)/N;plot(f,abs(F),'-*') % %% 求各采样点样本值 x% % %采样点数 采样时间终点给岀N 个采样时间ti(l=1:N) 采样周期采样频率 (Hz)计算 x 的快速傅立叶变换 X% F(k)=X(k)(k=1:N/2+1)% 使频率轴 f 从零开始 % 绘制振幅 - 频率图xlabel('Frequency');ylabel('|F(k)|')5.(1)p1=[1 2 0 0 7];p2=[1 -2];p3=[1 0 0 5 1];p12=conv(p1,p2);p=p12+[zeros(1,size(p12,2)-size(p3,2)),p3]; roots(p)(2)A=[-1,4,3;2,1,5;0,5,6];Px=polyval(p,A) Pxm=polyvalm(p,A)6.(1) z=fzero('3*x-sin(x)+1',0)(2)建立函数文件 function F=myfun(X)x=X(1);y=X(2); F(1)=x*x+y*y-9;F(2)=x+y-1; 在命令窗口中输入以下命令： x=fsolve(@myfun,[3,0]',optimset('Display','of f'))trapz(x,y);3.(1)fx=inline('-2*y+2*x*x+2*x');[t,y]=ode23(fx,[0,],1)(2)fx=inline('y-exp(x)*cos(x)');[t,y]=ode23(fx,[0,3],1)二、实验题1.for x=1:3fx=[x,xA2,xA3;1,2*x,3*x;0,2,6*x]; diff(fx)end2.(1 ) x=0::1;y=x.A10+10.Ax+1./log10(x);dy=diff(y)/;(2)x=0::1;y=log(1+x);dy=diff(y,2)/;plot(x(1:99),dy)3.(1 ) fx=inline( 'x.A2.*sqrt(2*x.*x+3)');quad(fx,1,5)(2)fx=inline( 'x./sin(x).A2' ); quad(fx,pi/4,pi/3)(3)fx=inline( 'abs(cos(x+y))' ); dblquad(fx,0,pi,0,pi)(4)syms x y;fx=x*y;int(int(fx,yA2,y+2),-1,2)x的积分区间为【0, 2】时fx=inline( 'x.*y' );dblquad(fx,0,2,-1,2)4.x=::;y=[,,,,,,];trapz(x,y)5.(1)yp=i nlin e( '-+si n( 10*x))*y' );[t,y]=ode23(yp,[0,5],1);(2)令x1y, x2y ,x3 y''，则可写出原方程的状态方程形式: x1x2X2 X3cos 1x3cost 5 X3 X2 X1(t 1)2 3 sin t0 1 0 X1 0=>> x20 0 1 X2 01 / 5cos2t X3 12X3 cost3 sint (t 1)2建立函数文件fun cti on y=ztfu n( t,x)b=[0;0;cos(t)];y=[0,1,0;0,0,1;-1/(3+si n( t)),-1,5*cos(2*t)/(t+1)A2]*x+b; 解微分方程组[t,y]=ode23(@ztfu n,[0,5],[1;0;2]);6.建立函数文件fun cti on yy=ztfu n( t,y)yy=[ y( 2)*y(3);-y(1)* y(3) ;*y(1)*y(2)];解微分方程组[t,y]=ode23(@ztfu n,[0,5],[0;1;1])第7章一、思考题3.(1 )数值积分fx=i nlin e('exp(x).*(1+exp(x)).A2'); quad(fx,0,log(2)) 符号积分f=sym('exp(x)*(1+exp(x))A2');v=in t(f,0,log (2));eval(v)(2 )略二、实验题1.A=sym('[1,2,3;x,y, z; 3,2,1]')rank(A)inv(A)det(A)2.(1)y=sym('sqrt(x+sqrt(x+sqrt(x)))');y1=diff(y)y2=diff(y,'x',2)(2)syms x y;fxy=si n( x A2*y)*exp(-x A2-y);diff(diff(fxy,x),y)3.(1)syms xin t(1/(1+xA4))(2)syms x tin t((-2*x*x+1)/(2*x*x-3*x+1)A2,x,cos(t),exp(2*t))4.syms n xsymsum(1/(2* n+1)/(2*x+1)A(2* n+1), n,0,i nf)symsum(1/(2* n+1)/(2*x+1)A(2* n+1), n, 0,5)5.(1)syms xtaylor((exp(x)+exp(-x))/2,5,0)(2)syms a xtaylor(exp(-5*x)*si n(3*x+pi/3),5,a)6.(1)x=solve(sym('xA3+a*x+1=0'))(2)[x y]=solve(sym('sqrt(xA2+yA2)-100=0,3*x+5*y-8=0'))y' y17.方程转化为：’，严y1' y 1 一符号解[y1,y11]=dsolve(,Dy=y1,Dy1+y=1-t A2/pi,,,y(-2)=5,y1(-2)=5',,t,)数值解编写函数文件fun cti on yy=ztfu n( t,y)yy=[y(2);1-tA2/pi-y(1)];在命令窗口输入以下命令[t,y]=ode45(@ztfu n,[-2,7],[-5;5]);t=li nspace(-2,7,49)y2=y8.[x,y]=dsolve('Dx=3*x+4*y,Dy=-4*x+3*y','x(0)=0,y(0)=1')第9章二、实验题1.(1 )新建一个Blank GUI。

matlab习题2答案MATLAB习题2答案在MATLAB习题2中，我们学习了如何使用MATLAB编程语言解决各种数学问题和数据分析任务。

在本篇文章中，我们将探讨一些常见的习题，并给出相应的答案和解释。

第一道习题是关于矩阵运算的。

我们需要编写一个MATLAB程序，计算两个矩阵的乘积。

在MATLAB中，可以使用矩阵乘法运算符“*”来实现这个功能。

例如，如果我们有两个矩阵A和B，我们可以使用以下代码来计算它们的乘积：```matlabC = A * B;```接下来，我们需要解决一个关于插值的问题。

我们需要编写一个程序，使用插值方法来估算给定点的函数值。

在MATLAB中，可以使用interp1函数来实现这个功能。

例如，如果我们有一个已知的函数y和一组离散的点x，我们可以使用以下代码来进行插值计算：```matlabx_interp = 0:0.1:10;y_interp = interp1(x, y, x_interp, 'spline');```最后，我们需要解决一个关于数据可视化的问题。

我们需要编写一个程序，绘制给定数据的散点图和拟合曲线。

在MATLAB中，可以使用scatter和plot函数来实现这个功能。

例如，如果我们有一组数据点x和y，我们可以使用以下代码来进行数据可视化：```matlabscatter(x, y);hold on;p = polyfit(x, y, 1);y_fit = polyval(p, x);plot(x, y_fit, 'r');```通过以上习题的练习，我们可以更好地掌握MATLAB编程语言的基本语法和常用函数。

希望本篇文章对大家有所帮助，也希望大家能够在日常工作和学习中充分利用MATLAB的强大功能，提高工作效率和解决问题的能力。

作业一4、写出完成下列操作的命令.(1）将矩阵A第2～5行中第1,3,5列元素赋给矩阵B。

>> M=［0：1：48]；〉〉A=reshape（M，7,7）A =0 7 14 21 28 35 421 8 15 22 29 36 432 9 16 23 30 37 443 10 17 24 31 38 454 11 18 25 32 39 465 12 19 26 33 40 476 13 20 27 34 41 48>> B=A（2:5,1：2：5）B =1 15 292 16 303 17 314 18 32（2）删除矩阵A的第七行元素。

〉〉A（7，:）=［］A =0 7 14 21 28 35 421 8 15 22 29 36 432 9 16 23 30 37 443 10 17 24 31 38 454 11 18 25 32 39 465 12 19 26 33 40 47 (3)将矩阵A的每个元素值加30.>〉A=A+30A =30 37 44 51 58 65 7231 38 45 52 59 66 7332 39 46 53 60 67 7433 40 47 54 61 68 7534 41 48 55 62 69 7635 42 49 56 63 70 77(4求矩阵A的大小和维素。

sizeA = size（A)dA = ndims（A）sizeA =6 7dA =2（5）将向量t的0元素用机器零来代替。

〉> t=[1 2 3 4 0 5]；t =1 2 3 4 0 5〉〉find（t==0）ans =5〉> t(5)=epst =1.0000 2。

0000 3。

0000 4。

0000 0.0000 5。

0000（6）将含有12个元素的向量x转换成3＊4矩阵.〉> x=［0：11]x =0 1 2 3 4 5 6 7 8 9 10 11 〉〉y=reshape（x，3，4）y =0 3 6 91 4 7 102 5 8 11（7）求一个字符串的ASCII。