t=0:0.1:30;
y=sinc(t);
subplot(2,1,1);
plot(t,y);
xlabel('t');ylabel('y');
title('时限信号');
axis([0 30 -1 1]);
N=1000;
k=0:N;
w1=3.5;
w=k*w1/N;
F=y*exp(-i*t'*w);
subplot(2,1,2);
plot(w,F);
xlabel('w');
ylabel('F');
title('时限信号的频谱图');
ws=7;
Ts=2*pi6;
n=-10:Ts:10;
x=sinc(n*Ts);
k=0:length(n)-1;
subplot(3,1,1);
stem(k,x);grid;
xlabel('k');ylabel('x');title('w=7临界采样');
X=x*exp(-i*n'*w);
subplot(3,1,2);
plot(w,X);grid;
xlabel('W');ylabel('X');title('w=7临界采样频谱');
axis([0 10 0 5]);
[B,A]=butter(10,3/pi); %n为滤波器阶数,Wn是自然频率
[H,m]=freqz(B,A,512,2000); %计算滤器频率响应函数,如果n没指定就是512
fa=filter(B,A,x); %滤除向量x中的数据,b是分子系数,a是分母系数
subplot(3,1,3);
plot(fa);
xlabel('t');ylabel('y');
title('重构时间信号');
axis([0 30 -1 1]);
ws=5;
Ts=2*pi;
n=-10:Ts:10;
x=sinc(n*Ts);
k=0:length(n)-1;
figure;subplot(2,2,1);
stem(k,x);grid;
xlabel('k');ylabel('x');title('w=5欠采样');
X=x*exp(-i*n'*w);
subplot(2,2,2);
plot(w,X);grid;
xlabel('W');ylabel('X');title('w=5欠采样频谱');
axis([0 10 0 5]);
[B,A]=butter(10,3/pi);
[H,m]=freqz(B,A,512,2000);
fa=filter(B,A,x);
subplot(2,2,3);
plot(fa);
xlabel('t');ylabel('y');
title('重构时间信号');
axis([0 30 -1 1]);
ws=10;Ts=2*pi;n=-10:Ts:10;x=sinc(n*Ts);k=0:length(n)-1;figure;subplot(3,2,1);stem(k,x);grid;xlabel('k');ylabel('x');title('w=10过采样');X=x*exp(-i*n'*w);subplot(2,2,2);plot(w,X);grid;xlabel('W');ylabel('X');title('w=10过采样频谱');axis([0 10 0 5]);[B,A]=butter(10,3/pi);[H,m]=freqz(B,A,512,2000);fa=filter(B,A,x);subplot(2,2,3);plot(fa);xlabel('t');ylabel('y');title('重构时间信号');axis([0 30 -1 1]);
【matlab求教】有学长能来解释一下各行代码是什么意思嘛
最新推荐文章于 2024-09-26 21:03:52 发布