MATLAB仿真 课后习题

第一章习题

3.请指出以下的变量名（函数名、M文件名）中，哪些是合法的？Abc 2004x lil-1 wu_2004 a&b qst.u _xyz 解：合法的变量名有：Abc wu_2004

4．指令窗操作

（1）求[12+2×（7-4）]÷32的运算结果

解：>> [12+2*(7-4)]/3^2

ans =

2

（2）输入矩阵A=[1，2，3；4，5，6；7，8，9]，观察输出。

解：>> A=[1,2,3;4,5,6;7,8,9]

A =

1 2 3

4 5 6

7 8 9

（3）输入以下指令，观察运算结果；

clear;x=-8:0.5:8;

y=x';

X=ones(size(y))*x;

Y=y*ones(size(x));

R=sqrt(X.^2+Y.^2)+eps;

Z=sin(R)./R;

mesh(X,Y,Z);

colormap(hot)

xlabel('x'),ylabel('y'),zlabel('z')

解：

7．指令行编辑

（1）依次键入以下字符并运行：y1=2*sin(0.3*pi)/(1+sqrt(5))

解：>>y1=2*sin(0.3*pi)/(1+sqrt(5))

y1 =

0.5000

(2)通过反复按键盘的箭头键，实现指令回调和编辑，进行新的计算；y2=2*cos(0.3*pi)/(1+sqrt(5))

解：>>y2=2*cos(0.3*pi)/(1+sqrt(5))

y2 =

0.3633

11.编写题4中（3）的M脚本文件，并运行之。

解：

第二章习题

1.在指令窗中键入x=1:0.2:2和y=2:0.2:1,观察所生成的数组。 解:>> x=1:0.2:2 x =

1.0000 1.2000 1.4000 1.6000 1.8000

2.0000 >> y=2:0.2:1 y =

Empty matrix: 1-by-0

2．要求在[0，2π]上产生50个等距采样数据的一维数组，试用两种不同的指令实现。

解： y1=0:2*pi/49:2*pi

y2=linspace(0,2*pi,50)

3.计算e -2t sint ，其中t 为[0，2π]上生成的10个等距采样的数组。

解：>> t=linspace(0,2*pi,10); x=exp(-2*t).*sin(t) x =

0 0.1591 0.0603 0.0131 0.0013 -0.0003 -0.0002 -0.0001 -0.0000 -0.0000

4.已知A=??????4321 , B=??

?

???8765,计算矩阵A 、B 乘积和点乘. 解：>> A=[1,2;3,4]; B=[5,6;7,8]; x=A*B x =

19 22 43 50 >> x=A.*B x =

5 12

21 32

5.已知A=??????05314320，B=??

?

???05314320，计算A&B, A|B, ~A, A==B, A>B. 解：>> A=[0,2,3,4;1,3,5,0]; B=[1,0,5,3;1,5,0,5]; a1=A&B

a2=A|B

a3=~A

a4=(A==B)

a5=(A>B)

a1 =

0 0 1 1

1 1 0 0

a2 =

1 1 1 1

1 1 1 1

a3 =

1 0 0 0

0 0 0 1

a4 =

0 0 0 0

1 0 0 0

a5 =

0 1 0 1

0 0 1 0

7.将题5中的A阵用串转换函数转换为串B，再size指令查看A、B的结构，有何不同？

解：>> A=[0,2,3,4;1,3,5,0]

B=num2str(A)

size(A)

size(B)

A =

0 2 3 4

1 3 5 0

B =

0 2 3 4

1 3 5 0

ans =

2 4

ans =

2 10

第三章习题

1.已知系统的响应函数为

)sin(1

1)(θββ

ε+-

=-t e t y t ,其中

???

?

??-=-=εεθεβ22

1arctan ,1 ,要求用不同线型或颜色,在同一张图上绘制ε取值分别为0.2、

0.4、0.6、0.8时,系统在t ∈[0,18] 区间内的响应曲线,并要求用ε=0.2和 ε=0.8对他们相应的两条曲线进行文字标志。

解： clc

close all clear all t=0:0.02:18;

xi=[0.2,0.4,0.6,0.8]'; sxi=sqrt(1-xi.^2); sita=atan(sxi./xi);

y=1-exp(-xi*t).*sin(sxi*t+sita*ones(1,901))./(sxi*ones(1,901))

plot(t,y(1), 'r-', t,y(2), ' b*', t,y(3), ' g+', t,y(4), ' k.')

text(4.2,1.4,'\xi =0.2') text(3.8,0.9,'\xi=0.8')

2.用plot3、mesh 、surf 指令绘制

()()2

22

2111

y x y x z +++

+-=

三维图(x,y 范围自定)。 解：

clc;close all ;clear all ; x=-5:0.1:5;y=-5:0.1:5; [X,Y]=meshgrid(x,y); a=sqrt((1-X).^2+Y.^2); b=sqrt((1+X).^2+Y.^2); Z=1./(a+b);

a1=sqrt((1-x).^2+y.^2); b1=sqrt((1+x).^2+y.^2); z=1./(a1+b1);

subplot(1,3,1),plot3(x,y,z),xlabel('x'),ylabel('y'),zlabel('z');box on ;

subplot(1,3,2),surf(X,Y,Z),xlabel('x'),ylabel('y'),zlabel('z');box

on ;

subplot(1,3,3),mesh(X,Y,Z),xlabel('x'),ylabel('y'),zlabel('z');box on ;

3.对向量t 进行以下运算可以构成三个坐标的值向量:x=sin(t),y=cos(t),z=t.利用指

令plot3,并选用绿色的实线绘制相应的三维曲线. 解：

t=(0:0.01:2)*pi;

x=sin(t);

y=cos(t);

z=t;

plot3(x,y,z,'b-');box on

1.请分别用for 和while 循环语句计算K=∑=63

02i i 的程序，再写出一种避免循环的

计算程序。(提示：可考虑利用MATLAB 的sum(X,n)函数，实现沿数组X 的第n 维求和。) 解：

1）K=0; for i=0:63; K=K+2^i; end K

K =1.8447e+019 2）i=0;K=0; while i<=63; K=K+2^i; i=i+1; end; K

K =1.8447e+019 3）i=0; X=0:63; for i=0:63;

X(i+1)=2^i; end

sum(X,2) ans =1.8447e+019

1.将下列系统的传递函数模型用MATLAB 语言表达出来。

(1))

170046842541254289()

1700109329135()(23452341+++++++++=s s s s s s s s s s G 解：

num=[1,35,291,1093,1700];

den=[1,289,254,2541,4684,1700]; sys=tf(num,den) (2))

15).(5).(1()

3(15)(2++++=

s s s s s G 解：

z=-3;

p=[-1,-5,-15]; k=15;

sys=zpk(z,p,k)

(3))

252).(1).(1()

23.()2.(.100)(23223+++-++++=

s s s s s s s s s s G

解：

z=[0,-2,-2]; p=[-1,1]; k=100;

sys1=zpk(z,p,k); num=[1,3,2]; den=[1,2,5,2]; sys2=tf(num,den); sys=series(sys1,sys2)

4.求题3中的系统模型的等效传递函数模型和零极点模型。 解：

A=[3,2,1;0,4,6;0,-3,-5]; B=[1,2,3]' ; C=[1,2,5]; D=0;

sys=ss(A,B,C,D); systf=tf(sys) syszpk=zpk(sys)

Transfer function: 20 s^2 - 83 s + 138 --------------------- s^3 - 2 s^2 - 5 s + 6 Zero/pole/gain:

20 (s^2 - 4.15s + 6.9) ----------------------- (s-3) (s-1) (s+2)

5.已知系统的动力学方程如下，试用MATLAB 语言写出它们的传递函数。 (1))(2)()(500)(50)(15)(.

..

.

..

)3(t r t r t y t y t y t y +=+++ 解：

num=[1,2,0];

den=[1,15,50,500]; sys=tf(num,den) Transfer function: s^2 + 2 s ------------------------- s^3 + 15 s^2 + 50 s + 500

(2) )(4)(4)(6)(3)(...

t r dt t y t y t y t y =+++? 解：

num=[4,0]; den=[1,3,6,4]; sys=tf(num,den) Transfer function: 4 s --------------------- s^3 + 3 s^2 + 6 s + 4

6.试用MATLAB 语言表示图5-13所示系统。当分别以y =x 2和f 为系统输出、输入时的传递函数模型和状态空间模型(图中k =7N/m,c 1=0.5N/m.s -1, c 2=0.2N/m.s -1,m 1=3.5kg, m 2=5.6kg)。

f

解：)(t

k=7;

c1=0.5;

c2=0.2;

m1=3.5;

m2=5.6;

num=[m1,c1,k];

den=[m1*m2,c1*m1+c2*m1+c1*m2,c1*c2+m2*k,c1*k+c2*k,0];

sys=tf(num,den)

Transfer function:

3.5 s^2 + 0.5 s + 7

--------------------------------------

19.6 s^4 + 5.25 s^3 + 39.3 s^2 + 4.9 s

7.试用MATLAB语言分别表示图5-14所示系统质量m1,m2的位移x1,x2对输入f 的传递函数X2(s)/F(s)和X1(s)/F(s),其中m1=12kg, m2=38kg,k=1000N/m, c=0.1N/m.s-1。

解：

m1=12;

m2=38;

k=1000;

c=0.1;

num=[c,k];

den=[m1*m2,m1*c+m2*c,m1*k+m2*k,0,0];

sys1=tf(num,den)

num=[m1,c,k];

den=[m1*m2,m1*c+m2*c,m1*k+m2*k,0,0];

sys2=tf(num,den)

Transfer function:

0.1 s + 1000

---------------------------

456 s^4 + 5 s^3 + 50000 s^2

Transfer function:

12 s^2 + 0.1 s + 1000

---------------------------

456 s^4 + 5 s^3 + 50000 s^2

补充题

求图示传递函数

sys1=tf([1,2],[1,3,4]);

sys2=tf([1,4,5] ,[1,6,7,8]);

sys3=tf([1,0],[1,2]);

sys4=tf([1],[1,3]);

sys5=parallel(sys3,sys4);

sys=feedback(sys1*sys2*sys5,1,-1)

结果

s^5 + 10 s^4 + 39 s^3 + 74 s^2 + 66 s + 20

-----------------------------------------------------------------

s^7 + 14 s^6 + 81 s^5 + 262 s^4 + 530 s^3 + 684 s^2 + 538 s + 212

第六章习题2.将例6-2中的微分方程改写为以下形式：

1

)0(

,0

)0(0

.) 1.(

..

2

..

=

==

+

-

-

y

y y

y

y

yμ

求μ分别为1、2时，在时间区间t=[0,20]微分方程的解。

解：

M函数文件

function dx=wffc(t,x,flag,ps)

dx=zeros(2,1);

dx(1)=x(2);

dx(2)=ps*(1-x(1)^2)*x(2)-x(1);

调用程序

clc;close all;clear all;

tspan=[0,20];

x0=[0,1];

ps=1;

[T1,X1]=ode45('wffc',tspan,x0,odeset,ps);

ps=2;

[T2,X2]=ode45('wffc',tspan,x0,odeset,ps);

plot(T1,X1(:,1),'r',T2,X2(:,1),'b-.')

X1(:,1)

X2(:,1)

3.对图6-18所示反馈系统进行单位阶跃响应和方波响应(方波周期为30s)仿真。要求：

(1)利用MATLAB模型连接函数求出系统闭环传递函数。

(2)利用step 函数求单位阶跃响应。

(3)利用gensig 函数产生方波信号，利用lsim 函数求方波响应。 解：

clc;close all ;clear all ; % (1)

sys1=tf([1,0.5],[1,0.1]); sys2=ZPK([],[0,-2,-10],20); sys3=series(sys1,sys2); sys4=feedback(sys3,1,-1); % (2)

subplot(1,2,1) step(sys4); % (3)

[u,t]=gensig('square',30,60); subplot(1,2,2) lsim(sys4,'r',u,t)

20 (s+0.5)

--------------------------------------------

(s+10.23) (s+0.8195) (s^2 + 1.052s + 1.193)

4.已知系统传递函数01

.12.01

)(2++=

s s s G ;

(1)绘制系统阶跃响应曲线。

(2)绘出离散化系统阶跃响应曲线，采样周期T s =0.3s 。 解：

clc;close all ;clear all ; % (1)

sys=tf([1],[1,0.2,1.01]); subplot(1,2,1)

step(sys) % (2)

sys=tf([1],[1,0.2,1.01]); sys1=c2d(sys,0.3,'zoh'); [num,den]=tfdata(sys1,'v'); subplot(1,2,2) dstep(num,den)

附加题

1、已知二阶微分方程0342=-+-y y y y y ，其初始条件为0)0(=y ，1)0(=y ，求在时间范围t=[0 5]内该微分方程的解。 M 函数为：

function dy=vdp(t,y)

dy=zeros(2,1); dy(1)= y(2); dy(2)= 4*y(2)-(y(1)^2)*y(2)+3*y(1); 调用函数为：

[T,Y]=ode45('vdp',[0 5],[0,1]); plot(T,Y(:,1),'r-',T,Y(:,2),'b:')

2、已知系统模型为

7

22)(3

+++=s s s s G ，计算系统在周期10s 的方波信号作用下5

个周期内的时间响应，并在同一图形窗口中绘制输入信号和时间响应曲线。 sys=tf([1,2],[1,0,2,7]);

[u,t]=gensig('square',10,50); %产生方波信号数据

lsim(sys,'r',u,t) , hold on %产生方波响应并绘曲线 plot(t,u,'-.') %在同一坐标系绘方波波形 hold off

第七章习题

1.绘制下列各单位反馈系统开环传递函数的Bode 图和Nyquist 图，并根据其稳定裕度判断系统的稳定性。 (1))

31).(21).(1(10

)(s s s s G k +++=

解：

clc;clear all ;close all ; % (1)

Gk=zpk([],[0,-0.5,-1/3],5/3); subplot(1,2,1) margin(Gk) grid on

subplot(1,2,2) nyquist(Gk)

由上图的稳定裕度知系统临界稳定。 (2))

101).(1.(10

)(s s s s G k ++=

解：

clc;clear all ;close all ; % (2)

Gk=zpk([],[0,-1,-0.1],1); subplot(1,2,1) margin(Gk) grid on

subplot(1,2,2) nyquist(Gk)

由上图的稳定裕度知系统不稳定。 (3))

2.01).(1.01.(10

)(2s s s s G k ++=

解：

clc;clear all ;close all ; % (3)

Gk=zpk([],[0,0,-10,-5],500); subplot(1,2,1) margin(Gk) grid on

subplot(1,2,2) nyquist(Gk)

由上图的稳定裕度知系统不稳定。 (4) )

101).(1.01.(2

)(2s s s s G k ++=

解：

clc;clear all;close all; % (4)

Gk=zpk([],[0,0,-10,-0.1],2); subplot(1,2,1) margin(Gk) grid on

subplot(1,2,2) nyquist(Gk)

由上图的稳定裕度知系统不稳定。

2.设单位反馈系统的开环传递函数为

)

12.()(22++=

n n

k w s w s s K

s G ξ,其中无阻尼固有频率w n =90rad/s ，阻尼比ξ=0.2,试确定使系统稳定的K 的范围。 解： 方法1

g=tf(1,[1/90^2 0.4/90 1 0]);%系统开环模型

w=logspace(0,3,1000); %生成频率向量 bode(g,w)

[mag,phase,w]=bode(g,w); %产生幅值（非分贝）和相位向量

mag1=reshape(mag,1000,1); %重构幅值向量（1000*1）

phase1=reshape(phase,1000,1);%重构相频向量（1000*1）

wc=interp1(phase1,w,-180) %插值求-180度所对应的频率——wc

gk=interp1(w,mag1,wc) %插值求wc所对应的增益

gkk=1/gk %该增益的倒数即为可增加的最大增益

wc =

90.0004

gk =

0.0278

gkk =

36.0033

方法2

wc=0;wg=0.01;k=1;

while wc

sys=tf(k,[1/(90*90),2*0.2/90,1,0]);

[gm,pn,wg,wc]=margin(sys);

k=k+0.1;

end

k-0.1

ans =

36.0000

方法3

xi=0.2;omega=90;w=90;

sys1=tf(1,[1,0]);

sys2=tf(1,[1/w^2,2*xi/w,1]);

sys=series(sys1,sys2);

[Gm,Pm,Wcg,Wcp]=margin(sys);

k=Gm

k =

36

3.设系统结构如图7-22所示，试用LTI Viewer分析系统的稳定性,并求出系统的稳定裕度及单位阶跃响应峰值。

clc;close all;clear all;

G11=0.5;

G12=zpk([0],[-0.5],1);

G1=G11-G12;

G2=tf(1,[1 2 0]);

Gk=G1*G2;

Gb=feedback(Gk,1,-1);

[Gm,Pm,Wcg,Wcp]=margin(Gb)

step(Gb)

[y,t]=step(Gb);

matlab基础练习题带答案

Matlab 基础练习题 常量、变量、表达式 1、 MATLAB 中，下面哪些变量名是合法的？（ ） （A ）_num （B ）num_ （C ）num- （D ）-num 2、 在MA TLAB 中，要给出一个复数z 的模，应该使用（ ）函数。 （A ）mod(z) （B ）abs(z) （C ）double(z) （D ）angle(z) 3、 下面属于MATLAB 的预定义特殊变量的是？（ ） （A ）eps （B ）none （C ）zero （D ）exp 4、 判断：在MA TLAB 的内存工作区中，存放一个英文字符 'a' 需要占用1个字节，存放 一个中文字符‘啊’需要占用2个字节。（ 错，都是2个字节 ） 5、 判断：MA TLAB 中，i 和j （ 对 ） 6、 判断：MA TLAB 中，pi 代表圆周率，它等于3.14。（ 错，后面还有很多位小数 ） 7、 在MA TLAB 中，若想计算的5 1)3.0sin(21+= πy 值，那么应该在MA TLAB 的指令窗中 输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。 8、 在MA TLAB 中，a = 1，b = i ，则a 占_8__个字节，b 占_16_个字节，c 占________字 节。 9、 在MA TLAB 中，inf 的含义是__无穷大__，nan 的含义是__非数（结果不定）___。 数组 1、 在MA TLAB 中，X 是一个一维数值数组，现在要把数组X 中的所有元素按原来次序 的逆序排列输出，应该使用下面的（ ）指令。 （A ）X[end:1] （B ）X[end:-1:1] （C ）X (end:-1:1) （D ）X(end:1) 2、 在MA TLAB 中，A 是一个字二维数组，要获取A 的行数和列数，应该使用的MATLAB 的命令是（ ）。 （A ）class(A) （B ）sizeof(A) （C ）size(A) （D ）isa(A) 3、 在MATLAB 中，用指令x=1:9生成数组x 。现在要把x 数组的第二和第七个元素都 赋值为0，应该在指令窗中输入（ ） （A ）x([2 7])=(0 0) （B ）x([2,7])=[0,0] （C ）x[(2,7)]=[0 0] （D ）x[(2 7)]=(0 0) 4、 在MA TLAB 中，依次执行以下指令：clear;A=ones(3,4); A(:)=[-6:5];这时， 若在指令窗中输入指令b=A(:,2)'，那么，MATLAB 输出的结果应该是（ ） （A ）b = -3 -2 -1 （B ）b = -2 -1 0 1 （C ）b = -5 -1 3 （D ）b = -5 -2 1 4 5、 在MA TLAB 中，A = 1:9，现在执行如下指令L1 = ~(A>5)，则MATLAB 的执行结果应 该是L1 =___ 1 1 1 1 1 0 0 0 0___。

控制系统的MATLAB仿真与设计课后答案

控制系统的MATLAB仿真与设计课后答案

>>z=-4*sqrt(2)*sin(t); >>plot3(x,y,z,'p'); >>title('Line in 3-D Space'); >>text(0,0,0,'origin'); >>xlabel('X'),ylable('Y'),zlable('Z');grid; 4>>theta=0:0.01:2*pi; >>rho=sin(2*theta).*cos(2*theta); >>polar(theta,rho,'k'); 5>>[x,y,z]=sphere(20); >>z1=z; >>z1(:,1:4)=NaN; >>c1=ones(size(z1)); >>surf(3*x,3*y,3*z1,c1); >>hold on >>z2=z; >>c2=2*ones(size(z2)); >>c2(:,1:4)=3*ones(size(c2(:,1:4))); >>surf(1.5*x,1.5*y,1.5*z2,c2); >>colormap([0,1,0;0.5,0,0;1,0,0]); >>grid on >>hold off 第四章 1>>for m=100:999 m1=fix(m/100); m2=rem(fix(m/10),10); m3=rem(m,10); if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end 2M文件：function[s,p]=fcircle(r) s=pi*r*r; p=2*pi*r; 主程序： [s,p]=fcircle(10) 3>>y=0;n=100; for i=1:n y=y+1/i/i; end >>y

matlab课后习题解答第二章

第2章符号运算 习题2及解答 1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度” 对象，还是“符号”符号对象 3/7+; sym(3/7+; sym('3/7+'); vpa(sym(3/7+) 〖目的〗 不能从显示形式判断数据类型，而必须依靠class指令。 〖解答〗 c1=3/7+ c2=sym(3/7+ c3=sym('3/7+') c4=vpa(sym(3/7+) Cs1=class(c1) Cs2=class(c2) Cs3=class(c3) Cs4=class(c4) c1 = c2 = 37/70 c3 = c4 = Cs1 = double Cs2 = sym Cs3 = sym Cs4 = sym 2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认 为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') 〖目的〗 理解自由符号变量的确认规则。 〖解答〗 symvar(sym('sin(w*t)'),1) ans = w symvar(sym('a*exp(-X)'),1)

ans = a symvar(sym('z*exp(j*th)'),1) ans = z 5求符号矩阵???? ??????=3332 31 232221 131211 a a a a a a a a a A 的行列式值和逆，所得结果应采用“子表达式置换”简洁化。 〖目的〗 理解subexpr 指令。 〖解答〗 A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]') DA=det(A) IA=inv(A); [IAs,d]=subexpr(IA,d) A = [ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA = a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs = [ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d = 1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31) 8（1）通过符号计算求t t y sin )(=的导数 dt dy 。（2）然后根据此结果，求- =0t dt dy 和2 π = t dt dy 。 〖目的〗 diff, limit 指令的应用。 如何理解运行结果。 〖解答〗 syms t

MATLAB仿真课后习题

第一章习题 3.请指出以下的变量名（函数名、M文件名）中，哪些是合法的？Abc 2004x lil-1 wu_2004 a&b qst.u _xyz 解：合法的变量名有：Abc wu_2004 4．指令窗操作 （1）求[12+2×（7-4）]÷32的运算结果 解：>> [12+2*(7-4)]/3^2 ans = 2 （2）输入矩阵A=[1，2，3；4，5，6；7，8，9]，观察输出。解：>> A=[1,2,3;4,5,6;7,8,9] A = 1 2 3 4 5 6 7 8 9 （3）输入以下指令，观察运算结果； clear;x=-8:0.5:8; y=x'; X=ones(size(y))*x; Y=y*ones(size(x)); R=sqrt(X.^2+Y.^2)+eps;

Z=sin(R)./R; mesh(X,Y,Z); colormap(hot) xlabel('x'),ylabel('y'),zlabel('z') 解： 7．指令行编辑 （1）依次键入以下字符并运行：y1=2*sin(0.3*pi)/(1+sqrt(5)) 解：>>y1=2*sin(0.3*pi)/(1+sqrt(5)) y1 = 0.5000 (2)通过反复按键盘的箭头键，实现指令回调和编辑，进行新的计算；y2=2*cos(0.3*pi)/(1+sqrt(5)) 解：>>y2=2*cos(0.3*pi)/(1+sqrt(5)) y2 =

0.3633 11.编写题4中（3）的M脚本文件，并运行之。解：

第二章习题 1.在指令窗中键入x=1:0.2:2和y=2:0.2:1,观察所生成的数组。 解:>> x=1:0.2:2 x = 1.0000 1.2000 1.4000 1.6000 1.8000 2.0000 >> y=2:0.2:1 y = Empty matrix: 1-by-0 2．要求在[0，2π]上产生50个等距采样数据的一维数组，试用两种不同的指令实现。 解： y1=0:2*pi/49:2*pi y2=linspace(0,2*pi,50) 3.计算e -2t sint ，其中t 为[0，2π]上生成的10个等距采样的数组。 解：>> t=linspace(0,2*pi,10); x=exp(-2*t).*sin(t) x = 0 0.1591 0.0603 0.0131 0.0013 -0.0003 -0.0002 -0.0001 -0.0000 -0.0000 4.已知A=??????4321 , B=?? ????8765,计算矩阵A 、B 乘积和点乘. 解：>> A=[1,2;3,4]; B=[5,6;7,8]; x=A*B x =

matlab课后答案解析完整版

ones表示1矩阵 zeros表示0矩阵 ones(4)表示4x4的1矩阵 zeros(4)表示4x4的0矩阵 zeros(4,5)表示4x5的矩阵 eye(10,10)表示10x10的单位矩阵 rand(4,5)表示4x5的伴随矩阵 det(a)表示计算a的行列式 inv(a)表示计算a的逆矩阵 Jordan(a)表示求a矩阵的约当标准块rank(a)表示求矩阵a的秩 [v,d]=eig(a)对角矩阵 b=a’表示求a矩阵的转置矩阵 sqrt表示求平方根 exp表示自然指数函数 log自然对数函数 abs绝对值 第一章 一、5（1） b=[97 67 34 10;-78 75 65 5;32 5 -23 -59]; >> c=[97 67;-78 75;32 5;0 -12]; >> d=[65 5;-23 -59;54 7]; >> e=b*c e = 5271 11574 -11336 664 1978 3112 （2）a=50:1:100 二、1 、x=-74; y=-27; z=(sin(x.^2+y.^2))/(sqrt(tan(abs(x+y)) )+pi) z = -0.0901 2、a=-3.0:0.1:3.0; >> b=exp(-0.3*a).*sin(a+0.3) 3、x=[2 4;-0.45 5]; y=log(x+sqrt(1+x.^2))/2

y = 0.7218 1.0474 -0.2180 1.1562 4、a*b表示a矩阵和b矩阵相乘 a.*b表示a矩阵和b矩阵单个元素相乘A(m,n)表示取a矩阵第m行，第n列 A(m,:)表示取a矩阵第m行的全部元素A(:,n)表示取a矩阵的第n列全部元素 A./B表示a矩阵除以b矩阵的对应元素， B.\A等价于A./B A.^B表示两个矩阵对应元素进行乘方运算 A.^2表示a中的每个元素的平方 A^2表示A*A 例：x=[1,2,3]; y=[4,5,6]; z=x.^y z= 1 3 2 729 指数可以是标量（如y=2）.底数也可以是标量（如x=2） 5、a=1+2i; >> b=3+4i; >> c=exp((pi*i)/6) c = 0.8660 + 0.5000i d=c+a*b/(a+b) d = 1.6353 + 1.8462i 第二章 二、4、（1） y=0;k=0; >> while y<3 k=k+1; y=y+1/(2*k-1); end >> display([k-1,y-1/(2*k-1)]) ans =

matlab例题

五、某公司投资2000万元建成一条生产线。投产后，在时刻t 的追加成本和追加收益分别为3/225)(t t t G ++=（百万元/年），3/218)(t t H -=（百万元/年）。试确定该生产线在何时停产可获最大利润？最大利润是多少？ 提示：利用函数?=T G H t R 0t 20-d ))t (-)t (()(（百万元），由于H （t ）-G （t ）单调 下降，所以H （t ）=G （t ）时，R （t ）取得最大利润。 5.解：构造函数f(t)=H(t)-G(t)=13-t-3t 2/3=0 ; 令t 1/3=x,则f(t)=-t 3-3t 2+13 可得矩阵P=[-1,-3,0,13] 求最佳生产时间的源程序如下： p=[-1,-3,0,13]; x=roots(p); t=x.^3 运行结果如下： t = 3.6768 +21.4316i 3.6768 -21.4316i 4.6465 再分别将t 的三个值带入函数f(t)，比较大小后，得到最大利润与最佳生产时间。 求最大利润的程序代码如下： ① t=3.6768 +21.4316i; x=0:0.01:t; y=13-x-3*x.^(2/3); trapz(x,y) 运行结果： ans = 25.2583 ② t=3.6768 -21.4316i; x=0:0.01:t; y=13-x-3*x.^(2/3); trapz(x,y) 运行结果： ans = 25.2583 ③ t=4.6465; x=0:0.01:t; y=13-x-3*x.^(2/3); trapz(x,y) 运行结果： ans = 26.3208 比较以上三组数据，可知最佳生产时间t=4.6465年，可获得的最大利润 26.3208（百万元/年）。 clear; close; fplot('18-t^(2/3)',[0,20]);grid on;hold on; fplot('5+t+2*t^(2/3)',[0,20],'r');hold off; %发现t 约为4

matlab程序设计第三章课后习题答案

1. p138 第6题在同一坐标轴中绘制下列两条曲线并标注两曲线交叉点。 >> t=0:0.01:pi; >> x1=t; >> y1=2*x1-0.5; >> x2=sin(3*t).*cos(t); >> y2=sin(3*t).*sin(t); >> plot(x1,y1,'r-',x2,y2,'g-') >> axis([-1,2,-1.5,1]) >> hold on >> s=solve('y=2*x-0.5','x=sin(3*t)*cos(t)','y=sin(3*t)*sin(t)'); >> plot(double(s.x),double(s.y),'*'); 截图：

p366 第4题绘制极坐标曲线，并分析对曲线形状的影响。 function [ output_args ] = Untitled2( input_args ) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here theta=0:0.01:2*pi; a=input('请输入a的值：'); b=input('请输入b的值：'); n=input('请输入n的值：'); rho=a*sin(b+n*theta); polar(theta,rho,'k'); end 下面以a=1，b=1，n=1的极坐标图形为基础来分析a、b、n的影响。

对a的值进行改变：对比发现a只影响半径值的整倍变化 对b的值进行改变：对比发现b的值使这个圆转换了一定的角度

对n的值进行改变：对比发现当n>=2时有如下规律 1、当n为整数时，图形变为2n个花瓣状的图形 2、当n为奇数时，图形变为n个花瓣状的图形 分别让n为2、3、4、5

matlab课后习题及答案详解

第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 在安装MATLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？

matlab典型例题

【例】水资源系统规划调度常应用系统分析方法处理，以一个水资源分配问题为例，讨论线性规划问题。例：有甲、乙两个水库同时给A、B、C三个城市供水，甲水库的日供水量为28万m3/d，乙水库的日供水量为35万m3/d，三个城市的日需水量分别为A≥10万m3/d，B≥15万m3/d，C≥20 万m3/d。由于水库与各城市的距离不等，输水方式不同，因此单位水费也不同。各单位水费分别为c11＝2000元/万m3、c12＝3000元/万m3、c13＝4000元/万m3、c21＝4500元/万m3、c22＝3500元/万m3、c23＝3000元/万m3。试作出在满足对三个城市供水的情况下，输水费用最小的方案。

设甲水库向三城市日供水量分别为x ll、x12、x13，乙水库向三城市日供水量分别为x2l、x22、x23。 建立约束条件： x11 + x21 ≥10 x12 + x22 ≥15 x13 + x23 ≥20 x11 + x12 + x13 ≤28 x21 + x22 + x23 ≤35 x11，x12，x13，x21，x22，x23，≥0 目标函数： fmin=c11x11+c12x12+c13x13+c21x21+c22x22+c23x23

这样的问题单纯求解是非常繁琐的，而MLTLAB求解是十分简单的，只要在命令行输入： 》f=[0.2 0.3 0.4 0.45 0.35 0.3]’； 》A=[-1 0 0 -1 0 0；0 -1 0 0 -1 0；0 0 -1 0 0 -1；1 1 1 0 0 0； 0 0 0 1 1 1]； 》B=[-10 -15 -20 28 35]； 》lb=zeros（6，1）； 》[X，Zmin]=linprog(f，A，B，[]，[]，lb,[]) 最后得出x=[10 15 0 0 0 20]万m3，Zmin=12.5万元。

Simulink 比值控制系统课后习题

Simulink 比值控制课后答案 作业题目： 在例一中如系统传递函数为43 ()151s G s e s -= +，其他参数不变，试对其进行单闭环比 值控制系统仿真分析，并讨论43 ()151 s G s e s -=+分母中“15”变化10%±时控制系统的鲁棒 性。 （1）分析从动量无调节器的开环系统稳定性。 由控制理论知，开环稳定性分析是系统校正的前提。系统稳定性的分析可利用Bode 图进行，编制MATLAB Bode 图绘制程序（M-dile ）如下： clear all close all T=15;K0=3;tao=4; num=[K0];den=[T,1]; G=tf(num,den,'inputdelay',tao); margin(G) 执行该程序得系统的Bode 图如图所示，可见系统是稳定的。幅值裕量为6.77dB ，对应增益为2.2。 -40-30-20-100 10M a g n i t u d e (d B )10 10101010 -2160 -1800-1440-1080-720 -3600P h a s e (d e g ) Bode Diagram Gm = 6.77 dB (at 0.431 rad/sec) , P m = 66.3 deg (at 0.189 rad/sec) Frequency (rad/sec) （2）选择从动量控制器形式及整定其参数。 根据工程整定的论述，选择PI 形式的控制器，即() I p K G s K s =+ 。本处采用稳定边界法整定系统。先让I K =0，调整p K 使系统等幅振荡（由稳定性分析图知在p K =2.2附近时系统震荡），即使系统处于临界稳定状态。

matlab习题

1、标点符号;可以使命令行不显示运算结果，%用来表示该行为注释行。 2、x为0～4pi，步长为0.1pi的向量，使用命令x=0:0.1*pi:4*pi创建。4、输入矩阵A= 错误！未找到引用源。 ，使用全下标方式用A(2,2)取出元素 “-5”，使用单下标方式用A(5)取出元素“-5”。 5、符号表达式sin(2*a+t)+m中独立的符号变量为t。 6、M脚本文件和M函数文件的主要区别是M脚本文件没有函数定义和M函数文件有函数定义。 7. 设x是一维数组，x的倒数第3个元素表示为 ______x(_end-2_)________________；要在x的第36和37个元素之间插入一个元素154，使用的命令(集)为_x=x(_1:36,[154],37:end)_；设y为二维数组，要删除y 的第34行和48列，可使用命令_y(34,:)=[];y(:,48)=[]_； 8. 将变量x以Ascii文本格式存储到文件fname.txt，应使用命令__save _x_；将Excel文件data.xls读入WorkSpace并赋值给变量x，可使用命令 _x=xlsread('data.xls')_； 9. 在while表达式,语句体,End 循环语句中，表达式的值__非零__时表示循环 条件为真，语句体将被执行，否则跳出该循环语句； 10. 打开Matlab的一个程序文件fname.m，以添加的方式进行读写，应当使用命 令_fid= fopen('fname.m','w+'); 11.要从键盘读入一个字符串并赋值给变量x，且给出提示“Who is she?”，应使用命令__x=input(…Who is she??,?s?)_；使用fprintf()函数以含有2位小数的6位浮点格式输出数值变量weight到文件wt.txt，使用的命令（集）为fid = fopen(…wt.txt?,??,_); fprintf () ；12．设A= 错误！未找到引用源。,和B= 错误！未找到引用源。 和 C=错误！未找到引用源。均为m*n矩阵，且存在于WorkSpace中，要产生矩阵D=

控制系统的MATLAB仿真和设计课后答案

第二章 1>>x=[15 22 33 94 85 77 60] >>x(6) >>x([1 3 5]) >>x(4:end) >>x(find(x>70)) 2>>T=[1 -2 3 -4 2 -3] ; >>n=length(T); >>TT=T'; >>for k=n-1:-1:0 >>B(:,n-k)=TT.^k; >>end >>B >>test=vander(T) 3>>A=zeros(2,5); >>A(:)=-4:5 >>L=abs(A)>3 >>islogical(L) >>X=A(L) 4>>A=[4,15,-45,10,6;56,0,17,-45,0] >>find(A>=10&A<=20) 5>>p1=conv([1,0,2],conv([1,4],[1,1]));

>>p2=[1 0 1 1]; >>[q,r]=deconv(p1,p2); >>cq='商多项式为'; cr='余多项式为'; >>disp([cq,poly2str(q,'s')]),disp([cr,poly2str(r,'s')]) 6>>A=[11 12 13;14 15 16;17 18 19]; >>PA=poly(A) >>PPA=poly2str(PA,'s') 第三章 1>>n=(-10:10)'; >>y=abs(n); >>plot(n,y,'r.','MarkerSize',20) >>axis equal >>grid on >>xlabel('n') 2>>x=0:pi/100:2*pi; >>y=2*exp(-0.5*x).*sin(2*pi*x); >>plot(x,y),grid on; 3>>t=0:pi/50:2*pi; >>x=8*cos(t); >>y=4*sqrt(2)*sin(t); >>z=-4*sqrt(2)*sin(t); >>plot3(x,y,z,'p');

matlab练习题和答案

matlab练习题和答案 控制系统仿真实验 Matlab 部分实验结果 目录 实验一MATLAB基本操作..................................................................... .......................1 实验二Matlab 编程..................................................................... ...............................5 实验三Matlab 底层图形控制..................................................................... .....................6 实验四控制系统古典分析...................................................................... .......................12 实验五控制系统现代分析..................................................................... . (15) 实验六PID 控制器的设计..................................................................... ......................19 实验七系统状态空间设计...................................................................... .......................23 实验九直流双闭环调速系统仿真..................................................................... . (25) 实验一MATLAB基本操作

最新MATLAB基础教程-薛山第二版-课后习题答案

《MATLAB及应用》实验指导书《MATLAB及应用》实验指导书 班级：T1243-7 姓名：柏元强 学号：20120430724 总评成绩： 汽车工程学院 电测与汽车数字应用中心

目录 实验04051001 MATLAB语言基础 (1) 实验04051002 MATLAB科学计算及绘图 (18) 实验04051003 MATLAB综合实例编程 (31)

实验04051001 MATLAB语言基础 1实验目的 1)熟悉MATLAB的运行环境 2)掌握MATLAB的矩阵和数组的运算 3)掌握MATLAB符号表达式的创建 4)熟悉符号方程的求解 2实验内容 第二章 1.创建double的变量，并进行计算。 (1)a=87，b=190，计算 a+b、a-b、a*b。 clear,clc a=double(87); b=double(190); a+b,a-b,a*b (2)创建 uint8 类型的变量，数值与(1)中相同，进行相同的计算。 clear,clc a=uint8(87); b=uint8(190); a+b,a-b,a*b 2．计算： (1) () sin 60

(2) e3 (3) 3cos 4??π ??? clear,clc a=sind(60) b=exp(3) c=cos(3*pi/4) 3．设2u =，3v =，计算： (1) 4 log uv v (2) () 2 2 e u v v u +- (3) clear,clc u=2;v=3; a=(4*u*v)/log(v) b=((exp(u)+v)^2)/(v^2-u) c=(sqrt(u-3*v))/(u*v) 4．计算如下表达式： (1) ()() 3542i i -+ (2) () sin 28i - clear,clc (3-5*i)*(4+2*i) sin(2-8*i) 5．判断下面语句的运算结果。 (1) 4 < 20 (2) 4 <= 20

matlab课后习题答案

第4章数值运算 习题 4 及解答 1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分diff 或数值梯度gradient指令计算)(t y'曲线绘制 y'，然后把)(t y和)(t 在同一张图上，观察数值求导的后果。（模拟数据从prob_data401.mat获得） 〖目的〗 ●强调：要非常慎用数值导数计算。 ●练习mat数据文件中数据的获取。 ●实验数据求导的后果 ●把两条曲线绘制在同一图上的一种方法。 〖解答〗 （1）从数据文件获得数据的指令 假如prob_data401.mat文件在当前目录或搜索路径上 clear load prob_data401.mat （2）用diff求导的指令 dt=t(2)-t(1); yc=diff(y)/dt; %注意yc的长度将比y短1 plot(t,y,'b',t(2:end),yc,'r') （3）用gradent求导的指令（图形与上相似） dt=t(2)-t(1); yc=gradient(y)/dt;

plot(t,y,'b',t,yc,'r') grid on 〖说明〗 ● 不到万不得已，不要进行数值求导。 ● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级 以上。 ● 求导会使数据中原有的噪声放大。 2 采用数值计算方法，画出dt t t x y x ? =0 sin )(在]10 ,0[区间曲线，并计算)5.4(y 。 〖提示〗 ● 指定区间内的积分函数可用cumtrapz 指令给出。 ● )5.4(y 在计算要求不太高的地方可用find 指令算得。 〖目的〗 ● 指定区间内的积分函数的数值计算法和cumtrapz 指令。 ● find 指令的应用。 〖解答〗 dt=1e-4; t=0:dt:10; t=t+(t==0)*eps; f=sin(t)./t; s=cumtrapz(f)*dt; plot(t,s,'LineWidth',3) ii=find(t==4.5); s45=s(ii) s45 =

matlab总练习题(完整版)

>> vpa('pi',20) ans = 3.1415926535897932385 >> vpa('exp(1)',20) ans = 2.7182818284590452354 >> x=linspace(-pi,pi,21); >> y=sin(x) y = Columns 1 through 6 -0.0000 -0.3090 -0.5878 -0.8090 -0.9511 -1.0000 Columns 7 through 12 -0.9511 -0.8090 -0.5878 -0.3090 0 0.3090 Columns 13 through 18 0.5878 0.8090 0.9511 1.0000 0.9511 0.8090 Columns 19 through 21 0.5878 0.3090 0.0000 >> k=1:1000; >> kk=1./k; >> kkk=kk./k; >> res=sum(kkk)-(pi^2)/6 res =

-9.9950e-04 %承接上题 >> sum(kk)-log(1000) ans = 0.5777 >> power(1+eps,1/eps) ans = 2.7183 >> a=rand(2,3) a =

0.8147 0.1270 0.6324 0.9058 0.9134 0.0975 >> x=a(1,:) x = 0.8147 0.1270 0.6324 >> y=a(2,:) y = 0.9058 0.9134 0.0975 >> norm(x) ans = 1.0391 >> norm(y) ans = 1.2900 >> acos(dot(x,y)/norm(x)/norm(y)) ans = 0.8189 >> rand(3,3) ans = 0.2785 0.9649 0.9572 0.5469 0.1576 0.4854 0.9575 0.9706 0.8003 >> det(ans) ans =

matlab课后习题答案第一章

第1章基础准备及入门 习题1及解答 1数字1.5e2，1.5e3中的哪个与1500相同吗？ 〖解答〗 1.5e3 2请指出如下5个变量名中，哪些是合法的？ abcd-2xyz_33chan 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 〖解答〗 （1）不同。具体如下 w1=a^(2/3)%仅求出主根 w2=(a^2)^(1/3)%求出(-8)^2的主根 w3=(a^(1/3))^2%求出(-8)主根后再平方 w1= -2.0000+3.4641i w2= 4.0000 w3= -2.0000+3.4641i （2）复数的多方根的，下面是求取全部方根的两种方法： （A）根据复数方根定义 a=-8;n=2;m=3; ma=abs(a);aa=angle(a); for k=1:m%m决定循环次数 sa(k)=(aa+2*pi*(k-1))*n/m;%计算各根的相角 end result=(ma^(2/3)).*exp(j*sa)%计算各根 result= -2.0000+3.4641i 4.0000-0.0000i-2.0000-3.4641i

（B ）利用多项式02 3=-a r 求根p=[1,0,0,-a^2];r=roots(p)r =-2.0000+3.4641i -2.0000-3.4641i 4.00005指令clear,clf,clc 各有什么用处？ 〖解答〗clear 清除工作空间中所有的变量。clf 清除当前图形。clc 清除命令窗口中所有显示。6以下两种说法对吗？（1）“MATLAB 进行数值的表达精度与其 指令窗中的数据显示精度相同。”（2） MATLAB 指令窗中显示的数值有效位数不超过7位。” 〖解答〗 （1）否；（2）否。 7想要在MATLAB 中产生二维数组???? ??????=987654321S ，下面哪些指令能实现目的？ （A ）S=[1,2,3;4,5,6;7,8;9] （B ）S=[123;456;789] （C ）S=[1，2，3；4，5，6；7，8，9] %整个指令在中文状态下输入 〖解答〗 前两种输入方法可以，后一种方法不行。8试为例1.3-5编写一个解题用的M 脚本文件？ 〖解答〗 直接点击新文件图标，出现M 文件编辑器窗口；在该M 文件编辑器中，输入例1.3-5中的全部指令；并另存为p109.m ，便得到所需的脚本文件。

matlab课后习题答案 (附图)

习题2.1 画出下列常见曲线的图形 y （1）立方抛物线3x 命令：syms x y; ezplot('x.^(1/3)') (2)高斯曲线y=e^(-X^2); 命令：clear syms x y; ezplot('exp(-x*x)') (3)笛卡尔曲线

命令：>> clear >> syms x y; >> a=1; >> ezplot(x^3+y^3-3*a*x*y) (4)蔓叶线 命令：>> clear >> syms x y; >> a=1 ezplot(y^2-(x^3)/(a-x)) (5)摆线:()()t sin- = , = - b y 1 t x cos t a 命令：>> clear >> t=0:0.1:2*pi; >> x=t-sin(t);

>>y=2*(1-cos(t)); >> plot(x,y) 7螺旋线 命令：>> clear >> t=0:0.1:2*pi; >> x=cos(t); >> y=sin(t); >> z=t; >>plot3(x,y,z) (8)阿基米德螺线

命令：clear >> theta=0:0.1:2*pi; >> rho1=(theta); >> subplot(1,2,1),polar(theta,rho1) (9) 对数螺线 命令：clear theta=0:0.1:2*pi; rho1=exp(theta); subplot(1,2,1),polar(theta,rho1) (12)心形线

命令：>> clear >> theta=0:0.1:2*pi; >> rho1=1+cos(theta); >> subplot(1,2,1),polar(theta,rho1) 练习2.2 1. 求出下列极限值 （1）n n n n 3 lim 3 +∞ → 命令：>>syms n >>limit((n^3+3^n)^(1/n)) ans = 3 （2）)121(lim n n n n ++-+∞ → 命令：>>syms n >>limit((n+2)^(1/2)-2*(n+1)^(1/2)+n^(1/2),n,inf) ans = 0 （3）x x x 2cot lim 0 → 命令：syms x ；

MATLAB 实用教程 课后习题答案

第二章 1.计算复数3+4i与5-6i的乘积。 a=3+4i b=5-6i c=a*b 2.构建结构体Students，属性包含Name、age和Email，数据包括｛’Zhang’，18，*‘Zhang@16 https://www.doczj.com/doc/5114962598.html,’，’Zhang@https://www.doczj.com/doc/5114962598.html,’+｝、｛’Wang’，21，[]｝和｛’Li’，[],[]｝，构建后读取所有Name属性值，并且修改’Zhang’的Age属性值为19。 Students(1).Age=18 Students(1).Email='zhang@https://www.doczj.com/doc/5114962598.html,','zhang@https://www.doczj.com/doc/5114962598.html,' Students(2).Name='Wang' Students(2).Age=21 Students(2).Email=[] Students(3).Name='Li' Students(3).Age=[] Students(3).Email=[] https://www.doczj.com/doc/5114962598.html, Student(1).Age(1)=19 Student.Age 3.用满矩阵和稀疏矩阵存储方式分别构造下属矩阵： A=[0 1 0 0 0;1 0 0 0 0;0 0 0 0 0;0 0 0 1 0] A=[0 1 0 0 0;1 0 0 0 0;0 0 0 0 0;0 0 0 1 0] S=sparse(A) S=sparse([2,1,4],[1,2,4],[1,1,1],4,5) 4.采用向量构造符得到向量[1,5,9....,41]. A=1:4:41 5.按水平和竖直方向分别合并下述两个矩阵：A=[1 0 0;1 1 0;0 0 1],B=[2 3 4;5 6 7;8 9 10] A=[1 0 0;1 1 0;0 0 1] B=[2 3 4;5 6 7;8 9 10] C=[A B] D=[A;B] 6.分别删除第五题两个结果的第2行。 A=[1 0 0;1 1 0;0 0 1] B=[2 3 4;5 6 7;8 9 10] C=[A B] D=[A;B] C(2,:)=[] D(2,:)=[] 7.分别将第5题两个结果的第2行最后3列的数值改为[11 12 13]。 A=[1 0 0;1 1 0;0 0 1] B=[2 3 4;5 6 7;8 9 10] C=[A B] D=[A;B]