|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我在simulink里面建模,并且取出来了一组数,想对这组数进行离散傅里叶变换,下面是我的程序,运行的时候出现问题,我是按照一本书的一个例题改的,下面是程序,请各位帮我看看是什么问题,多谢了啊
立体程序是正确的,x(n)={1,2,3,4,5}
n=-1:3;x=1:5;
k=0:500;w=(pi/500)*k;
X=x*exp((-j*pi/500)).^(n'*k);
magX=abs(X);angX=angle(X);
realX=real(X);imagX=imag(X);
subplot(2,2,1);plot(k/500,magX);grid
xlabel('frequency in pi units');title('Magnitude Part')
subplot(2,2,3);plot(k/500,angX);grid
xlabel('frequency in pi units');title('Angle Part')
subplot(2,2,2);plot(k/500,realX);grid
xlabel('frequency in pi units');title('real part')
subplot(2,2,4);plot(k/500,imagX);grid
xlabel('frequency in pi units');title('Imaginary Part')
我编的程序,我从simulink里面读出数据,使用了to workspace的模块,读出的数据存为out1,out1=[0.4012 0.3630 0.3285 0.2972 0.2689],要对这组数进行傅里叶变换,我的程序是n=-1:3;x=0.4012 0.3630 0.3285 0.2972 0.2689;
k=0:500;w=(pi/500)*k;
X=x*exp((-j*pi/500)).^(n'*k);
magX=abs(X);angX=angle(X);
realX=real(X);imagX=imag(X);
subplot(2,2,1);plot(k/500,magX);grid
xlabel('frequency in pi units');title('Magnitude Part')
subplot(2,2,3);plot(k/500,angX);grid
xlabel('frequency in pi units');title('Angle Part')
subplot(2,2,2);plot(k/500,realX);grid
xlabel('frequency in pi units');title('real part')
subplot(2,2,4);plot(k/500,imagX);grid
xlabel('frequency in pi units');title('Imaginary Part')但是出现了错误,所以修改第一行为n=-1:3;x=0.4012 ,0.3630 ,0.3285 ,0.2972 ,0.2689;
但是执行完程序以后出现了四幅图,就是说他对这四个数分别进行了傅里叶变换,我不知道问题在哪里,请各位指点迷津。。。多谢,O(∩_∩)O~ |
|