UC知道时域采样理论的验证程序清单,运行后无图像
来源：百度知道 编辑：UC知道 时间：2024/06/05 21:06:29
《数字信号处理》，西安电子科技大学出版社，实验二时域采样理论的验证程序，修改后得出主要程序
Tp=64/1000;
Fs=1000;T=1/Fs;
N=Tp*Fs;
n=0:N-1;
A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;
xnt=A*exp(-alph*n*T).*sin(omega*n*T);
xk=T*fft(xnt,N);
wk=2*pi*k/N;
fk=k/Tp;
yn='xa(n*T)';
subplot(2,1,1);
stem(xnt,yn);
k=0:N-1;fk=k/Tp;
subplot(2,1,2);
plot(fk,abs(xk));
运行后除了边框，什么图像也没有，也没有报错提示。哪里出错了？
Tp=64/1000;
Fs=1000;T=1/Fs;
N=Tp*Fs;
n=0:N-1;
A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;
xnt=A*exp(-alph*n*T).*sin(omega*n*T);
xk=T*fft(xnt,N);
wk=2*pi*k/N;
fk=k/Tp;
yn='xa(n*T)';
subplot(2,1,1);
stem(xnt,yn);
k=0:N-1;fk=k/Tp;
subplot(2,1,2);
plot(fk,abs(xk));
运行后除了边框，什么图像也没有，也没有报错提示。哪里出错了？
yn='xa(n*T)'; 这句不对，用下面的：
Tp=64/1000;
Fs=1000;T=1/Fs;
N=Tp*Fs;
n=0:N-1;
A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;
xnt=A*exp(-alph*n*T).*sin(omega*n*T);
xk=T*fft(xnt,N);
wk=2*pi*n/N;
fk=n/Tp;
subplot(2,1,1);
stem(fk,xnt);
k=0:N-1;fk=k/Tp;
subplot(2,1,2);
plot(fk,abs(xk));