院党委书记
2009-08-18, 19:49
>> t=linspace(-2*pi,2*pi,512);
w1=100;
w2=3000;
RC=[0.01:0.04:1];
x=cos(100*t)+cos(3000*t);
for k=1:length(RC)
H1(k)=1/(RC(k)*j*w1+1);
H2(k)=1/(RC(k)*j*w2+1);
y(k)=abs(H1(k))*cos(100*t+angle(H1(k)))+abs(H2(k))*cos(w2*t+angle(H2(k)));
figure(k);
subplot(2,1,1);
plot(t,y(k),'r');
xlabel('\omega');
ylabel('|H(j\omega)|');
title(['\bf当RC取=',num2str(RC(k))],'fontsize',15);
subplot(2,1,2);
plot(t,x ,'g');
xlabel('\omega');
ylabel('x(t)');
end
??? In an assignment A(I) = B, the number of elements in B and
I must be the same.
x为原波形 y为经滤波器后的波 题目让找合适的RC值 使x中的高频虑掉 我找了RC从0.01开始到1之间的数
编出后程序没有提示error 我估计是逻辑错误 可找了半天找不出来
w1=100;
w2=3000;
RC=[0.01:0.04:1];
x=cos(100*t)+cos(3000*t);
for k=1:length(RC)
H1(k)=1/(RC(k)*j*w1+1);
H2(k)=1/(RC(k)*j*w2+1);
y(k)=abs(H1(k))*cos(100*t+angle(H1(k)))+abs(H2(k))*cos(w2*t+angle(H2(k)));
figure(k);
subplot(2,1,1);
plot(t,y(k),'r');
xlabel('\omega');
ylabel('|H(j\omega)|');
title(['\bf当RC取=',num2str(RC(k))],'fontsize',15);
subplot(2,1,2);
plot(t,x ,'g');
xlabel('\omega');
ylabel('x(t)');
end
??? In an assignment A(I) = B, the number of elements in B and
I must be the same.
x为原波形 y为经滤波器后的波 题目让找合适的RC值 使x中的高频虑掉 我找了RC从0.01开始到1之间的数
编出后程序没有提示error 我估计是逻辑错误 可找了半天找不出来