MATLAB第五章习题参考答案

clear;

x=-2:0.01:2;

for ii=1:401

if(x(ii)<0)

y(ii)=x(ii)+1;

elseif (x(ii)<1)&(x(ii)>=0)

y(ii)=1;

elseif(x(ii)>=1)

y(ii)=x(ii)^3;

end

end

plot(x,y);

2:

function y=fenduan(x)

if(x<0)

y=x+1;

elseif (x<1)&(x>=0)

y=1;

elseif(x>=1)

y=x^3;

end

3:

建立文件arcsin.m如下:

function y=arcsin(x)

if abs(x)>1

disp('输入参数的绝对值必须小于1');

return;

end

n=0;

y=0;

u=x;

while u>eps

u=(jiecheng(2*n)*x^(2*n+1))/(2^(2*n)*(jiecheng(n)^2*(2*n+1)));

y=y+u;

n=n+1;

end

建立文件jiecheng.m如下:

function y=jiecheng(n)

if(n<0)

disp('输入参数必须为大于或等于0的整数');

return;

if n==0 | n==1

y=1;

else

y=n*jiecheng(n-1);

end

4:

建立shuliang.m文件如下:

function [tu,ji]=shuliang(tou,jiao)

tu=jiao/2-tou;

ji=2*tou-jiao/2;

在命令窗口中执行:>>[a,b]=shuliang(36,100)

5:

for i=2:999

ge=rem(i,10);

shi=rem(fix(i/10),10);

bai=fix(i/100);

if rem((ge+shi+bai),2) & isprime(i)

disp(i);

end

end

6:

function y=issushu(n)

y=1;

for ii=2:fix(n/2)

if(rem(n,ii)==0)

y=0;

end

end

MATLAB第五章习题参考答案

clear; x=-2:0.01:2; for ii=1:401 if(x(ii)<0) y(ii)=x(ii)+1; elseif (x(ii)<1)&(x(ii)>=0) y(ii)=1; elseif(x(ii)>=1) y(ii)=x(ii)^3; end end plot(x,y); 2: function y=fenduan(x) if(x<0) y=x+1; elseif (x<1)&(x>=0) y=1; elseif(x>=1) y=x^3; end 3: 建立文件arcsin.m如下: function y=arcsin(x) if abs(x)>1 disp('输入参数的绝对值必须小于1'); return; end n=0; y=0; u=x; while u>eps u=(jiecheng(2*n)*x^(2*n+1))/(2^(2*n)*(jiecheng(n)^2*(2*n+1))); y=y+u; n=n+1; end 建立文件jiecheng.m如下: function y=jiecheng(n) if(n<0) disp('输入参数必须为大于或等于0的整数'); return;

if n==0 | n==1 y=1; else y=n*jiecheng(n-1); end 4: 建立shuliang.m文件如下: function [tu,ji]=shuliang(tou,jiao) tu=jiao/2-tou; ji=2*tou-jiao/2; 在命令窗口中执行:>>[a,b]=shuliang(36,100) 5: for i=2:999 ge=rem(i,10); shi=rem(fix(i/10),10); bai=fix(i/100); if rem((ge+shi+bai),2) & isprime(i) disp(i); end end 6: function y=issushu(n) y=1; for ii=2:fix(n/2) if(rem(n,ii)==0) y=0; end end

matlab课后习题及答案

第一章 5题 已知a=4.96,b=8.11,计算 ) ln(b a e b a +-的值。 解:clear clc a=4.96; b=8.11; exp(a-b)/log(a+b) ans = 0.0167 6题 已知三角形的三边a=9.6,b=13.7, c=19.4,求三角形的面积。提示:利用海伦公式area =))()((c s b s a s s ---计算,其中S=(A+B+C)/2. 解:clear clc a=9.6; b=13.7; c=19.4; s=(a+b+c)/2 area=sqrt(s*(s-a)*(s-b)*(s-c)) s = 21.3500 第二章 8题 已知S=1+2+2^2+2^3+……+2^63,求S 的值 解:clear clc S=0; for i=0:1:63 S=S+2^i; end S S = 1.8447e+019 9题 分别用for 和while 循环结构编写程序,计算∑ =-100 1 n 1n 2) (的值。 解:clear clc s=0;

for n=1:100 s=s+(2*n-1); end s s = 10000 clear clc n=1; s=0; while n<=100 s=s+(2*n-1); n=n+1; end s s = 10000 第三章 2题 在同一坐标下绘制函数 x , ,2 x -,2 x xsin(x)在()∏∈,0x 的曲线。 解:clear clc x=0:0.2:pi; y1=x; y2=x.^2; y3=-(x.^2); y4=x.*sin(x); plot(x,y1,'-' ,x,y2,'-' ,x,y3,'-' ,x,y4,'-')

Get清风MATLAB教程a第5章习题解答张志涌

MATLAB教程2012a第5章习题解答-张志涌

第5章 数据和函数的可视化 习题5及解答 1 椭圆的长、短轴2,4==b a ,用“小红点线〞画椭圆⎩⎨⎧==t b y t a x sin cos 。〔参 见图p5-1〕 〖解答〗 clf a=4;b=2; t=0:pi/80:2*pi; x=a*cos(t); y=b*sin(t); plot(x,y,'r.','MarkerSize',15) axis equal xlabel('x') ylabel('y') shg -4 -3 -2 -1 01 2 3 4 -3-2 -1 1 2 3 x y 2 根据表达式θρcos 1-=绘制如图p5-2的心脏线。〔提示:采用极 坐标绘线指令polar 〕 〖解答〗 clf theta=0:pi/50:2*pi; rho=1-cos(theta); h=polar(theta,rho,'-r'); %极坐标绘线指令。h 是所画线的图柄。

set(h,'LineWidth',4) %利用set 设置h 图形对象的“线宽〞 axis square %保证坐标的圆整性 0.5 1 1.5 2 30 21060 240 90270 120 300 150 330 180 0ρ=1-cos θ 3 A,B,C 三个城市上半年每个月的国民生产总值如见表p5.1。试画 出如图p5-3所示的三城市上半年每月生产总值的累计直方图。 表p5.1 各城市生产总值数据〔单位:亿元〕 城市 1月 2月 3月 4月 5月 6月 A 170 120 180 200 190 220 B 120 100 110 180 170 180 C 70 50 80 100 95 120 〖目的〗 ● 借助MATLAB 的帮助系统,学习直方图指令polar 的使用。 ● bar 指令常用格式之一:bar(x,Y,'style') 。x 是自变量列向量;Y 是与x 行数相同的矩阵,Y 的每一行被作为“一组〞数据;style 取stacked 时,同一组数据中每个元素对应的直方条被相互层叠。 ● 在本例中,Y 中的一列代表一个城市。绘图时,各列的颜色,取自“色图〞;图形窗的默认色图是jet 。

第5章MATLAB绘图_习题答案

第5章MATLAB绘图 习题5 一、选择题 1.如果x、y均为4×3矩阵,则执行plot(x,y)命令后在图形窗口中绘制()条曲线。D A.12B.7C.4D.3 2.下列程序的运行结果是()。A x=0:pi/100:2*pi; forn=1:2:10 plot(n*sin(x),n*cos(x)) holdon end axissquare A.5个同心圆B.5根平行线 C.一根正弦曲线和一根余弦曲线D.5根正弦曲线和5根余弦曲线3.命令text(1,1,'{\alpha}+{\beta}')执行后,得到的标注效果是()。C A.{\alpha}+{\beta}B.αβ}C.α+βD.αβ 4.subplot(2,2,3)是指()的子图。A A.两行两列的左下图B.两行两列的右下图 C.两行两列的左上图D.两行两列的右上图 x的曲线绘制成直线,应采用的绘图函数是()。C 5.要使函数y=2e A.polarB.semilogxC.semilogyD.loglog 6.下列程序的运行结果是()。B [x,y]=meshgrid(1:5); surf(x,y,5*ones(size(x))); A.z=x+y平面B.与xy平面平行的平面 C.与xy平面垂直的平面D.z=5x平面 7.下列函数中不能用于隐函数绘图的是()。D A.ezmeshB.ezsurfC.ezplotD.plot3 8.下列程序运行后,看到的图形()。C t=0:pi/20:2*pi; [x,y]=meshgrid(-8:0.5:8); z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);

matlab课后习题答案(1-9章)

1 数字1.5e2,1.5e3 中的哪个与1500相同吗?1.5e3 2 请指出如下5个变量名中,哪些是合法的? abcd-2 xyz_3 3chan a 变量 ABCDefgh 2、5是合法的。 3 在MATLAB 环境中,比1大的最小数是多少? 1+eps 4 设 a = -8 , 运行以下三条指令,问运行结果相同吗?为什么? w1=a^(2/3) w2=(a^2)^(1/3) w3=(a^(1/3))^2 w1 = -2.0000 + 3.4641i ;w2 = 4.0000 ;w3 =-2.0000 + 3.4641i 5 指令clear, clf, clc 各有什么用处? clear 清除工作空间中所有的变量。 clf 清除当前图形。clc 清除命令窗口中所有显示。 第二章 1 说出以下四条指令产生的结果各属于哪种数据类型,是“双精度”对象,还是“符号”符号对象? 3/7+0.1双; sym(3/7+0.1)符; sym('3/7+0.1') 符;; vpa(sym(3/7+0.1)) 符; 2 在不加专门指定的情况下,以下符号表达式中的哪一个变量被认为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') symvar(sym('sin(w*t)'),1) w a z 3 (1)试写出求三阶方程05.443 =-x 正实根的程序。注意:只要正实根,不要出现其他根。 (2)试求二阶方程022=+-a ax x 在0>a 时的根。 (1)reset(symengine) syms x positive solve(x^3-44.5) ans = (2^(2/3)*89^(1/3))/2 (2)求五阶方程02 2 =+-a ax x 的实根 syms a positive %注意:关于x 的假设没有去除 solve(x^2-a*x+a^2) Warning: Explicit solution could not be found. > In solve at 83 ans = [ empty sym ] syms x clear syms a positive solve(x^2-a*x+a^2) ans = a/2 + (3^(1/2)*a*i)/2 a/2 - (3^(1/2)*a*i)/2 4 观察一个数(在此用@记述)在以下四条不同指令作用下的异同。 a =@, b = sym( @ ), c = sym( @ ,' d ' ), d = sym( '@ ' ) 在此,@ 分别代表具体数值 7/3 , pi/3 , pi*3^(1/3) ;而异同通过vpa(abs(a-d)) , vpa(abs(b-d)) , vpa(abs(c-d))等来观察。 ● 理解准确符号数值的创建法。 ● 高精度误差的观察。 (1)x=7/3 x=7/3;a=x,b=sym(x),c=sym(x,'d'),d=sym('7/3'), a = 2.3333 b = 7/3 c =

MATLAb与数学实验 第五章习题解答

第五章习题 1.创建矩阵A= 100 021 013 ⎡⎤ ⎢⎥ - ⎢⎥ ⎢⎥ - ⎣⎦ ,B= 357 010 ⎡⎤ ⎢⎥ ⎣⎦ ,C= 11 11 11 ⎡⎤ ⎢⎥ ⎢⎥ ⎢⎥ ⎣⎦ ,D= 000 000 ⎡⎤ ⎢⎥ ⎣⎦ 答 A=[1 0 0;0 2 -1 ;0 -1 3] A = 1 0 0 0 2 -1 0 -1 3 >> C=[1 1;1 1;1 1] C = 1 1 1 1 1 1 >> B=[3 5 7;0 1 0] B = 3 5 7 0 1 0 >> D=[0 0 0;0 0 0] D = 0 0 0 0 0 0 2.随机生成:(1)一个含有五个元素的列向量. (2)一个数值在0~100之间的三行四列的矩阵. 答 (1) rand(5,1) ans = 0.9501

0.2311 0.6068 0.4860 0.8913 (2) 0+(100-0)*rand(3,4) ans = 76.2097 82.1407 79.1937 17.6266 45.6468 44.4703 92.1813 40.5706 1.8504 61.5432 73.8207 93.5470 3,生成一个5阶魔方矩阵 答 a=magic(5) a = 17 24 1 8 15 23 5 7 14 16 4 6 13 20 22 10 12 19 21 3 11 18 25 2 9 4.生成如下三对角矩阵A= 21000 22300 04250 00627 00082 -⎡⎤⎢⎥- ⎢⎥⎢⎥ - ⎢⎥ - ⎢⎥⎢⎥ -⎣⎦ 答 a1=-2*ones(1,5) a2=[1 3 5 7] a3=[2 4 6 8] A=diag(a1,0)+diag(a2,1)+diag(a3,-1) a1 = -2 -2 -2 -2 -2 a2 = 1 3 5 7

matlab课后习题答案1到6章

欢迎共阅 习题二 1. 如何理解“矩阵是MATLAB 最基本的数据对象”? 答:因为向量可以看成是仅有一行或一列的矩阵,单个数据(标量)可以看成是仅含一个元素的矩阵,故向量和单个数据都可以作为矩阵的特例来处理。 (2) 删除矩阵A 的第7号元 素。 答:A(7)=[] (3) 将矩阵A 的每个元素值 加30。 答:A=A+30; (4) 求矩阵A 的大小和维 数。 答:size(A); ndims(A); (5) 将向量 t 的0元素用 机器零来代替。 答:t(find(t==0))=eps; (6) 将含有12个元素的向 量 x 转换成34⨯矩阵。 答:reshape(x,3,4); (7) 求一个字符串的ASCII 码。 答:abs(‘123’); 或E 。 答:B=A(1:3,:); C=A(:,1:2); D=A(2:4,3:4); E=B*C; (2) 分别求E=10&A<25)。 答:E

B{1,2}='Brenden'; B{2,1}=reshape(1:9,3,3); B{2,2}={12,34,2;54,21,3;4 ,23,67}; (1)size(B)和ndims(B)的 值分别是多少? 答:size(B) 的值为2, 2。 ndims(B) 的值为 2。 (2)B(2)和B(4)的值分别是 (2)建立5×6随机矩阵A, 其元素为[100,200] 范围内的随机整数。答: round(100+(200-100)*ra nd(5,6)); (3)产生均值为1,方差为 0.2的500个正态分 布的随机数。 答: 1+sqrt(0.2)*randn(5,100 ); (4)产生和A同样大小的 幺矩阵。 答:ones(size(A)); (5)将矩阵A对角线的元 素加30。 答:A+eye(size(A))*30; (6)从矩阵A提取主对角 方程组的解为 x= 6.0000 26.6667 27.3333 -⎡⎤ ⎢⎥ ⎢⎥ ⎢⎥ ⎣⎦ 5.求下列矩阵的主对角线元素、 上三角阵、下三角阵、秩、范 数、条件数和迹。 (1) 1123 5142 3052 111509 A - ⎡⎤ ⎢⎥ - ⎢⎥ = ⎢⎥ ⎢⎥ ⎣⎦ (2) 0.43432 B ⎡⎤ =⎢⎥ ⎦ 答:[V,D]=eig(A); 习题四 1.从键盘输入一个4位整数,按 如下规则加密后输出。加密规 则:每位数字都加上7,然后 用和除以10的余数取代该数 字;再把第一位与第三位交 换,第二位与第四位交换。 答:

第5章 MATLAB绘图_习题答案

第5章 MATLAB绘图 习题5 一、选择题 1.如果x、y均为4×3矩阵,则执行plot(x,y)命令后在图形窗口中绘制()条曲线。D A.12 B.7 C.4 D.3 2.下列程序的运行结果是()。A x=0:pi/100:2*pi; for n=1:2:10 plot(n*sin(x),n*cos(x)) hold on end axis square A.5个同心圆B.5根平行线 C.一根正弦曲线和一根余弦曲线D.5根正弦曲线和5根余弦曲线 3.命令text(1,1,'{\alpha}+{\beta}')执行后,得到的标注效果是()。C A.{\alpha}+{\beta} B.{\α}+{\β} C.α+βD.\α+\β 4.subplot(2,2,3)是指()的子图。A A.两行两列的左下图B.两行两列的右下图 C.两行两列的左上图D.两行两列的右上图 5.要使函数y=2e x的曲线绘制成直线,应采用的绘图函数是()。C A.polar B.semilogx C.semilogy D.loglog 6.下列程序的运行结果是()。B [x,y]=meshgrid(1:5); surf(x,y,5*ones(size(x))); A.z=x+y平面B.与xy平面平行的平面 C.与xy平面垂直的平面D.z=5x平面 7.下列函数中不能用于隐函数绘图的是()。D A.ezmesh B.ezsurf C.ezplot D.plot3 8.下列程序运行后,看到的图形()。C t=0:pi/20:2*pi; [x,y]=meshgrid(-8:0.5:8); z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);

MATLAB语言与控制系统仿真-参考答案-第5章

控制系统的时域响应MATLAB 仿真实训 实训目的 1. 学会利用MATLAB 绘制系统的单位阶跃响应曲线,掌握读取系统动态性能指标的方法; 2. 学会利用MATLAB 绘制系统的单位脉冲响应曲线的方法; 3. 掌握利用MATLAB 绘制系统的零输入响应曲线的方法; 4. 掌握利用MATLAB 绘制系统的一般输入响应曲线的方法; 5. 学会通过仿真曲线读取相关信息,并依据有关信息进行系统的时域分析。 实训内容 1.编写程序求取下列各系统的单位阶跃响应,完成表5-5并记录相关曲线。 162.316)(21++= s s s G 16 4.216 )(22 ++=s s s G 166.116)(2 3++=s s s G 1616 )(24++=s s s G 解:>> n1=16; >> d1=[1,,16]; >> sys1=tf(n1,d1); >> step(sys1) >> n2=16; >> d2=[1,,16]; >> sys2=tf(n2,d2); >> step(sys2)

>> n3=16; >> d3=[1,,16]; >> sys3=tf(n3,d3); >> step(sys3) >> n4=16; >> d4=[1,1,16]; >> sys4=tf(n4,d4); >> step(sys4)

表5-5 序号 ξ n ω m ax c p t s t (%5=?) 计算值 实验 计算值 实验 计算值 实验值 1 4 2 4 3 4 4 4 w=4; cmax1=1+exp(-z1*pi/sqrt(1-z1^2)); tp1=pi/(w*sqrt(1-z1^2)); ts1=(z1*w); [cmax1,tp1,ts1] ans = >> z2=;

第5章 MATLAB绘图_习题答案

第5章MATLAB绘图 习题5 一、选择题 1.如果x、y均为4×3矩阵,则执行plot(x,y)命令后在图形窗口中绘制( )条曲线。D A。12 B。7 C.4 D.3 2.下列程序得运行结果就是()。A x=0:pi/100:2*pi; for n=1:2:10 plot(n*sin(x),n*cos(x)) hold on end axissquare A。5个同心圆B.5根平行线 C.一根正弦曲线与一根余弦曲线 D.5根正弦曲线与5根余弦曲线3.命令text(1,1,'{\alpha}+{\beta}’)执行后,得到得标注效果就是( )。C A.{\alpha}+{\beta}B。{\α}+{\β} C.α+βD.\α+\β 4。subplot(2,2,3)就是指( )得子图。A A.两行两列得左下图B.两行两列得右下图 C.两行两列得左上图D.两行两列得右上图 5。要使函数y=2ex得曲线绘制成直线,应采用得绘图函数就是().C A.polarB。semilogxC。semilogy D。loglog 6.下列程序得运行结果就是( )。B [x,y]=meshgrid(1:5); surf(x,y,5*ones(size(x))); A.z=x+y平面 B.与xy平面平行得平面 C。与xy平面垂直得平面 D.z=5x平面 7.下列函数中不能用于隐函数绘图得就是()。D A.ezmesh B.ezsurf C.ezplotD.plot3 8.下列程序运行后,瞧到得图形().C

t=0:pi/20:2*pi; [x,y]=meshgrid(-8:0、5:8); z=sin(sqrt(x、^2+y、^2))、/sqrt(x、^2+y、^2+eps); surf(x,y,z) view(0,90);axis equal A.像墨西哥帽子 B.就是空心得圆 C。边界就是正方形D.就是实心得圆 9。下列程序运行后得到得图形就是( ).A [x,y]=meshgrid(-2:2); z=x+y; i=find(abs(x)<1 & abs(y)<1); z(i)=NaN; surf(x,y,z);shadinginterp A。在一个正方形得正中心挖掉了一个小得正方形 B.在一个正方形得正中心挖掉了一个小得长方形 C.在一个正方形得上端挖掉了一个小得正方形 D。在一个正方形得下端挖掉了一个小得正方形 10.在使用MA TLAB“绘图”选项卡中得命令按钮绘图之前,需要().B A.在命令行窗口中输入绘图命令B。在工作区中选择绘图变量 C。打开绘图窗口D.建立M文件 二、填空题 1.执行以下命令: x=0:pi/20:pi; y=sin(x); 以x为横坐标、y为纵坐标得曲线图绘制命令为,给该图形加上“正弦波”标题得命令为,给该图形得横坐标标注为“时间",纵坐标标注为“幅度”得命令分别为与.plot(x,y),title(‘正弦波'),xlabel(‘时间’),ylabel(‘幅度') 2。在同一图形窗口中绘制y1与y2两条曲线,并对图形坐标轴进行控制,请补充程序. x=-3:0、1:3; y1=2*x+5; y2=x、^2-3*x; plot(x,y1)%绘制曲线y1 ①; plot(x,y2) %绘制曲线y2 m1=max([y1,y2]); m2=min([y1,y2]);

MATLAB习题及参考答案

习题: 1, 计算⎥⎦ ⎤⎢⎣⎡=572396a 与⎥⎦⎤ ⎢⎣⎡=864142b 的数组乘积。 2, 对于B AX =,如果⎥⎥ ⎥⎦ ⎤ ⎢⎢⎢⎣⎡=753467294A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ,求解X 。 3, 已知:⎥⎥ ⎥⎦ ⎤ ⎢⎢⎢⎣⎡=987654321a ,分别计算a 的数组平方和矩阵平方,并观察其结果。 4, 角度[]604530=x ,求x 的正弦、余弦、正切和余切。(应用sin,cos,tan.cot) 5, 将矩阵⎥⎦ ⎤⎢ ⎣⎡=7524a 、⎥⎦⎤⎢⎣⎡=3817b 和⎥⎦⎤ ⎢⎣⎡=2695c 组合成两个新矩阵: (1)组合成一个4⨯3的矩阵,第一列为按列顺序排列的a 矩阵元素,第二列为按列顺序排列的b 矩阵元素,第三列为按列顺序排列的c 矩阵元素,即 ⎥⎥ ⎥⎥⎦ ⎤⎢⎢⎢ ⎢⎣⎡237 912685 574 (2)按照a 、b 、c 的列顺序组合成一个行矢量,即 []296531877254 6, 将(x -6)(x -3)(x -8)展开为系数多项式的形式。(应用poly,polyvalm) 7, 求解多项式x 3-7x 2+2x +40的根。(应用roots) 8, 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。(应用poly,polyvalm) 9, 计算多项式9514124234++--x x x x 的微分和积分。(应用polyder,polyint ,poly2sym) 10, 解方程组⎥⎥⎥⎦ ⎤ ⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡66136221143092x 。(应用x=a\b) 11, 求欠定方程组⎥⎦ ⎤ ⎢⎣⎡=⎥ ⎦⎤⎢ ⎣⎡5865394742x 的最小范数解。(应用pinv)

matlab期末大作业题目及答案

matlab 期末大作业(30分,每题6分) 1. 积分运算(第四数值和五章符号) (1)定积分运算:分别采用数值法(quad ,dblquad )和符号运算(syms, int ) 一重定积分 π⎰ 1. 数值法(quad ) a) 运行代码: b) 运行结果: 2. 符号运算(syms ) a) 运行代码: b) 运行结果:

二重定积分11 22 00 () x y dxdy + ⎰⎰ 1.数值法(dblquad): a)运行代码: b)运行结果: 2.符号运算(syms): a)运行代码: b)运行结果:

(2) 不定积分运算 sin dxdy ⎰⎰ ((x/a)+b/y) i. 运行代码: ii. 运行结果: 2. 用符号法和数值法求解线性代数方程 (第五章和第二章) ⎩⎨ ⎧=+=+12 *22x *213 *12x *a11y a a y a (1) 用syms 定义待解符号变量x,y 和符号参数a11,a12,a21,a22,用符号 solve 求x,y 通解 1. 运行代码: 2. 运行结果:

(2) 用subs 带入a11=2,a12=4,a21=6,a22=8,求x 和y 特解,用vpa 输出有 效数值4位的结果 1. 运行代码: 2. 运行结果: (3) 采用左除(\)和逆乘法求解符号参数赋值后的方程 ⎩⎨ ⎧=+=+12*8x *63*4x *2y y 1. 运行代码: 2. 运行结果:

3.数值法和符号法求解非线性方程组(第四数值和五章符号 ) (1)采用数值法(fsolve )求解初始估计值为x0 = [-5; -5]的数值解 1. 运行代码: 2. 运行结果: 2 1x 21x 21e x 2x e x x 2--=+-= -

matlab第五章部分课后答案

第五章1. (1) x=-10:0.1:10; y=100./(1+x.^2); plot(x,y) (2) x=-10:0.1:10; y=1/(2*pi)*exp(-x.^2/2); plot(x,y) (3)ezplot('x^2+y^2=1') (4) t=-10:0.1:10; x=t.^2; y=5*t.^3; plot(x,y) 2. (1) theta=0:0.01:2*pi; rho=5*cos(theta)+4; polar(theta,rho) (2) theta=0.001:0.1:2*pi; rho=12./(sqrt(theta); polar(theta,rho) (3)theta=0:0.1:2*pi; rho=5./(cos(theta)+eps)-7; polar(theta,rho) (4) theta=0:0.1:2*pi;

rho=pi/3.*theta.^2; polar(theta,rho) 3. (1) t=0:pi/100:2*pi; x=cos(t); y=sin(t); z=t; plot3(x,y,z) (2) u=0:pi/100:2*pi; v=0:pi/100:2*pi; x=(1+cos(u)).*cos(v); y=(1+cos(u)).*sin(v); z=sin(u); plot3(x,y,z) (3) z=ones(50,50)*5; surf(z); (4) [x,y,z]=sphere(); surf(10*x,10*y,10*z); 4. clc clear all subplot(2,2,1) t=-pi/4:pi/2:(2*pi-pi/4); x=cos(t); y=sin(t); plot(x,y);

MATLAB习题及参考答案

2 9 0 13 10 , 解方程组3 4 11 x 6。(应用x=a\b) 2 2 6 6 11,求欠定方程组2 9 4 7 3 5 4 x 6 8 的最小范数解。(应用 pinv) 5 习题: 6 9 3 2 4 1 与b 的数组乘积。 2 7 5 4 6 8 1,计算a 2, 对于AX B,如果A 37 26,求解X。 28 3,已知:a 1 2 3 4 5 6,分别计算a的数组平方和矩阵平方,并观察其结果。 7 8 9 4,角度x 30 45 60,求x的正弦、余弦、正切和余切。(应用sin,cos, 4 2 7 1 5 9 5,将矩阵a 、b 和c 组合成两个新矩阵: 5 7 8 3 6 2 (1)组合成一个4 3的矩阵,第一列为按列顺序排列的a矩阵元素,第二列为按列顺序排列的元素,第三列为按列顺序排列的c矩阵元素,即b矩阵 4 7 5 5 8 6 2 1 9 7 3 2 (2)按照a、b、c的列顺序组合成一个行矢量,即 452778135692 6,将(x-6)(x-3)(x-8)展开为系数多项式的形式。(应用 poly,polyvalm) 7,求解多项式X3-7X2+2X+40的根。(应用roots) 8,求解在x=8 时多项式(x-1)( x-2) ( x-3)( x-4)的值。(应用poly,polyvalm) 9, 计算多项式4x412x314x25x 9的微分和积分。(应用polyder,polyint ,poly2sym)

2 2 计算表达式z 10 x 3 y 5 e x y 的梯度并绘图。(应用meshgrid, gradient, con tour, hold on, quiver) 15,用符号函数法求解方程 a t 2+b*t +c=0。(应用solve ) 16,用符号计算验证三角等式:(应用syms,simple ) 用 syms,ezplot) 用 plot,title,text,legend) 24, x= [66 49 71 56 38] ,绘制饼图,并将第五个切块分离出来。 25, 用 sphere 函数产生球表面坐标,绘制不通明网线图、透明网线图、表面图和带剪孔的表面图。 ( 应 用 sphere, mesh, hidden off, surf, NaN ) 17, 求矩阵A an 312 的行列式值、逆和特征根。 (应用syms,det,inv,eig ) a 21 a 22 18, 因式分解: 4 x 5x 3 5x 2 5x 6 (应用 syms, factor) 19, f a 2 x 1 x ,用符号微分求 df/dx 。(应用 syms,diff) ax e log(x) sin (x) 20, 符号函数绘图法绘制函数 x=s in( 3t)cos(t) ,y=sin (3t )sin (t ) 的图形,t 的变化范围为[0,2 ]。(应 23,表中列出了 4个观测点的6次测量数据,将数据绘制成为分组形式和堆叠形式的条形图。 12,矩阵a 4 2 6 7 5 4 ,计算a 的行列式和逆矩阵。 3 49 (应用 det,inv) 13, y =sin( x ) , x 从 0 到 2 x =,求y 的最大值、最小值、均值和标准差。(应用max,min,mean,std ) 14,参照课件中例题的方法, 21, 绘制曲线y 3 x x 1,x 的取值范围为[-5,5] 。(应用 plot ) 22,有一组测量数据满足 y e -at ,t 的变化范围为0~10,用不同的线型和标记点画出 a=、 a=禾口 a=三种 情况下的曲线,在图中添加标题y e -at ,并用箭头线标识出各曲线 a 的取值,并添加标题y e -at 和图例框。(应

MATLAB语言基础与应用(第二版)第5章 习题答案

第5章习题与答案 5.1用矩阵三角分解方法解方程组 1231231 23214453186920 x x x x x x x x x +-=⎧⎪ -+=⎨⎪+-=⎩ 解答: >>A=[2 1 -1;4 -1 3;6 9 -1] A = 2 1 -1 4 -1 3 6 9 -1 >>b=[1 4 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 分解方法解方程组 123121 33235223 3127 x x x x x x x ++=⎧⎪ +=⎨⎪+=⎩ 解答: >> A=[3 2 3;2 2 0;3 0 12] A =

3 2 3 2 2 0 3 0 12 >> b=[5;3;7] b = 5 3 7 >> L=chol(A) L = 1.7321 1.1547 1.7321 0 0.8165 -2.4495 0 0 1.7321 >> y=backsub(L,b) y = -11.6871 15.7986 4.0415 >> x=backsub(L',y) x = -6.7475 28.8917 49.9399 5.3 解答: 观察数据点图形 >> x=0:0.5:2.5 x = 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)

MATLAB教程课后实验报告题目及解答[第一至第五章]

MATLAB教程实验报告 实验项目名称实验一 Matlab基本操作 学生姓名汪德旺专业班级 09数教(1)班学号 0301090131 实验成绩日期 一. 实验目的和要求 1、了解MATLAB 的开发环境。 2、熟悉Matlab的基本操作。 3、掌握建立矩阵的方法。 4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。 5、填写实验报告,实验报告文件取名为report1.doc。 6、于邮件附件形式将实验报告文件report1.doc 发到邮箱*******************,邮件主题为班级学号姓名,如:09数教1班15号张三。 二、实验内容 1、先求下列表达式的值,然后显示MATLAB工作空间的使用情况并保存全部变量。 (1) 0 12 2sin85 1 z e = + (2 ) 2 212 1 ln( 0.455 2 i z x x + ⎡⎤==⎢⎥ -⎣⎦ 其中 (3) 0.30.3 3 0.3 sin(0.3)ln, 22 a a e e a z a - -+ =++a=-3.0,-2.9,-2.8,…, 2.8, 2.9,3.0 (4) 2 2 4 2 ,01 1,12,0:0.5:2.5 21,23 t t z t t t t t ⎧≤< ⎪ =-≤<= ⎨ ⎪-+≤< ⎩ 其中t 2.已知: 12344131 34787,203 3657327 A B --⎡⎤⎡⎤ ⎢⎥⎢⎥== ⎢⎥⎢⎥ ⎢⎥⎢⎥ - ⎣⎦⎣⎦ 求下列表达式的值: (1)A+6*B和A-B+I(其中I为单位矩阵) (2)A*B和A.*B (3)A^3和A.^3 (4)A/B和B\A (5)[A,B]和[A([1,3],:);B^2]

matlab习题解答

上机练习题一 班级: 姓名: 学号: 1.建立起始值=3,增量值=,终止值=44的一维数组x 答案: x=(3::44) 2.写出计算 Sin(30o )的程序语句. 答案: sin(pi*30/180) 或 sin(pi/6) 3.矩阵 ⎥⎥ ⎥⎦ ⎤ ⎢⎢⎢⎣⎡=187624323A ,矩阵 ⎥⎥ ⎥⎦ ⎤ ⎢⎢⎢⎣⎡=333222111B ;分别求出B A ⨯及A 与 B 中对应元素之间的乘积. 答案:A = [3,2,3; 4,2,6; 7,8,1] B = [1,1,1; 2,2,2; 3,3,3] A*B ;A.*B 4计算行列式的值1 876 24 323=A 。答案:det(A)

5对矩阵 ⎥⎥ ⎥⎦ ⎤ ⎢⎢⎢⎣⎡=187624323A 进行下述操作。 (1)求秩。答案:rank(A) (2)求转置。答案:A' (3) 对矩阵求逆,求伪逆。答案:inv(A) ,pinv(A) (4) 左右反转,上下反转。答案:fliplr(A),flipud(A) (5) 求矩阵的特征值. 答案:[u,v]=eig(A) (6) 取出上三角和下三角. 答案:triu(a) tril(a) (7)以A 为分块作一个3行2列的分块矩阵。答案:repmat(a) 6 计算矩阵 ⎥⎥⎥⎦ ⎤⎢⎢⎢⎣⎡897473535与 ⎥⎥⎥⎦ ⎤ ⎢⎢⎢⎣⎡638976242之和。 >> a=[5 3 5;3 7 4;7 9 8]; >> b=[2 4 2;6 7 9;8 3 6]; >> a+b 7 计算⎥⎦ ⎤⎢⎣⎡=572 396 a 与⎥⎦ ⎤⎢⎣⎡=864 142 b 的数组乘积。 >> a=[6 9 3;2 7 5]; >> b=[2 4 1;4 6 8]; >> a.*b ans = 12 36 3

相关主题
相关文档
最新文档