習(xí)題答案-MATLAB和LabVIEW仿真技術(shù)及應(yīng)用實例(第2版)-聶春艷

1章MATLAB語言突出的特點(diǎn)是什么？MATLAB放性強(qiáng)等特點(diǎn)。MATLAB為獨(dú)立窗口？又如何將脫離出去的窗口重新放置到桌面上？MATLAB5CloseUndockUndock按鈕就可以使該窗口脫離桌desktopDock……菜單項就可以將獨(dú)立的窗口重新放置的桌面上。2章4327i 175i 2.1求x379344i >>x=[4+8i3+5i2-7i1+4i7-5i;3+2i7-6i9+4i3-9i4+4i];>>y=x’y=4.0000-8.0000i 3.0000-2.0000i3.0000-5.0000i 7.0000+6.0000i2.0000+7.0000i 9.0000-4.0000i1.0000-4.0000i 3.0000+9.0000i7.0000+5.0000i 4.0000-4.0000i2.2計算a

6 9 3與b2 4 1的數(shù)組乘積。2 7 5

4 6 8 >>a=[693;275];>>b=[241;468];>>a.*bans=12 36 38 42 40AXB

4 9 2A7 6 4 B

37262-3對于

，如果

， X。3 5 7

28A=[492;764;357]; AX=B A-1AX=A-1B X=A-1*B=A\BA/B=AB*-1 、B=[37;26;28];X=A\BX=-0.51184.04271.33181 2 3 2-4已知：a4 5 6，分別計算a 7 8 9觀察其結(jié)果。>>a=[123;456;789];1 2 31 2 3 >>a.^2 4 5 6 7 8 97 8 9ans=1 4 916 25 3649 64 81>>a^2ans=3036426681961021261502.5 a

2 5 b

7 46 2ab之間的六種關(guān)系運(yùn) 算的結(jié)果。>>a=[125;36-4];>>b=[874;362];>>a>bans=011000>>a>=b>>a<b....>>a<=b....>>a==b....>>a~=b....3章3.1繪制曲線yx3x1，x的取值范圍為[-5,5]。解：>>x=-5:0.2:5;>>y=x.^3+x+1;>>plot(x,y)題圖yeat，t[0，10]的線型和標(biāo)記點(diǎn)畫出a=0.1、a=0.2、a=0.5三種情況下的曲線。t=0:0.5:10;a1=0.1;a2=0.2;a3=0.5;y1=exp(-a1*t);y2=exp(-a2*t);y3=exp(-a3*t);plot(t,y1,'-ob',t,y2,':*r',t,y3,'-.^g')subpolt怎樣編程？t=0:0.5:10;subplota1=0.1;y1=exp(-a1*t);plot(t,y1,'-ob')subplot(312)a2=0.2;y2=exp(-a2*t);plot(t,y2,':*r')subplot(313)a3=0.5;y3=exp(-a3*t);plot(t,y3,'-.^g')題圖3-2yeata加圖例。直接在Insert中的title中加y=e^{-\itat}---a,t都是斜體e^{-a\itt}----，t是斜體，不能用小括號3.3題圖添加圖例直接在Insert中點(diǎn)擊legend3-4 有[5，6，8，0，2，3，7，4，1，2]10個數(shù)據(jù)，編制程序并畫出數(shù)據(jù)點(diǎn)變化的曲線，要求標(biāo)出圖中X軸，Y軸標(biāo)識。解：>>x=[5,6,8,0,2,3,7,4,1,2,3];>>plot(x)>>plot(x,'o')>>plot(x,'-ob')8765432101 2 3 4 5 6 7 8 9 10 113-4題圖35正弦信號為sin(100 t),噪聲為0.5randn(size(t)分別繪制信號噪聲和疊加后的3條曲線，要求在同一坐標(biāo)里。解：t=0:0.005:1;a=sin(100*pi*t);n=0.5*randn(size(t));s=a+n;plot(t,a,'r',t,n,'b',t,s,'g')3210-1-2-30 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 13-5題圖4章用符號方法求解非線性方程組x22xyy23的數(shù)值解。x24x30解：>>syms xy>>[x,y]=solve('x^2-2*x*y+y^2=3','x^2-4*x+3=0');>>solution=[x,y]solution=[1,3^(1/2)+1][3,3^(1/2)+3][3,3-3^(1/2)][1,1-3^(1/2)]用ezplot符號簡易繪圖法繪制表達(dá)式5ex(sinxcosx)的二維圖形。解：>>syms xy>>y=5*exp(-x)*(sin(x)-cos(x));>>ezplot(y)-(5(cos(x)-sin(x)))/exp(x)4002000-200-6 -4 -2 0x

2 4 64-2題圖用ezmeshzx3y3其中x和y[-，。解：>>symsxy>>z=x^3+y^3;>>ezmesh(z,[-5,5])x3+y32000-200550 0y -5 -5 x4-3題圖5章命令文件與函數(shù)文件的主要區(qū)別是什么？命令文件是一系列命令的組合，函數(shù)文件的第一行必須用function說明。函數(shù)文件可以接受輸入?yún)?shù)，也可以返回參數(shù)，相當(dāng)于子程序。y=x2+cos(x)量，調(diào)用該函數(shù)文件，求運(yùn)行結(jié)果。解：在編輯器里：functiony=hansu(x)y=x.^2+cos(x);在命令窗口調(diào)用：a=3;b=hansu(a)b=8.01005-3編制函數(shù)文件實現(xiàn)直角坐標(biāo)（x,y）與極坐標(biāo)（,）之間的轉(zhuǎn)換。

，極坐標(biāo)的極角： y。x2yx2y2arctan(在編輯器里：function[p,o]=zj(x,y)p=sqrt(x.^2+y.^2);o=atan(y./x);在命令窗口調(diào)用：調(diào)用：a=6;d=9;[m,n]=zj(a,d)m=10.8167n=0.982867章7-1155倍的結(jié)果合成一個二維數(shù)組，用數(shù)字顯示器顯示出來。題圖x已知數(shù)學(xué)模型示x的波形。x

0.5sint，建立仿真模型，并在示波器中顯=0.5

+sint7-2題圖 已知數(shù)學(xué)模型x2x5x30.9sin(50t，利用Simulink 立仿真模型，并顯示x變化曲線以及x與

之間關(guān)系曲線，同時將x和數(shù)據(jù)分別送入工作空間。+2

-5x3

=0.9sin50t

=0.9sin50t-2

+5x3題圖圖中是xdot與x）已知控制系統(tǒng)的傳遞函數(shù)為 1 用Simulink建立s28系統(tǒng)仿真模型，并對系統(tǒng)的階躍響應(yīng)進(jìn)行仿真。77-4題圖2V、頻率0.5Hz0.1V號顯示在示波器中并將結(jié)果傳送到工作空間。7-5題圖第8章略。略。象作為將創(chuàng)建的子系統(tǒng)部分，要求創(chuàng)建子系統(tǒng)并進(jìn)行子系統(tǒng)封裝。8-3題圖第9章如圖所示的電路，已知Us=4V，Is=4A，R1=3Ω，R2=5Ω，求電路中的電壓。R2R2+IsUR1Us-9-1題圖9-36所示電路，已知壓U。R2R2+IsUR1U-ContinuousContinuouspowerguiSeriesRLCBranchpowerguiSeriesRLCBranchACCurrentSource+ControlledVoltageSourceSeriesRLCBranch1s --VoltageMeasurement+v2題圖二階低通電路如圖所示，若以U2

為輸出時，求頻率響應(yīng)函數(shù)，并畫出幅頻特性和相頻特性。+L+L+U1CRU2-解：1)建立模型11U R 1H(j) 2 1 1(jQ1LC式中 LCn

1 n Rn CLQCL幅頻響應(yīng)為：G20logH(j)相頻特性為：()arg[H(j)]橫坐標(biāo)取

0.1, , 1, , 10Q11

1 畫圖當(dāng) 1時，22Q22n 32稱為最平幅度特性，即在通帶內(nèi)其幅頻特性最平坦。2)根據(jù)模型編寫M文件進(jìn)行仿真根據(jù)對電路的分析編寫M文件如下：clear,formatcompactn=50;%設(shè)置數(shù)據(jù)數(shù)目ww=logspace(-1,1,n);組ww=w/wnforQ=[1/3,1/2,1/sqrt(2),1,2,5]H=1./(1-(ww).^2+j*ww/Q);%設(shè)置頻率數(shù)%求復(fù)頻率響應(yīng)figure(1)subplot(2,1,1),plot(ww,abs(H)),holdon%繪制幅頻特性subplot(2,1,2),plot(ww,angle(H)),holdon%繪制相頻特性figure(2)subplot(2,1,1),semilogx(ww,20*log10(abs(H))),holdon%繪制對數(shù)頻率特性subplot(2,1,2),semilogx(ww,angle(H)),holdon%繪制相頻特性endfigure(1),subplot(2,1,1),grid,xlabel('ww'),ylabel('abs(H)') %對坐標(biāo)軸進(jìn)行標(biāo)注subplot(2,1,2),grid,xlabel('ww'),ylabel('angle(H)')figure(2),subplot(2,1,1),grid,xlabel('ww'),ylabel('分貝')subplot(2,1,2),grid,xlabel('ww'),ylabel('angle(H)')線性頻率特性 (b)對數(shù)頻率特性頻率特性9-3題圖A0A1A2Y7Y6Y5Y4Y3Y2Y1Y0000A0A1A2Y7Y6Y5Y4Y3Y2Y1Y000010011010101111001101111011111利用Simulink81提示：真值表如下：8選1數(shù)據(jù)選擇器真值表輸入輸出DEA2A1A0Y1××X0D000001D100011D200101D300111D401001D501011D601101D701111邏輯表達(dá)式如下：Y=

210D0+

21A0D1。。。。利用Simulink模塊創(chuàng)建一個半加器的仿真模型。1.真值表輸入輸出ABY0000111011102.邏輯表達(dá)式Y(jié)=A(-)B+AB(-)9-610章10-1/輸出微分方程為dt激響應(yīng)和階躍響應(yīng)。解：a=[2];b=[1,3];figure(1),impulse(a,b);figure(2),step(a,b)ImpulseResponse2ppmA

1.510.500 0.5 1Time(sec)

1.5 2StepResponseppmA

0.60.40.200 0.5 1 1.5 210-1題圖/輸出微分方程為d2

若x(t)6(t)，求其零狀態(tài)響應(yīng)。解：

dt2 dt>>a=[1];>>b=[143];>>t=1:0.1:4;>>sin(t).*ones(1,length(t))

x=6*ones(1,length(t));!!!!! 11111>>lsim(a,b,x,t)ppmA

LinearSimulationResults65432101 1.5 2 2.5Time(sec)10-2題圖

3 3.5 4略

H(s) s3s44s316s212s轉(zhuǎn)換成狀態(tài)空間形式和零極點(diǎn)增益模型形式。解：>>a=[13];>>b=[1416120];>>[A,B,C,D]=tf2ss(a,b)A=-4 -16 -121 0 00 1 0B=100C=0 1 3

AxBuyCxDuD=0y5y7y'3yu''3u'2u求出它的傳遞函數(shù)模型、零極點(diǎn)增益模型和狀態(tài)空間模型。解：兩 邊 拉 氏 變 換S35S2Y(s)S22U(s)H(s)

Y(s)U(s)

所以，

H(s)S3

S23S25S27S3零極點(diǎn)增益求法：a=[132];b=[1573];[Z,P,K]=tf2zp(a,b)狀態(tài)空間模型求法a=[132];b=[1573];[A,B,C,D]=tf2ss(a,b)11章11-1已知x(n)2n5 n1，y(n)n

n0

()與()之和。

0 n

3(2n)

n0 xn yn解：clear,formatcompactn=-10:10;x=(0).*(n<-1)+(5+2.^-n).*(n>=-1);y=((2.^n).*3).*(n<0)+(n+2).*(n>=0);sum=x+y;subplot(3,1,1),stem(n,x,'.');axis([-5,10,-1,8]);ylabel('X');subplot(3,1,2),stem(n,y,'.');axis([-5,10,--1,13]);ylabel('Y');subplot(3,1,3),stem(n,sum,'.');axis([-5,10,-1,21]);ylabel('sum=X+Y');86X 420-512108Y 642-5Y 20+=X 10=0mus0

0 5 100 5 10-5 0 5 10題11-1圖已知序列x(n)2n5 n1，求x(n)的反褶序列。0 n10解：clear,formatcompactn=-10:10;x=(0).*(n<1)+(2.^-n+5).*(n>=1);y=fliplr(x);n1=-fliplr(n);subplot(2,1,1);stem(n,x);xlabel('n'),ylabel('x(n)');subplot(2,1,2);stem(n1,y);xlabel('n'),ylabel('y(n)=x(-n)');64x 20-106) 4n 2y0

-5 0n

5 10-10 -5 0n

5 1011-2圖用Hanning窗設(shè)計一個滿足下列指標(biāo)的 FIR低通濾波器：s p， ，Rp3dB，R 40s 解：clear,formatcompactwp=0.2*pi;ws=0.5*pi;ww=ws-wp;M=ceil(8*pi/ww);wc=(wp+ws)/2;window=hanning(M+1);b=fir1(M,wc/pi,window);freqz(b,1,512)M

500-50-100-1500 0.2 0.4 0.6 0.8 1NormalizedFrequency(rad/sample)e(eaP

0-500-1000-1500

0 0.2 0.4 0.6 0.8 1NormalizedFrequency(rad/sample)11-3圖24階FIR0.2~0.5解：clear,formatcompactw=[0.2,0.5];b=fir1(24,w);freqz(b);figure;stem(b,'.');line([0,25],[0,0]);xlabel('n');ylabel('h(n)');M

1000-100-2001000

0 0.2 0.4 0.6 0.8 1NormalizedFrequency(rad/sample)e(ea

0-1000P -2000

0 0.2 0.4 0.6 0.8 1NormalizedFrequency(rad/sample)0.60.4h 0.20-0.20 5 10 15 20 25n11-4圖12章12-1在［-5,5］區(qū)間上分別用10個數(shù)據(jù)點(diǎn)和60個數(shù)據(jù)點(diǎn)，采用樣條插值函數(shù)繪制函數(shù)的圖形。解：x=linspace(-5,5,10);y=5+x.^3;xi=linspace(-5,5,60);yi=interp1(x,y,xi,'spline');plot(x,y,'or',xi,yi)150100500-50-100-150-5 0 5題圖12-2對數(shù)據(jù)x=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1]和y=[-0.232 0.647 1.877 3.565 5.134 7.443 9.221 10.01111.678 12.566 13.788]線。解：x=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1];y=[-0.232 0.647 1.877 3.565 5.134 7.443 9.22110.01111.678 12.566 13.788];p=polyfit(x,y,2);x1=linspace(0,1,50);y1=polyval(p,x1);plot(x,y,'o',x1,y1)151050-50 0.2 0.4 0.6 0.8 1題圖410二維條形圖方法繪制數(shù)據(jù)分析圖。第1天第2天第3天第4天第5天第6天第7天第8天第9天10天第1臺20232722283032252024第2臺22242524202830242125第3臺23252431293031322940第4臺27263035323340363329解：Y=[20232722283032252024;22242524202830242125;23252431293031322940;27263035323340363329];bar(Y,'group')4030201001 2 3 412-3題圖13章13-1G(s)

3s42s35s24

，繪制其零極點(diǎn)圖，并判斷系統(tǒng)的穩(wěn)定性。解：num=[32504];den=[134586];pzmap(num,den)

s53s44s35s28s6Pole-ZeroPole-ZeroM ap10.5doe(AI

0-0.5-1-1.5--RealAxis(seconds )13-1此系統(tǒng)有位于s右半平面，所以系統(tǒng)不穩(wěn)定。13-2

G(s)

s22s4s310s25s

，繪制單位階躍響應(yīng)曲線圖。解：G1=tf([124],[1