信號與系統(tǒng)-MATLAB-實驗-課后習題答案

第三章練習題

1、

a=[l,1,1];

b=[l,1];

sys=tf(b,a);

t=[0:0.01:10];

figure;

subplot(2,2,1);

step(sys);

subplot(2,2,2);

x_step=zeros(size(t));

x_step(t>0)=l;

x_step(t==0)=l/2;

lsim(sys,x_step,t);

subplot(2,2,3);

impulse(sys,t);

titleCImpulseResponse));

xlabelCTime(sec),)：

ylabel('Amplitude,);

subplot(2,2,4);

x_delta=zeros(size(t));

x_delta(t=0)=100;

[yl,t]=lsim(sys,x_delta,t);

y2=yl;

plot(t,y2);

titleC*ImpulseResponse5);

xlabel('Time(sec)’);

ylabel('Amplitude');

LinearSlmjlatlonResults

hmc(see)Hmo(see)

ImpulseResponse

InpulseResponse

05

Tme(sec)(sec；Time(sec)

2、

函數inti如下：

function[F,tF]=intl(f,tf,a)

T=tf(2)-tf(l);

F=zeros(size(tf));

tF=zeros(size(tf));

tF=tf;

forn=l:length(tf)-l;

F(n+l)=F(n)+T*f(n);

end

驗證如下：

t=[-l:0.01:4];

e=zeros(size(t));

e=(t>-l/2&t<l);

[z,zz]=intl(e,t,-1);

figtire：

plot(zz,z);

1.5

1

0.5

0

-1-0.500511.522.533.54

第四章練習題

1、

Tl=l：

Nl=10000;

tl=linspace(0,T1-T1/N1,Nl)*;

fl=l-2*tl;

OMG=32*pi;

Kl=100;

omg=linspace(-0MG/2,0MG/2-0MG/K1,KI)';

X1=T1/Nl*exp(-j*kron(omg,tl.*))*fl;

fsl=0MG/2/pi/Kl*exp(j*kron(tl,omg.1))*X1;

T2=5;

N2=10000;

t2=linspace(0,T2-T2/N2,N2)';

fs2=O*t2;

f2=sawtooth(t2*2*pi,0);

X2=T2/N2*exp(-j*kron(omg,t2.'))*f2：

fs2=fs2+0MG/2/pi/Kl*exp(j*kron(t2,omg.*))*X2;

figure;

subplot(2,2,1);

plot(omg,abs(XI),'r');

xlabel('Frequency*),ylabel('Amplitude1)

titleC單個鋸齒周期幅頻特性曲線');

subplot(2,2,2);

plot(tl,fsl,'r');

xlabel('Time*),ylabel('Amplitude))

titleCFunctionafterrecovered');

subplot(2,2,3);

plot(omg,abs(X2),'r');

xlabel('Frequency*),ylabel('Amplitude*)

titleC五個鋸齒周期幅頻特性曲線')；

subplot(2,2,4);

plot(t2,fs2,frf);

xlabel('Time*),ylabel('Functionafterrecovered,)

titleCFunctionafterrecovered*);

subplot(3,1,1);

plot(t,f);

titlef驗證原函數')

N=25;

Fl=fsana(t,f,N);

subplot(3,1,2);

stem(abs(Fl),*s*);

titleC前N項傅立葉級數系數幅度曲線')

f2=fssyn(Fl,t);

subplot(3,1,3);

plot(t,f2);

xlabel('time[s]>),ylabel('Amplitude,);

titleC傅立葉逆變換后時域函數');

驗證原函數

1

0

-1

00.10.20.30.40.50.60.70.80.91

前N項傅立葉級數系數幅度曲線

0.4r

□□

0.2?

□□

JJiIilrLlhh??一一：…?

0102030405060

傅立葉逆變換后時域函數

2

①

p

ir

d三

E0

<

-2

00.10.20.30.40.50.60.70.80.91

timefsl

第五章練習題

1、

(a)Residue計算

al=[l/5/6];

bl=[4z5];

[rl/pl,kl]=residue(bLzal);

t=[0:0.01:10];

el=zeros(size(t));

forn=l:size(rl);

el=el+rl(n)*exp(pl(n)*t);

end;

figure;

subplot(1,2,1);

plot(t,el);

title「residue計算,);

xlabel('t/s');

ylabeK'u/v');

lism仿真

sysl=tf([4,5],[1,5,6；);

t=[0:0.01:10];

delta=zeros(size(t));

delta(t==0)=100;

hl=lsim(sysl,delta,t);

subplot(1,2,2)；

plot(t,hl);

title(*lism仿真，)；

xlabel('t/s');

ylabel('u/v');

residue計算lism仿真

Vs

Residue計算和lism仿真結果相同

(b)Residue計算

t=[OzO.01,10];

delta=zeros(size(t));

delta(t==0)=100;

e2=sin(t);

figure;

subplot(2,1,1);

plot(t,e2);

title(*residue計算,);

xlabel('t/s');

ylabel('u/v')；

lism仿真

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

t=[0:0.01:10];

delta=zeros(size(t));

delta(t==0)=100;

h2=lism(sysl,delta,t);

subplot(2,1,2);

plot(h,t2);

axis([0,10,1,1]);

title(*lism仿真');

xlabel('t/s');

ylabel('u/v')；

Residue計算和lism仿真結果相同

2、

S=isstable(sys)函數：

Functions=isstable(sys);

X=ploe(sys);

S=l；

Forn=l:Size(x)

Ifx(n)>0

S=0;break;

End;

End;

穩(wěn)定系統(tǒng)：

Sys=tf(1,[1,2]);

S=isstable(sys);

S=

1

不穩(wěn)定系統(tǒng)：

Sys=tf(1,[1,-2]);

S=isstable(sys);

S=

0

第七章練習題

1、

a=[lzO.5,-0.2,-0.1];

b=[l,-0.3];

n=[0:10]

[hi,t]=impz(b,a,n);

subplot(1,2,1);

stem(n,hi);

u-(n>-0);

hn=filter(b,a,u);

subplot(1,2,2);

stem(n,hn);

2、

nl=[0:9]

n2=[10:19]';

xl=(nl>=0);

x2="(n2>=10);

al=