matlab基础与应用教程课后答案
Matlab编程与应用习题和一些参考答案
Matlab编程与应用习题和一些参考答案Matlab 上机实验一、二3.求下列联立方程的解⎪⎪⎩⎪⎪⎨⎧=+-+-=-+=++-=--+41025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10];>> b=[4;4;9;4];>> c=a\b4．设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=81272956313841A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=793183262345B ，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。
>> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8];>> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7];>> C1=A*B'>> C2=A'*B>> C3=A.*B>> inv(C1)>> inv(C2)>> inv(C3)5．设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。
>> x=linspace(0,2*pi,101);>> y=cos(x)*(0.5+(1+x.^2)\3*sin(x));>> plot(x,y,'r')6．产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。
并求该矩阵全体数的平均值和均方差。
（mean var ）a=randn(8,6)mean(a)var(a)k=mean(a)k1=mean(k)i=ones(8,6)i1=i*k1i2=a-i1i3=i2.*i2g=mean(i3)g2=mean(g)10.利用帮助查找limit 函数的用法，并自己编写，验证几个函数极限的例子。
MATLAB语言基础与应用(第二版)第5章 习题答案
第5章习题与答案5.1用矩阵三角分解方法解方程组123123123214453186920x x x x x x x x x +-=⎧⎪-+=⎨⎪+-=⎩ 解答：>>A=[2 1 -1;4 -1 3;6 9 -1] A =2 1 -1 4 -13 6 9 -1 >>b=[14 18 20]; b =14 18 20 >> [L, U, P]=lu(A) L =1.0000 0 0 0.6667 1.0000 0 0.3333 0.2857 1.0000 U =6.0000 9.0000 -1.0000 0 -7.0000 3.6667 0 0 -1.7143 P =0 0 1 0 1 0 1 0 0 >> y=backsub(L,P*b’) y =20.0000 4.6667 6.0000 >> x=backsub(U,y) x =6.5000 -2.5000 -3.5000 5.2 Cholesky 分解方法解方程组123121332352233127x x x x x x x ++=⎧⎪+=⎨⎪+=⎩ 解答：>> A=[3 2 3;2 2 0;3 0 12] A =3 2 32 2 03 0 12>> b=[5;3;7]b =537>> L=chol(A)L =1.7321 1.1547 1.73210 0.8165 -2.44950 0 1.7321>> y=backsub(L,b)y =-11.6871 15.7986 4.0415>> x=backsub(L',y)x =-6.7475 28.8917 49.93995.3解答:观察数据点图形>> x=0:0.5:2.5x =0 0.5000 1.0000 1.5000 2.0000 2.5000 >> y=[2.0 1.1 0.9 0.6 0.4 0.3]y =2.0000 1.1000 0.9000 0.6000 0.4000 0.3000 >> plot(x,y)图5.1 离散点分布示意图从图5.1观察数据点分布，用二次曲线拟合。
matlab基础教程课后答案
matlab基础教程课后答案【篇一：matlab基础练习题及答案】xt>1.4 matlab操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？答：在matlab操作桌面上有五个窗口。
在每个窗口的右上角有两个小按钮，一个是关闭窗口的close按钮，一个是可以使窗口成为独立窗口的undock按钮，点击undock按钮就可以使该窗口脱离桌面成为独立窗口。
在独立窗口的desktop菜单中选择dock...项就可以将独立的窗口重新放置到桌面上。
1.5 如何启动m文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，m文件编辑/调试器将被启动。
在命令窗口中键入edit命令时也可以启动m文件编辑/调试器。
1.6 存储在工作空间中的数组能编辑吗？如何操作？1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到m文件中。
1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的set path菜单项来完成。
在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被matlab运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。
1.9 在matlab中有几种获得帮助的途径？答：（1）帮助浏览器：选择view菜单中的help菜单项或选择help菜单中的matlab help菜单项可以打开帮助浏览器。
（2）help命令：在命令窗口键入“help” 命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息。
matlab程序设计与应用第二版习题答案
matlab程序设计与应用第二版习题答案matlab程序设计与应用第二版习题答案【篇一：matlab程序设计与应用(第二版)实验答案】%实验一 matlab运算基础%第1题%（1）z1=2*sin(85*pi/180)/(1+exp(2))%（2）x=[2,1+2i;-0.45,5];z2=0.5*log(x+sqrt(1+x.^2))%(3)a=-3.0:0.1:3.0;z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)%(4)t=0:0.5:2.5;z4=t.^2.*(t=0t1)+(t.^2-1).*(t=1t2)+(t.^2-2*t+1).*(t=2t3)%第2题a=[12 34 -4;34 7 87;3 65 7];b=[1 3 -1;2 0 3;3 -2 7];a+6*ba-b+eye(size(a))a*ba.*ba^3a.^3a/bb\a[a,b][a([1,3],:);b^2]%第3题a=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25] b=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]c=a*bf=size(c)d=c(f(1)-2:f(1),f(2)-1:f(2))whos%第4题%(1):a=100:999;b=rem(a,21);c=length(find(b==0))%(2):a=lsdhksdlkklsdkl;k=find(a=aa=z);a(k)=[]%实验二 matlab矩阵分析与处理 %第1题e=eye(3);r=rand(3,2);o=zeros(2,3);s=diag([2,3]);a=[e,r;o,s];a^2b=[e,(r+r*s);o,s^2]%第2题h=hilb(5)p=pascal(5)hh=det(h)hp=det(p)th=cond(h)tp=cond(p)%第3题a=fix(10*rand(5))h=det(a)trace=trace(a)rank=rank(a)norm=norm(a)%第4题a=[-29,6,18;20,5,12;-8,8,5][v,d]=eig(a)%数学意义略%第5题方法一%(1):a=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; b=[0.95,0.67,0.52]; x=inv(a)*b%(2):b=[0.95,0.67,0.53];x=inv(a)*b%(3):cond(a)%第5题方法二a=hilb(4)a(:,1)=[]a(4,:)=[]b=[0.95,0.67,0.52];x=inv(a)*bb1=[0.95,0.67,0.53];x1=inv(a)*b1n=cond(b)n1=cond(b1)na=cond(a) %矩阵a为病态矩阵%第6题a=[1,4,9;16,25,36;49,64,81]b=sqrtm(a)c=sqrt(a) %sqrtm函数是以矩阵为单位进行计算，sqrt函数是以矩阵中的元素进行计算%实验三选择程序结构设计%第1题程序一x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];y=[]; %建立存放所有y值的矩阵for x0=xif x00x0~=-3y=[y,x0*x0+x0-6];elseif x0=0x05x0~=2x0~=3y=[y,x0*x0-5*x0+6];elsey=[y,x0*x0-x0-1];endendx%输出所有xy%输出所有y%第1题程序二x=[-5,-3,1,2,2.5,3,5];y=[];for a=1:7if x(a)0x(a)~=-3y=[y,(x(a))^2+x(a)-6];elseif x(a)=0x(a)5x(a)~=2x(a)~=3y=[y,(x(a))^2-5*x(a)+6];elsey=[y,x(a)*x(a)-x(a)-1];endend%第2题程序一x=input(请输入一个百分制成绩：);if x100|x0disp(您输入的成绩不是百分制成绩，请重新输入。
MATLAB基础及其应用教程-周开利-邓春晖课后答案,第三章
第三章习题答案1．代码：a=[1 -1 -1]; roots(a)结果：ans =-0.61801.6180２．代码：x=0:10;y=sin(x);xi=0:0.15:10; %选取了67个插值点，要增加n，只需减小步长即可y0=sin(xi); %算精确值y1=interp1(x,y,xi); %分段线性插值y2=interp1(x,y,xi,'spline'); %三次样条插值plot(xi,y0,'o',xi,y1,xi,y2,'-.')legend('精确值','分段线性插值','三次样条插值')结果：３．理论公式为：p=1.0332*exp(-(x+500)/7756)，所以拟合模型可写为：p=a*exp(-k*x+b) 式中，a, k, b为常数，两边同时取自然对数，得：log(p)=-k*x+b+log(a)问题转化为线性模型。
注意：自然对数是log(x), 以10为底的对数是log10(x)代码：clear;x=[0 300 600 1000 1500 2000];p=[0.9689 0.9322 0.8969 0.8519 0.7989 0.7491];lnp=log(p); %转化为 p 的自然对数值,模型转化为线性模型pk=polyfit(x,lnp,1); % 线性拟合,得到模型的斜率pk(1)和常数pk(2) 模型为: p=exp(pk(1)*x)*exp(pk(2))xi=0:50:2000;p0=1.0332*exp(-(xi+500)/7756); %理论值p1=exp(pk(1)*xi+pk(2)); %拟合模型值p2=interp1(x,p,xi,'spline'); %三次样条插值plot(x,p,'p',xi,p0,xi,p1,'--',xi,p2,'-.');legend('测量值','理论值','拟合值','三次样条值');format long % 数据显示格式为15位有效数字x2=0:200:2000 % 取10个点，比较差异pp1=1.0332*exp(-(x2+500)/7756) %理论值pp2=exp(pk(1)*x2+pk(2)) % 拟合值pp3=interp1(x,p,x2,'spline') % 样条插值err1=sum(abs(pp2-pp1).^2) % 拟合值的误差绝对值总和err2=sum(abs(pp3-pp1).^2) % 样条值的误差绝对值总和结果：0200400600800100012001400160018002000从图像上，都符合得很好，但很难看出差异。
matlab程序设计与应用第二版习题答案
matlab程序设计与应用第二版习题答案Matlab程序设计与应用第二版习题答案Matlab是一种强大的数学软件，广泛应用于科学计算、数据分析和工程设计等领域。
《Matlab程序设计与应用》是一本经典的教材，对于学习和掌握Matlab编程语言具有重要的意义。
本文将为大家提供《Matlab程序设计与应用第二版》中部分习题的答案，帮助读者更好地理解和应用Matlab。
第一章：Matlab基础1.1 基本操作1. a = 3; b = 4; c = sqrt(a^2 + b^2); disp(c);2. x = linspace(-pi, pi, 100); y = sin(x); plot(x, y);3. A = [1 2 3; 4 5 6; 7 8 9]; B = [9 8 7; 6 5 4; 3 2 1]; C = A + B; disp(C);1.2 控制结构1. for i = 1:10disp(i);end2. n = 0; sum = 0; while sum < 100n = n + 1;sum = sum + n;enddisp(n);3. x = 5; if x > 0disp('x is positive');elseif x < 0disp('x is negative');elsedisp('x is zero');end第二章：向量和矩阵运算2.1 向量运算1. A = [1 2 3]; B = [4 5 6]; C = A .* B; disp(C);2. A = [1 2 3]; B = [4 5 6]; C = A ./ B; disp(C);3. A = [1 2 3]; B = [4 5 6]; C = dot(A, B); disp(C);2.2 矩阵运算1. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = A * B; disp(C);2. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = B * A; disp(C);3. A = [1 2 3; 4 5 6]; B = [7 8; 9 10; 11 12]; C = A .* B; disp(C); 第三章：函数和脚本文件3.1 函数1. function y = myfunc(x)y = x^2 + 3*x + 2;end2. function [y1, y2] = myfunc(x1, x2)y1 = x1^2 + 3*x1 + 2;y2 = x2^2 + 3*x2 + 2;end3. function [y1, y2] = myfunc(x)y1 = x^2 + 3*x + 2;y2 = sin(x);end3.2 脚本文件1. x = linspace(0, 2*pi, 100); y = sin(x); plot(x, y);2. x = linspace(-10, 10, 100); y = x.^2 + 3*x + 2; plot(x, y);3. x = linspace(0, 2*pi, 100); y1 = sin(x); y2 = cos(x); plot(x, y1, x, y2);通过以上习题的答案，读者可以对Matlab程序设计的基本语法和常用函数有一个初步的了解。
matlab基础与应用教程课后答案
x=3*exp(-t); % 求各采样点样本值x
dt=t(2)-t(1); f=1/dt;
% 采样周期 % 采样频率(Hz)
X=fft(x);
% 计算x的快速傅立叶变换X
F=X(1:N/2+1); f=f*(0:N/2)/N;
% F(k)=X(k)(k=1:N/2+1) % 使频率轴f从零开始
plot(f,abs(F),'-*')
（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]
8. x=linspace(-3,3,100);y=linspace(-3,3,100); [x y]=meshgrid(x,y); fxy=-5./(1+x.^2+y.^2); i=find(abs(x)<=0.8 & abs(y)<=0.5);
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.^2; x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.^2; surf(x,y,z);shading interp; light('position',[1,0,1]); 10.
MATLAB基础及其应用教程周开利邓春晖课后答案
第三章习题及参考答案解答：>> p=[1 -1 -1];>> roots(p)ans =-0.61801.6180解答：取n=5，m=61>> x=linspace(0,2*pi,5); y=sin(x);>> xi=linspace(0,2*pi,61);>> y0=sin(xi);>> y1=interp1(x,y,xi);>> y2=interp1(x,y,xi,'spline');>> plot(xi,y0,'o',xi,y1,xi,y2,'-.');>> subplot(2,1,1); plot(xi,y1-y0);grid on >> subplot(2,1,2); plot(xi,y2-y0);grid on分段线性和三次样条插值方法与精确值之差取n=11，m=61>> x=linspace(0,2*pi,11); y=sin(x);>> xi=linspace(0,2*pi,61);>> y0=sin(xi);>> y1=interp1(x,y,xi);>> y2=interp1(x,y,xi,'spline');>> plot(xi,y0,'o',xi,y1,xi,y2,'-.');>> subplot(2,1,1); plot(xi,y1-y0);grid on>> subplot(2,1,2); plot(xi,y2-y0);grid on分段线性和三次样条插值方法与精确值之差解答：>> x=[0,300,600,1000,1500,2000];>> y=[0.9689,0.9322,0.8969,0.8519,0.7989,0.7491]; >> xi=0:100:2000;>> y0=1.0332*exp(-(xi+500)/7756);>> y1=interp1(x,y,xi,'spline');>> p3=polyfit(x,y,3);>> y3=polyval(p3,xi);>> subplot(2,1,1);plot(xi,y0,'o',xi,y1,xi,y3,'-.'); >> subplot(2,1,2);plot(xi,y1-y0,xi,y3-y0);grid on插值和拟合方法相比较，都合理，误差也相近。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
402
591
12489
A.*B
ans =
3 -108 134
68 -360 518
783 270 0
（2） A^3
ans =
-28917
240246
-4368
137883 -259101
27669
171333
252504
38673
A.^3
ans =
27
157464
8
39304
-91125
343
658503
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])
8. x=linspace(-3,3,100);y=linspace(-3,3,100); [x y]=meshgrid(x,y); fxy=-5./(1+x.^2+y.^2); i=find(abs(x)<=0.8 & abs(y)<=0.5);
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.^2; x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.^2; surf(x,y,z);shading interp; light('position',[1,0,1]); 10.
stem(t,y) subplot(2,2,4); semilogy(t,y); 5. 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);
6.
y=fun(40)/(fun(30)+fun(20)) （1）函数文件fun.m function f=fun(n) f=n+log(n^2+5); （2）函数文件fun.m function f=fun(n) a=1:n; f=sum(a.*(a+1));
第 3章 一、思考题
4. t=0:0.01:1.5; y=sqrt(3)/2*exp(-4*t).*sin(4*sqrt(3)*t+pi/3); 5. x=-10:0.01:10; y=linspace(-6,6,size(x,2)) z=x.^3+3*x.*y.^2; 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:0.01:10;
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;
exp(-0.3*a).*sin(a+0.3)
3.x=[2,4;-0.45,5];
log(x+sqrt(1+x.*x))/2
4. A=[3,54,2;34,-45,7;87,90,15];B=[1,-2,67;2,8,74;9,3,0];
（1）A*B
ans =
129
432
4197
7
-407-1052来自plot(x,x-x.^3/3/2) （2）plot(x,exp(-x.*x/2)/2/pi) （3）x=-8:0.01: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*x1-0.5; 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(0.4,1,'y=2x-0.5'); text(-0.4,0.5,'x=sin(3t)cos(t)'); text(-0.4,0.3,'y=sin(3t)sin(t)'); （2）
subplot(1,2,1); scatter(x1,y1,10)
title('y=2x-0.5'); subplot(1,2,2); scatter(x,y,10) 3.
subplot(1,21); 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);
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(0.2)*randn(1,500); （4）E=ones(size(A)); （5）A=A+30*eye(size(A)); （6）B=diag(diag(A)) 二、实验题
3 54 2 1 -2 67 34 -45 7 2 8 74 87 90 15 9 3 0 [A([1,3],:);B^2] ans =
3 54 2 87 90 15 600 183 -81 684 282 726 15 6 825 5.a=1+2i;b=3+4i; c=exp(pi*i/6) c= 0.8660 + 0.5000i c+a*b/(a+b) ans = 1.6353 + 1.8462i 第 2章 一、思考题 3.s=0; for n=0:63 s=s+2^n; end disp(s)
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,4.5e-4;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);
[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);
729000
3375
（3）A/B
ans =
-5.8845 5.3549 -0.2028
6.3554 -5.6596 4.3293
-6.6325 6.2078 9.0241
B/A
ans =
30.2855 19.2643 -8.5614
33.4394 21.1547 -9.3974 -0.7443 -0.3938 0.2830 （4）[A,B] ans =
第 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]
x=3*exp(-t); % 求各采样点样本值x
dt=t(2)-t(1); f=1/dt;
% 采样周期 % 采样频率(Hz)
X=fft(x);
% 计算x的快速傅立叶变换X
F=X(1:N/2+1); f=f*(0:N/2)/N;
% F(k)=X(k)(k=1:N/2+1) % 使频率轴f从零开始
plot(f,abs(F),'-*')
% 绘制振幅-频率图
xlabel('Frequency');
n=0:63; s=sum(2.^n) 二、实验题 1. x=input('输入一个 四位整 数：'); y=[fix(x/1000),mod(fix(x/100),10),mod(fix(x/10),10),mod(x,10)] 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>120