| 查看: 1371 | 回复: 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'); |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助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
- 帖子: 4699
- 在线: 4158.7小时
- 虫号: 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
