24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1049  |  回复: 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');

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

陈昌荣

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by Toauto at 2013-04-08 13:08:55
请楼主贴一下出错信息,另外fftk=repmat(k,200,1);的右括号成了全角的了

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

liusir3836

木虫 (职业作家)

猎手

优秀版主

帮你联系下Toauto,看能不能帮到你。
Toauto你好!能不能帮忙看下楼主的问题,谢谢啦!@Toauto
态度确定行为
2楼2013-04-08 10:48:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Toauto

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★
liusir3836: 金币+1, 感谢过来应助!谢谢! 2013-04-08 13:33:55
wolfghost: 金币+2, 欢迎参加有奖问答版块活动~金币奖励 2013-04-10 16:18:04
请楼主贴一下出错信息,另外fftk=repmat(k,200,1);的右括号成了全角的了
low-key
3楼2013-04-08 13:08:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 272分材料子求调剂 +40 Loy0361 2026-04-10 52/2600 2026-04-12 19:23 by laoshidan
[考研] 求调剂 +16 张番茄不炒蛋 2026-04-10 17/850 2026-04-12 13:58 by 熬夜成!
[考研] 297求调剂 +17 ORCHID1 2026-04-10 18/900 2026-04-12 09:59 by zhouwenxian
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 材料工程日语考生求调剂 +7 0856?调剂 2026-04-10 7/350 2026-04-11 21:33 by 蓝云思雨
[考研] 求调剂 +6 电气300求调剂不 2026-04-08 6/300 2026-04-11 20:14 by 逆水乘风
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 346,工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 342电子信息专硕求调剂 +9 你让我怎么荔枝 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 本科西工大 0856 324求调剂 +10 wysyjs25 2026-04-09 11/550 2026-04-10 08:37 by 5268321
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 318求调剂 +5 李青山山山 2026-04-07 5/250 2026-04-07 18:24 by 蓝云思雨
信息提示
请填处理意见