信號與系統(tǒng)-MATLAB-實驗-課后習題答案_第1頁
信號與系統(tǒng)-MATLAB-實驗-課后習題答案_第2頁
信號與系統(tǒng)-MATLAB-實驗-課后習題答案_第3頁
信號與系統(tǒng)-MATLAB-實驗-課后習題答案_第4頁
信號與系統(tǒng)-MATLAB-實驗-課后習題答案_第5頁
已閱讀5頁,還剩4頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領

文檔簡介

第三章練習題

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=

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容,請與我們聯系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論