| 查看: 891 | 回复: 3 | ||
[求助]
求助matlab编程
|
|
下面matlab编程中，经过fft2变换后，k和fftk应该怎么表示呢？下面中的k和fftk是有问题的，但是我不清楚怎么去正确表示，请指教。 n1=100;n2=100; for x=1:2*n1 for y=1:2*n2 test(x,y)=exp(sin(2*(x-1)*pi/n1+3*(y-1)*pi/n2)); testdx(x,y)=2*cos(2*(x-1)*pi/n1).*exp(sin(2*(x-1)*pi/n1+3*(y-1)*pi/n2)); end end testk=fft2(test,2*n1,2*n2); k=[0:n1-1,-n1:-1];%问题在此处以及下一行 fftk=repmat(k,200,1）； testdxk = fftk.*testk*i; testdx2 = ifft2(testdxk,2*n1,2*n1); figure; mesh(real(testdx2)); title('real part'); figure; mesh(testdx); title('anaylysis part'); |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有49人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
三无产品还有机会吗
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助MATLAB做三维图问题
已经有4人回复
- 求助一个matlab关于分类讨论的问题
已经有9人回复
- matlab画图求助
已经有11人回复
- 求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟
已经有10人回复
- 求助matlab编程
已经有13人回复
- 【求助】求助 matlab编程：可否实现清除文本文档数据
已经有5人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】matlab中如何将程序的结果输出到一个文件中
已经有11人回复
- 【求助】多目标pso算法的MATLAB程序
已经有3人回复
- 【求助】请教把java程序转化成matlab程序
已经有9人回复
- 【求助】初学matlab
已经有11人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
liusir3836
木虫 (职业作家)
猎手
- 博学EPI: 64
- 应助: 54 (初中生)
- 贵宾: 0.783
- 金币: 5493.3
- 散金: 17900
- 红花: 60
- 沙发: 2
- 帖子: 4698
- 在线: 4158.3小时
- 虫号: 708030
- 注册: 2009-02-25
- 专业: 数论
2楼2013-04-08 10:48:49
Toauto
至尊木虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 13024.8
- 散金: 1904
- 红花: 25
- 帖子: 1925
- 在线: 417.2小时
- 虫号: 339744
- 注册: 2007-04-07
- 专业: 文献学
3楼2013-04-08 13:08:55
|
n1=100;n2=100; for x=1:2*n1 for y=1:2*n2 test(x,y)=exp(sin(2*(x-1)*pi/n1+3*(y-1)*pi/n2)); testdx(x,y)=2*cos(2*(x-1)*pi/n1).*exp(sin(2*(x-1)*pi/n1+3*(y-1)*pi/n2)); end end testk=fft2(test,2*n1,2*n2); k=[0:n1-1,-n1:-1]; fftk=repmat(k,200,1); testdxk = fftk.*testk*i; testdx2 = ifft2(testdxk,2*n1,2*n1); figure; mesh(real(testdx2)); title('real part'); figure; mesh(testdx); title('anaylysis part'); error =abs(testdx-testdx2); figure; mesh(error); 右边的括号是我打错了。这个程序可以运行出来，但问题在于结果不对，就是误差error很大，正确的误差在10的负13次方左右，而上面运行出来的结果在0到10之间，我经过分段测试后，发现问题在k和fftk上，我不知道怎么去定义这个200行200列的矩阵fftk，请指教。 |
4楼2013-04-08 21:53:34