信号与系统课程设计实验三
⑴ 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,'.');
最后那个看不懂,因为没有图;其思路应该是
先求方波的各个谐波分量,再求每个分量产生的响应,再加起来,自己参考上面的