信號與系統課程設計實驗三

⑴ matlab中信號與系統的實驗三（數字信號部分1）

f=zeros(1,12)
f(2)=1,f(12)=-1
k=0:100;
h=0.9.^k;
y=conv(f,h);
figure
subplot(311)
stem(0:11,f,'filled');
axis([0 100 -2 2])
title('f')
subplot(312)
stem(k,h,'filled');
title('h')
subplot(313)
stem(k,y(1:length(k)),'filled');
title('y')

⑵ 自動控制/計算機控制/信號與系統三合一教學實驗系統

這么復雜的問題.....
要是懸賞個兩百分說不定才會有人下功夫幫你找了....

⑶ 信號與系統實驗

%
f0=50;N=[3,9,15,50];T=1/f0;t=[0:T/100:2*T];
for i=1:length(N);y=0;
for m=1:N(i);
y=y+1/(2*m-1)*cos((2*m-1)*f0*2*pi*t);
end
figure(i);
plot(t,y);
end
%三角波
K=[3,9,15,50]; T=2; tor=T;w0=2*pi/T;syms t k;
for L=1:length(K);N=K([L]);
xet= (1-abs(t))*exp(-j*k* w0*t);
a=int(xet,t,-tor/2,tor/2)/T;
n=[-N:N]; n=n+eps; ak=subs(a,n); ak([N+1])=tor/2/T;
expkk=exp(n*j* w0*t); xN=simple(ak*expkk');
figure(L); ezplot(xN,[-T/2,T/2]);end
%方波
K=[3,9,15,50]; tor=1;T=2; w0=2*pi/T;syms t k;
for L=1:length(K);N=K([L]);
xet=exp(-j*k* w0*t);
a=int(xet,t,-tor/2,tor/2)/T;
n=[-N:N]; n([N+1])=eps; ak=subs(a,n);
expkk=exp(n*j* w0*t); xN=simple(ak*expkk');
hold on; ezplot(xN,[-T/2,T/2]);end
%hold on改為figure(L)即可畫在不同的圖上；
%方波的頻譜圖:
N=80;T=20;tor=1;syms t k;
xet=exp(-j*k*2*pi/T*t);
a=int(xet,t,-tor/2,tor/2)/T;
n=[-N:N]; n=n+eps;ak=subs(a,n);
stem(n,ak,'.');

最後那個看不懂,因為沒有圖;其思路應該是
先求方波的各個諧波分量,再求每個分量產生的響應,再加起來,自己參考上面的