MATLAB基础及其应用教程-周开利-邓春晖课后答案 第六章

合集下载

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课后习题集附标准答案

MATLAB课后习题集附标准答案

第2章MATLAB概论1、与其他计算机语言相比较,MA TLAB语言突出的特点是什么?答:起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富,可扩展性强.2、MA TLAB系统由那些部分组成?答:开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB时,在选择组件窗口中哪些部分必须勾选,没有勾选的部分以后如何补安装?答:在安装MATLAB时,安装内容由选择组件窗口中各复选框是否被勾选来决定,可以根据自己的需要选择安装内容,但基本平台(即MATLAB选项)必须安装.第一次安装没有选择的内容在补安装时只需按照安装的过程进行,只是在选择组件时只勾选要补装的组件或工具箱即可.4、MATLAB操作桌面有几个窗口?如何使某个窗口脱离桌面成为独立窗口?又如何将脱离出去的窗口重新放置到桌面上?答:在MATLAB操作桌面上有五个窗口,在每个窗口的右下角有两个小按钮,一个是关闭窗口的Close 按钮,一个是可以使窗口称为独立的Undock按钮,点击Undock按钮就可以使该窗口脱离桌面称为独立窗口,在独立窗口的view菜单中选择Dock,菜单项就可以将独立的窗口重新防止的桌面上.5、如何启动M文件编辑/调试器?答:在操作桌面上选择“建立新文件”或“打开文件”操作时,M文件编辑/调试器将被启动.在命令窗口中键入edit命令时也可以启动M文件编辑/调试器.6、存储在工作空间中的数组能编辑吗?如何操作?答:存储在工作空间的数组可以通过数组编辑器进行编辑:在工作空间浏览器中双击要编辑的数组名打开数组编辑器,再选中要修改的数据单元,输入修改内容即可.7、命令历史窗口除了可以观察前面键入的命令外,还有什么用途?答:命令历史窗口除了用于查询以前键入的命令外,还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中.8、如何设置当前目录和搜索路径,在当前目录上的文件和在搜索路径上的文件有什么区别?答:当前目录可以在当前目录浏览器窗口左上方的输入栏中设置,搜索路径可以通过选择操作桌面的file菜单中的Set Path菜单项来完成.在没有特别说明的情况下,只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用,如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件,如果没有特别说明,数据文件将存储在当前目录上.9、在MA TLAB中有几种获得帮助的途径?答:(1)帮助浏览器:选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器.(2)help命令:在命令窗口键入“help”命令可以列出帮助主题,键入“help 函数名”可以得到指定函数的在线帮助信息.(3)lookfor命令:在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数.(4)模糊查询:输入命令的前几个字母,然后按Tab键,就可以列出所有以这几个字母开始的命令和函数.注意:lookfor和模糊查询查到的不是详细信息,通常还需要在确定了具体函数名称后用help命令显示详细信息.第3章 MATLAB 数值运算3.1在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375,并将其赋予变量a ? 答:在Command Window 窗口输入操作:>> a=[5 7 3;4 9 1]3.2有几种建立矩阵的方法?各有什么优点? 答:(1)直接输入法,如a=[1 2 3 4],优点是输入方法方便简捷;(2)通过M 文件建立矩阵,该方法适用于建立尺寸较大的矩阵,并且易于修改; (3)由函数建立,如y=sin(x),可以由MATLAB 的内部函数建立一些特殊矩阵; (4)通过数据文件建立,该方法可以调用由其他软件产生数据.3.3在进行算术运算时,数组运算和矩阵运算各有什么要求?答:进行数组运算的两个数组必须有相同的尺寸.进行矩阵运算的两个矩阵必须满足矩阵运算规则,如矩阵a 与b 相乘(a*b )时必须满足a 的列数等于b 的行数.渗釤呛俨匀谔鱉调硯錦。

MATLAB基础及其应用教程-周开利-邓春晖课后答案,第三章

MATLAB基础及其应用教程-周开利-邓春晖课后答案,第三章

第三章习题答案1.代码:a=[1 -1 -1]; roots(a)结果:ans =-0.61801.61802.代码: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('精确值','分段线性插值','三次样条插值')结果:3.理论公式为: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编程与应用习题和一些参考答案
5.利用 rand 产生 10 个随机数,利用 for 循环对其进行排序(从大到小)。 答:function pailie
n=rand(1,10) for j=1:10
for i=j+1:10 if n(j)<n(i) a=n(j); n(j)=n(i); n(i)=a; end
end end disp('重新排列后:') n
9.建立一个字符串向量,然后对该向量做如下处理: (1) 取第 1~5 个字符组成的子字符串。 (2) 将字符串倒过来重新排列。 (3) 统计字符串中小写字母的个数。
>> a='abcdef123456'; >> subch=a(1:5) subch = abcde >> revch=a(end:-1:1) revch = 654321fedcba>> k=find(a>='a'&a<='z'); >> a(k)=a(k)-('a'-'A'); >> char(a)
>> z=0:0.01:10;
>> x=z.*sin(3*z);
>> y=z.*cos(3*z);
>> plot3(x,y,z)
12.设 z x e2 (x2 y2 ) ,画出定义域 x=[-2,2],y=[-2,2]内的曲面图。
>> [x,y]=meshgrid([-2:0.1:2]); >> z=x.^2.*exp(-x.^2-y.^2); >> mesh(x,y,z) 13.设 z=0.05x-0.05y+0.1; 画出 z 的曲面(平面)图。 >> [X,Y]=meshgrid(1:10:200,1:10:200); >> Z=0.05*X-0.05*Y+0.1; >> surf(X,Y,Z)

MATLAB基础及其应用教程周开利邓春晖课后答案

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插值和拟合方法相比较,都合理,误差也相近。

matlab课后习题答案

matlab课后习题答案

第2章 MATLAB 矩阵运算基础2.1 在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375,并将其赋予变量a ? >> a=[5 7 3;4 9 1]2.5 计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。

>> a=[5 3 5;3 7 4;7 9 8]; >> b=[2 4 2;6 7 9;8 3 6];>> a+bans =7 7 7 9 14 13 15 12 142.6 求⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。

>> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i]; >> x’ans =4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 -5.0000i 7.0000 +6.0000i 2.0000 +7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i2.7 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。

>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8]; >> a.*b ans =12 36 3 8 42 402.9 对于B AX =,如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ,求解X 。

>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’;-0.5118 4.0427 1.33182.10 已知:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ,分别计算a 的数组平方和矩阵平方,并观察其结果。

matlab基础与应用教程课后答案

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之间实施加减乘除及幂运算。

MATLAB基础及其应用教程-周开利-邓春晖课后答案 第六章

MATLAB基础及其应用教程-周开利-邓春晖课后答案 第六章

第6章习题解答1. 在命令窗口中输入,一次只能执行一行命令,脚本文件是命令的集中,只需输入脚本文件名,就可以按顺序一次执行所有命令。

2.脚本文件是命令的集中,只需输入脚本文件名,就可以一次执行所有命令,没有输入参数,没有输出参数。

文件名可任意取,只需符合变量的定义要求即可。

脚本文件中生成的变量是全局变量。

函数文件可以接受参数,也可以返回参数,一般需要由其他语句调用。

首行语句必须以function开始,有函数名,文件名必须与函数名一致。

函数空间中生成的变量是局部变量。

3.function result=project1(year)if mod(year,400)==0 || (mod(year,4)==0 && mod(year,100)~=0)result='是闰年';elseresult='不是闰年';end4.function min=comp(x1,x2)min=x1;if x1>x2min=x2;end5.(1)在命令窗口输入:k=0;for i=-1000:1000;k=k+1;end;k 即可得出:2001步。

(2)方法同上,10步6.(1)循环进行的条件是var不是10的倍数,循环结束时,var必是10的倍数,即10(2)每次循环的结果是:4,16,256,此时,不满足循环条件而退出循环,即var为256(3)循环条件是var>100,而var的值为3,不满足条件,循环不执行,所以var值不变。

7.方法1:clear;A=[1 2 3;4 5 6];fid=fopen('c:\a1.txt','wt');fprintf(fid,'%d ',A(1,:));fprintf(fid,'\n');fprintf(fid,'%d ',A(2,:));fclose(fid);方法2:clear;a=[1 2];A=[1 2 3 4 5 6]B=reshape(A,2,2) dlmwrite('c:\a2.txt',B)。

Matlab基础与应用教程答案

Matlab基础与应用教程答案

第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基础及其应用教程》笔记

《matlab基础及其应用教程》笔记

《matlab基础及其应用教程》笔记前提:已经学过matlab很久了,但是将近半年没有再碰,打算接下去继续学习matlab,并在帮忙增加国内关于matlab资料较少的simmechanics物理建模仿真部分,因此继续开始学matlab,打算重新看基本基础的书回顾一下,并做简单笔记。

plus:括号内是自己加的,自己的经验以及想法等第一章:基础1、b='Hi,Miss Black';可以直接用于表示字符串2、cd可以设置路径;和DOS系统下一样,但是不够直观,推荐直接使用GUI的路径设置3、clf:清除图像内容;cla:清除坐标内容(GUI会用到);clc,clear4、lookfor:寻找函数名称;doc:直接打开某个函数的帮助文档(doc quad:打开积分的帮助文档)第二章:matlab语言基础1、matlab主要数据类型.pdf(matlab的数据类型,十分重要)2、&与&&区别:&左右两边都比较;&&只要左边不对,就不进行比较(表示右边有错误也不会管)3、各种函数比较的最好加上(),以防止因为优先级的问题而被吞4、线性等分:linspace;对数等分:logspace5、矩阵标志:①,单下标;②,多下标;③,:表示整行或者整列;④,删除直接使用[];6、cat:行/列拼接;reshape:矩阵变换;repmat:行列拼接;7、一般的函数:matlab函数(自己总结).xls里面有8、多维数组(三维,图片数组)第三章:matlab运算1、polyadd:相加;conv:相乘/卷积;deconv:相除2、polyval:多项式带入求解;polyvalm:矩阵带入求解3、roots:多项式求根;4、poly2sym:多项式转为字符形式;poly:构造多项式5、interpX:插值与拟合(和《数值计算方法》对应,X不同而采取不同的差值拟合方法)6、polyfit:生成多项式拟合(偶尔两端会出现Roung现象,使整个设计出现偏差)7、diff:微分;quad:积分8、cumsum:积分矩形计算法;trapz:积分梯形计算法;quad:辛普森积分法;quad1:科茨积分法;9、rref:解矩阵方程10、jacobi:雅克比迭代法;gseidel:Gauss-Seidel 迭代法;11、sparse:稀疏矩阵;full:全矩阵;spdiags:对角稀疏矩阵12、spalloc:分配存储空间(对稀疏矩阵而言)13、Euler:欧拉法解微分(1:前插;2:后插)14、ode23:二三阶龙科库塔公式解偏微分;ode45:三四阶龙科库塔公式第四章:结构数组与细胞数组(感觉暂时用到比较少,所以只是简单浏览了一下)1、结构数组和C语言的一样,后面加.2、或者使用struct构建结构数组3、结构数组相关函数.pdf4、rmfield:删除域;gerfield:获得结构数组中的值(GUI中十分有用)5、setfield:重新设定值6、cell:构建细胞数组7、cell2struc/struc2cell:结构数组与细胞数组的互相转换8、mat2cell/cell2mat:细胞数组与普通数组的转换第五章:符号运算1、sym():定义单个符号常量;syms:定义多个符号常量2、sym:也可以定义矩阵符号3、合并同类项(collect)、多项式展开(expand)、因式分解(factor)、一般化简(simplify)、不定化简(simple)、通分(numden)和书写格式美化(pretty)。

MATLAB基础及其应用教程-周开利-邓春晖课后标准答案

MATLAB基础及其应用教程-周开利-邓春晖课后标准答案

MATLAB基础及其应⽤教程-周开利-邓春晖课后标准答案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分段线性和三次样条插值⽅法与精确值之差解答:>> 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插值和拟合⽅法相⽐较,都合理,误差也相近。

Matlab 第6章习题A习题具有题解.docx

Matlab 第6章习题A习题具有题解.docx

A 习题(具有题解)A 6-1试判断下述系统的可控性及可观性。

y(k) = [2 —4]x 伙)能否找到一组控制序列,使系统从原点到达[1 1 If,解释为什么。

x(k +1)=0.5 0-0.5 0.25x 伙)+u(k)4C CFT0.5-0.5"[00.25_ 42 -4 1 -2rank= 1,系统不可观。

A 6-2下述连续系统被采样,求离散传递函数,并确定T 为何值时系统不可控,试说明之。

宀、 2(S + 5)G (s ) = ---------------?+105 + 29解:A (5)= 52 4-105 + 25 + 4 = 0 ; 512 =-5± j2;所以込一$2 = /4。

依要求可知,若S\-S2=jk 〒,采样系统不可控。

故有T = k7i/2时系统不可控。

Gd) =2[z~ ~ze 5T cos 2T] z 2-2ze-5T cos2T + e-10r,如当T =冗丨2时心(f 二吕'发生零极对消。

A 6-3给定下述系统「0 1 2「兀1(切'o'x(k + l) = 0 0 3x 2(k) + 10 0 0 x 2(k)(1) (2) 该控制序列最少步数是多少。

解:= [FG G] = <rank=2,系统可控试确定一组控制序列,使系统从%(o )= [i 1 if 达到原点。

_0 1 2_■f'0_'3_'0_解:1) x(l)=0 0 3 1 + 1 w(0) = 3 + 1 M (0)0 0 1如取u(0)=-3,则如取u(l)=-O,则x(2) = [0] o 表明u(0)=3、u ⑴=0,可使系统从x(0) = [l 1『达到原点。

2) 显然最少步数N=2。

_0 1 03) 因为W R =[F-G FG G]= 0 0 10 0 0移矩阵F 可见,X3(k)不受u(k)影响,且与其他状态无关,所以不能通过u(k)改变其状态。

中南大学matlab课后答案-第六章

中南大学matlab课后答案-第六章

实验指导二1,(1)>> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];B=[,,]';x=A\B (2)>> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];B=[,,]';x=A\B (3)>> cond(A)2,(1)建立函数文件,命令如下;Function fx=funx(x)fx=x^41+x^3+1;$调用fzero函数求根,命令如下;>> z=fzero(@funx,-1)(2) 建立函数文件,命令如下;Function fx=sin(x)fx=x-(sin(x))/x;调用fzero函数求根,命令如下;>> z=fzero(@sin,(3) 建立函数文件,命令如下;function q=myfun(p)x=p(1);、y=p(2);z=p(3);q(1)=sin(x)+y^2+log(z)-7;q(2)=3*x+2^y-z^3+1;q(3)=x+y+z-5;调用fsolve函数求根,命令如下;>> options=optimset('Display','off');x=fsolve(@myfun,[1,1,1]',options) 3,(1) 建立函数文件,命令如下;function yp=funt(t,y)>yp=-y*+sin(10*t));求微分方程,程序如下:>> t0=0;tf=5;y0=1;[t,y]=ode23(@funt,[t0,tf],y0)(2) 建立函数文件,命令如下;function yp=funr(t,y)yp=cos(t)-(y)/(1+t^2);求微分方程,程序如下:>> t0=0;tf=5;y0=1;[t,y]=ode23(@funr,[t0,tf],y0)4,建立函数文件命令如下:{function fx=mymax(x)fx=-1*(1+x^2)/(1+x^4);求最大值,程序如下:>> [x,y]=fminbnd(@mymax,0,2)5,编写目标函数M文件,命令如下:function f=fop(x)f=-1*(x(1)^(1/2)+x(2)^(1/2)+x(3)^(1/2)+x(4)^(1/2));设定约束条件,并调用fmincon函数求解此约束最优化问题,程序如下:>>x0=[200,200,200,200];%A=[1,0,0,0;,1,0,0;,,1,0;,,,1];b=[400,440,484,];Ib=[0,0,0,0];options=optimset('Display','off');[x,y]=fmincon(@fop,x0,A,b,[],[],Ib,[],[],options)思考练习1,(1)矩阵求逆:>> A=[2,3,5;3,7,4;1,-7,1];b=[10,3,5]';x=inv(A)*b矩阵除法:>> A=[2,3,5;3,7,4;1,-7,1];b=[10,3,5]';x=A\b!矩阵分解:>> A=[2,3,5;3,7,4;1,-7,1];b=[10,3,5]';[L,U]=lu(A);x=U\(L\b) (2)矩阵求逆:>>A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];b=[-4,13,1,11]';x=inv(A)*b矩阵除法:>>A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];b=[-4,13,1,11]';x=A\b矩阵分解:>> A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];b=[-4,13,1,11]'; [L,U]=lu(A);x=U\(L\b)2,(1)建立函数文件,命令如下:function fx=fun1(x)fx=3*x+sin(x)-exp(x);调用fzero函数求根,命令如下;>> y=fzero(@fun1,、(2) 建立函数文件,命令如下:function fx=fun2(x)fx=1-(1/x)+5;调用fzero函数求根,命令如下;>> y=fzero(@fun2,1)(3) 建立函数文件,命令如下;function q=fun3(p)x=p(1);y=p(2);q(1)=x^2+y^2-9;'q(2)=x+y-1;调用fsolve函数求根,命令如下;>> options=optimset('Display','off');x=fsolve(@fun3,[1,1]',options) 3,(1)建立函数文件,命令如下:function ydot=fun5(t,y)ydot(1)=(2-3*y(2)-2*t*y(1))/(1+t^2);ydot(2)=y(1);ydot=ydot';求解微分方程,命令如下:'>> t0=0;tf=5;x0=[0,1];[t,y]=ode45(@fun5,[t0,tf],x0);[t,y](2) (1)建立函数文件,命令如下:function ydot=fun6(t,y)ydot(1)=cos(t)+(5*y(1)*cos(2*t))/(t+1)^2-y(2)-y(3)/(3+sin(t)); ydot(2)=y(1);ydot(3)=y(2);ydot=ydot';求解微分方程,命令如下:>> t0=0;tf=5;x0=[0,1];[t,y]=ode45(@fun5,[t0,tf],x0);[t,y]4,】建立函数文件命令如下:function fx=max(x)fx=-1*(sin(x)+cos(x^2));求最大值,程序如下:>> [x,y]=fminbnd(@mymax,0,pi)5,编写目标函数M文件,命令如下:function f=topm(x)f=-1*x*(3-2*x)^2;设定约束条件,并调用fmincon函数求解此约束最优化问题,程序如下:>>x0=[0];A=[1];b=[];Ib=[0];options=optimset('Display','off');[x,y]=fminco n(@top,x0,A,b,[],[],Ib,[],[],options)。

唐昌建(MATLAB编程基础及应用)第六章

唐昌建(MATLAB编程基础及应用)第六章

唐昌建(MATLAB编程基础及应⽤)第六章第六章数值计算⽅法与数据第⼀节⾮线性数值计算⼀、⾮线性微分⽅程的求解ordinary differential equation (ODE ) 1.基本问题(1)⾮线性问题。

(2)基本命令:⾮刚性常微分⽅程(组)ode23,ode45,ode113 2.基本解法:(1)建⽴标准微分⽅程(组)(,)d yF y t dx =⽬的:将所有不同变量归为同⼀变量。

⾼阶微分⽅程归为⼀阶微分⽅程。

例如:y''+2y'-y=0 令 y 1=y;y 2=y'原式变为: y 2'=y 1-2y 2;y 1'=y 2(2)建⽴ode 相应的函数m ⽂件,格式如function ff=fun(t,y ) ff=[(,)F y t ](3)调⽤fun ⽤ode 函数求解可在命令窗进⾏也可以在m ⽂件中进⾏。

ode 的调⽤格式:[t,y ]=ode45('fun',[⾃变量的范围],[初值列阵]) (4)例⼦练习 6-1(1)2222dyy x xdx =-++ y(0)-1(fun1 exno19)(2)捕⾷者与被捕⾷者函数关系(x 是被捕⾷者)0.01dx x xy dt =+ 0.02dyx xy dt =-+ x(0)=30 y(0)=20(fun2exno20)(3)222(1)0d y dy y y dx dx µ=-+=(fun3 exno21)(0)0dydx = y(0)=1变换为:12dy y dx = 22121(1)dy y y y dx µ=-+y 1(0)=1 y 2(0)=0fun1function dy=fun1(x,y) dy=[-2*y+2*x.^2+2*x]e19[x,y]=ode23('fun1',[0,0.5],1)%⼀步⼀步的解,相当占⽤资源 plot(x,y)fun2function f=fun2(t,y)f=[y(1)-0.01*y(1)*y(2);-y(2)+0.02*y(1)*y(2)]e20ts=0:0.1:20 y0=[30,20][t,y]=ode45('fun2',ts,y0) [t,y]subplot(221)plot(t,y(:,1),'r',t,y(:,2),'b'),gtext('x'),gtext('y') subplot(222)plot(y(:,1),y(:,2)),gtext('x'),gtext('y')fun3function f=fun3(t,y)f=[y(2);2*(1-y(1)^2)*y(2)-y(1)]e21(⾃编)ts=0:0.1:20y0=[1,0][t,y]=ode113('fun3',ts,y0)[t,y]subplot(221)plot(t,y(:,2),'r',t,y(:,1),'b')⼆、⾮线性函数的最⼩值1.单变量(可代参数)a)[x,y,po1,po2]=fminbnd(fun,x1,x2,options,p1,p2,...)b)fun:可带参数的⽬标函数c)x2,x2:⾃变量的求值范围d)options:输⼊的优化参数,包括四个域display,maxfunevals,maxlter,tolx.这四个域由函数optimset()来传递。

matlab课后习题答案到6章

matlab课后习题答案到6章

(6) 将含有 12 个元素的向量 x 转换成3 4 矩阵。

11习 题二答: reshape(x,3,4);E|D=1 1 ,11 1.如何理解“矩阵是 MATLAB 最基 (7) 求一个字符串的 ASCII 码。

本的数据对象”?答: abs(‘123'); 或00~E|~D=1 0答:因为向量可以看成是仅有一行或一列 double( ‘ 123 ' );的矩阵, 单个数据 (标量) 可以看成是仅 (8) 求一个 ASCII 码所对应的字含一个元素的矩阵, 故向量和单个数据都 符。

00可以作为矩阵的特例来处理。

答: char(49);find(A>=10&A<25)=[1; 5] 。

因此,矩阵是 MATLAB 最基本、4. 下列命令执行后, L1 、L2 、 L3 、6.当 A=[34, NaN, Inf, -Inf, -pi, eps, 最重要的数据对象。

L4 的值分别是多少?0]时,分析下列函数的执行结果:2.设A 和B 是两个同维同大小的矩 A=1:9;B=10-A;...all(A) , any(A) , isnan(A) , isinf(A)阵,问:L1=A==B;isfinite(A) 。

(1) A*B 和 A.*B 的值是否相L2=A<=5;答 : all(A) 的值为 0等?L3=A>3&A<7;any(A) 的值为 1答: 不相等。

L4=find(A>3&A<7);isnan(A) 的值为[ 0, 1, 0, 0, 0, 0,(2) A./B 和 B.\A 的值是否相答: L1 的值为 [0, 0, 0, 0, 1, 0, 0, 0, 0]等?0]isinf(A) 的值为[ 0, 0, 1, 1, 0, 0,答:相等。

L2 的值为[1, 1, 1, 1, 1, 0,0](3) A/B 和 B\A 的值是否相0, 0, 0]isfinite(A) 的值为 [1, 0, 0, 0, 1,等?L3 的值为[0, 0, 0, 1, 1, 1, 1, 1]答: 不相等。

MATLAB8.5教程第6章 程序设计

MATLAB8.5教程第6章 程序设计

6.4.1 for循环
• for 语句可完成指定次重复的循环,这是广泛应用的语句。 • 1.for 语句的简单形式
• 格式:for variable=expression • statements • end • 说明:每当循环变量variable取expression中的一个值时, 就执行循环体statements一次,直到取完expression中的 值。expression常用数组形成表示。 • 2.多重嵌套的for 循环
6.1.4 运算符优先级
6.2 向量化和预分配
• 1.向量化 • 为使MATLAB高效运行,在M文件中最好把算法向量化。即将程序语 言中的for循环,可用向量或矩阵运算来代替。 • 譬如在求之和时,我们分别使用for循环语句和向量化编写程序如下: • %for循环程序 • s=0; • for k=1:50 • s=s+1/(2*k-1); • k=k+1; • end • s • %向量化程序 • >> x=1:2:99; • >> s=sum(1./x)
6.6 试探结构
try- catch语句是一种错误处理语句,提供一种捕获错误的方 法。 格式:try statements1 catch statements2 end 说明:通常情况下,只执行statements1,但当在执行 statements1 语句发生错误时,catch控制块就可捕获它,会 执行statements2,这样可以在statements2 中对错误作适当 处理。可调用lasterr函数查询出错原因。try和catch控制块中 的语句之间用逗号隔开。
6.3 顺序语句
• 顺序语句是指依次按顺序执行程序的各条语句,它不需要 任何特殊的流控制。 • 格式:expression • variable=expression %执行表达式命令,显示表达式值 %将表达式赋值给变量variable

MATLAB基础应用案例教程第六章

MATLAB基础应用案例教程第六章
• (6) BackgroundColor 和 EdgeColor 属性 设置文本对象的背景颜色和边框线的颜色,可取值 none ( 默认 值) 或颜色字母。
• (7) HorizontalAlignment 属性设置文本与指 定点的相对位置,可取值 left ( 默认值)、center 或 ri ght。
mal ( 默认值)、 bold、 light 或 demi。 • (4) FontAngle 属性设置斜体文字模式,取值可以是 no
rmal ( 默认值)、 italic 或 oblique。
上一页 下一页 返回
6.1 二维绘图功能
• (5) Rotation 属性设置字体旋转角,取值是数值量,默认 值为 0,取正值表示逆时针旋转,取负值表示顺时针旋转。
是网格点上的高度矩阵,c用于指定在不同高度下的颜色范围。 • (3) 省略 c 时,c = z,即颜色的设定正比于图形的高度。 • (4) 当x、 y 是向量时,x 的长度必须等于 z 矩阵的列,y 的长
度必须等于必须等于 z 的行,x、 y 向量元素的组合构成网格点的 坐标 (x。y),z 坐标则取自矩阵z,然后绘制三维曲线。
上一页 下一页 返回
6.1 二维绘图功能
• 2 函数图函数 • 语法格式:
• 说明: • (1) fun 必须是 m 文件的函数名或是独立变量为 x 的字符串。 • (2) fplot 函数不能画参数方程和隐函数图形,但在一个图上
可以画多个图形。
上一页 下一页 返回
6.1 二维绘图功能
• 【例 6 - 11】 建立函数文件 myfun1 m, 在 [ - 1。2] 区间上绘制 y = e2x + sin(3x2 ) 曲线。
有的允许属性。 • 2 获取图形对象属性

《MATLAB程序设计教程(第二版)》第6章__MATLAB解方程与最优化问题求解

《MATLAB程序设计教程(第二版)》第6章__MATLAB解方程与最优化问题求解

例6-7 分别用Jacobi迭代和Gauss-Serdel迭代法求解 下列线性方程组,看是否收敛。 命令如下: a=[1,2,-2;1,1,1;2,2,1]; b=[9;7;6]; [x,n]=jacobi(a,b,[0;0;0]) [x,n]=gauseidel(a,b,[0;0;0])
6.2 非线性方程数值求解 6.2.1 单变量非线性方程求解 在MATLAB中提供了一个fzero函数,可以用来求单变量非 线性方程的根。该函数的调用格式为: z=fzero('fname',x0,tol,trace) 其中fname是待求根的函数文件名,x0为搜索的起点。一个 函数可能有多个根,但fzero函数只给出离x0最近的那个根。 tol控制结果的相对精度,缺省时取tol=eps,trace指定迭代 信息是否在运算中显示,为1时显示,为0时不显示,缺省 时取trace=0。
例6-9 求下列非线性方程组在(0.5,0.5) 附近的数值解。 (1) 建立函数文件myfun.m。 function q=myfun(p) x=p(1); y=p(2); q(1)=x-0.6*sin(x)-0.3*cos(y); q(2)=y-0.6*cos(x)+0.3*sin(y); (2) 在给定的初值x0=0.5,y0=0.5下,调用fsolve函数求方程的 根。 x=fsolve('myfun',[0.5,0.5]',optimset('Display','off')) x= 0.6354 0.3734
例6-10 设有初值问题,试求其数值解,并与精确解 相比较。 (1) 建立函数文件funt.m。
function yp=funt(t,y) yp=(y^2-t-2)/4/(t+1); (2) 求解微分方程。 t0=0;tf=10; y0=2; [t,y]=ode23('funt',[t0,tf],y0); %求数值解 y1=sqrt(t+1)+1; t' y' y1' %求精确解

matlab答案第六章

matlab答案第六章

matlab答案第六章第六章1. 利用MATLAB提供的randn函数生成符合正态分布的10×5随机矩阵A,进行如下操作:(1) A各列元素的均值和标准方差。

(2) A的最大元素和最小元素。

(3) 求A每行元素的和以及全部元素之和。

(4) 分别对A的每列元素按升序、每行元素按降序排列。

答:clear all; close all; clc;A=randn(10, 5);meanA=mean(A); %(1)A各列元素的均值stdA=std(A); %(1)A各列元素的标准方差maxA=max(max(A)); %(2)A的最大元素minA=min(min(A)); %(2)A的最小元素rowsumA=sum(A, 2); %(3)A每行元素的和sumA=sum(rowsumA); %(3)A全部元素之和sort1=sort(A); %(4)A的每列元素按升序排列sort2=sort(A, 2, 'descend'); %(4)A的每行元素按降序排列2. 按要求对指定函数进行插值和拟合。

(1) 按表6.1用3次样条方法插值计算0~90D范围内整数点的正弦值和0~75D范围内整数点的正切值,然后用5次多项式拟合方法计算相同的函数值,并将两种计算结果进行比较。

表6.1 特殊角的正弦和正切值表α度0 15 30 4560 75 90sinα0 0.2588 0.50000.7071 0.8660 0.9659 1.0000 tanα0 0.2679 0.57741.0000 1.7320 3.7320(2) 按表6.2用3次多项式方法插值计算1~100之间整数的平方根。

表6.2 1~100内特殊值的平方根表N 1 4 9 16 25 36 49 64 81 100N的平方根 1 2 3 4 5 6 7 8 9 10答:(1) 程序设计:clear all; close all; clc;alpha1=0:15:90;sin_alpha1=sin(alpha1*pi/180); %精确正弦值plot(alpha1, sin_alpha1, 'k:p'); hold on; %绘精确正弦曲线alpha2=0:90;sin_Y1=interp1(alpha1, sin_alpha1, alpha2, 'spline'); %3次样条正弦插值plot(alpha2, sin_Y1, 'r-*'); hold on; %绘3次样条插值正弦曲线P1=polyfit(alpha1, sin_alpha1, 5); %5次多项式拟合sin_Y2= polyval(P1, alpha2); %5次多项式求值plot(alpha2, sin_Y2, 'b-o'); %绘5次多项式插值正弦曲线legend('精确正弦值', '3次样条正弦插值', '5次多项式正弦插值'); title('正弦值比较'); alpha3=0:15:75;tan_alpha3=tan(alpha3*pi/180); %精确正切值figure, plot(alpha3, tan_alpha3, 'k:p'); hold on; %绘精确正切曲线alpha4=0:75;tan_Y1=interp1(alpha3, tan_alpha3, alpha4, 'spline'); %3次样条正切插值plot(alpha4, tan_Y1,'r-*'); hold on; %绘3次样条正切曲线P2=polyfit(alpha3, tan_alpha3, 5); %5次多项式拟合tan_Y2= polyval(P2, alpha4); %5次多项式求值plot(alpha4, tan_Y2, 'b-o'); %绘5次多项式插值正弦曲线legend('精确正切值', '3次样条正切插值', '5次多项式正切插值'); title('正切值比较');(2)程序设计:clear all; close all; clc;X=[1, 4, 9, 16, 25, 36, 49, 64, 81, 100]; Y=1:10;X1=1:100; Y1=interp1(X, Y, X1,'cubic');plot(X, Y, 'r:o'); hold on; %绘精确平方根曲线plot(X1, Y1, 'k-x'); %绘3次多项式插值平方根曲线legend('精确平方根', '3次多项式插值');3. 已知一组实验数据如表6.3所示。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
相关文档
最新文档