应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台» 网络生活区» 有奖问答 » 有奖问答完结子版 » 求助matlab编程
41/1返回列表
查看: 1017  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

陈昌荣

新虫 (初入文坛)

[求助] 求助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');
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-04-08 10:09:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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）；的右括号成了全角的了
赞一下(2人) 回复此楼
low-key
3楼2013-04-08 13:08:55
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 陈昌荣 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[考研] 材料学求调剂 +6 Stella_Yao 2026-03-20 6/300 2026-03-25 00:37 by baoball
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 材料学硕，求调剂 6+3 糖葫芦888ll 2026-03-22 7/350 2026-03-24 17:11 by hello七七
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 一志愿重庆大学085700资源与环境，总分308求调剂 +7 墨墨漠 2026-03-23 8/400 2026-03-23 20:36 by Creta
[考研] 一志愿东华大学化学070300，求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 289求调剂 +7 怀瑾握瑜l 2026-03-20 7/350 2026-03-22 15:57 by ColorlessPI
[考研] 一志愿南大，0703化学，分数336，求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
信息提示
关闭
请填处理意见
关闭