哈工大信號與系統(tǒng)

KKb0問題2已知傳遞函數(shù)同 19問題3已知系統(tǒng)的零極點(diǎn)z=-1;p=[-2–3–4];[num,den]=zp2tf(-1,[-2–4],5)（回車原理：兩個信號卷 對于兩個不規(guī)則波形的卷積，依靠手算是很 的， 種則變得十分簡單。 n1單位階躍序列u(n)1ones函數(shù)實(shí)現(xiàn)

n

這一函數(shù)可以用nn1:x=ones(1,N); 也可以表示成

x[(nn0)0]f1(t)u(t1)u(tf2(t)u(t2)u(tu(t1)u(t u(t2)u(t求其卷積C(t)f1(t)*f2(t)的 t1=1:0.01:2;f1=ones(size(t1)).*(t1>1);（表示一個高度為1的門函數(shù)，時(shí)間從t=1t=2）t2=2:0.01:3;f2=ones(size(t2)).*(t2>2); （表示一個高度為1的門函數(shù)，時(shí)間從t=2t=3）c=conv(f1,f2);（卷積t3=3:0.01:5;其結(jié)果如圖所示c(n1)

N1a(k1)b(nkkNakcconv(a,c 1f1(t)t tf2(t)

t 計(jì)算卷積c(t) f1(t)*f2 程序如下：t2=-1:0.01:2;f2=t2.*exp(-t2).*(t2>0)+exp(t2).*(t2<0);t3=-） 3（1）門脈沖信號：symstut=sym(‘Heaviside(t+0.5)-Heaviside(t-0.5)’);Fw=fourier(ut,t,w);FFw=maple(‘convert’,Fw,’piecewise’);FFFw=FFw;FFP=abs(FFFwEzplot(FFP,[-10*pi10*pi])在10,10]fft函數(shù)的波形Axis([-10*pi10*pi01]);x軸在10,10]y軸在[01]//繪制函數(shù)的外觀和尺度FFFw=FFw*FFw;Axis([-10*pi10*pi01]);問題：問題：1）f(t)et(t信號f(t)et2幅度譜。解答：symstwft=exp(-t)*ut;fwp=abs(fw)ezplot(fwp)symsw實(shí)驗(yàn)二系統(tǒng)的時(shí)間響應(yīng)分析二階系統(tǒng)不同阻尼比ξ時(shí)的階躍響應(yīng)二階系統(tǒng)的傳遞函數(shù)為H(s) ，其階躍響應(yīng)可用下n2 2nsn2列程序做出其單位階躍和沖激響應(yīng)的曲線（為簡單起見，令n 1如圖（ab）所示。CloseHoldonzeta=[0.10.20.40.71.0];num=[1]t=0:.01:12;fork=1:5den=[12*zeta(k)1];printsys(num,den,’s’);[y(:,k),x]=step(num,den,1);plot(t,y(:,k));endfork=1:5den=[1Zeta(k)[y(:,k),x]=impulse(num,den,t);plot(t,y(:,k));holdonats；b.tr；c.σ；closeall; t=0.03:0.03:3y=step(num,den,t);fork=1:99;tt=k+1if(((y(k)-0.1)*(y(tt)-0.1))<=0trb=t(k)+0.03*((0.1-y(k))/(y(tt)-y(k)));else (((y(k)-0.9)*(y(tt)-0.9))<=0tre=t(k)+0.03*((0.9-y(k))/(y(tt)-y(k)));break;endendtr(n)=tre-trbendholdonplot(zeta,ts);xlabel(‘Ts’);holdoff（b）5

s(5s

5(1s(5s3

s(5s

num=5;den=[510];55s^2

subplot(2,1,1);step(num1,den1);holdonnum=[45];den=[510];printsys(num,den,’s’)(運(yùn)行結(jié)果：sign=-1

num/den=

4s5s^2[num2,den2]=cloop(num,den,sign);printsys(num2,den2,’s’)(運(yùn)行結(jié)果：num/den=

4s5s^25snum=[05];den=[510];(運(yùn)行結(jié)果：55s^2 1];denf=[01];sign=-1;(

5s^25sstep(num3,den3);axis([01001.8]);subplot(2,1,2);holdonimpulse(num1,den1);impulse(num2,den2);impulse(num3,den3);axis([010–12]);holdoff4．不同阻尼比下二階系統(tǒng)的波特圖closezeta=[0.10.20.40.71.0];fork=1:5num=[9];den=[12*zeta(k)9];printsys(num,den,'s');w=1.5:.1:10;iw=i*w;H=p