| 查看: 854 | 回复: 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'); |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有6人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有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
- 金币: 5499.3
- 散金: 17900
- 红花: 60
- 沙发: 2
- 帖子: 4697
- 在线: 4158.2小时
- 虫号: 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