版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

lhl851128

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 25
帖子: 6
在线: 2.3小时
虫号: 1349442
注册: 2011-07-18
专业: 光学

[求助] 维数问题，帮忙看看错在哪里

clc;clear;
F = zeros(100000,100000);
t=0:0.0000005:0.0000005*99999;
f=0:0.5*pi/99999:0.5*pi;
[T1,Y1]=meshgrid(t,f);
i=sqrt(-1);
po=[.5,0,0,.5;0,0,0,0;0,0,0,0;.5,0,0,.5];
for x = 0:1:99999
pt= [ 1/2*(cos(411/10*t(x+1))-i*sin(411/10*t(x+1)))*exp(i*f(x+1))*conj(cos(411/10*t(x+1))-i*sin(411/10*t(x+1))), 0,0， 1/2*(cos(411/10*t(x+1))i*sin(411/10*t(x+1)))*exp(i*f(x+1))*conj(cos(1193/10*t(x+1))+i*sin(1193/10*t(x+1)));
0, 0, 0, 0;
0, 0, 0, 0;
1/2*(cos(1193/10*t(x+1))+i*sin(1193/10*t(x+1)))*exp(i*f(x+1))*conj(cos(411/10*t(x+1))-i*sin(411/10*t(x+1))), 0,0,1/2*(cos(1193/10*t(x+1))+i*sin(1193/10*t(x+1)))*exp(i*f(x+1))*conj(cos(1193/10*t(x+1))+i*sin(1193/10*t(x+1)))];
b(x+1)=[trace((((po)^.5*pt*((po)^.5)))^.5)]^2;
end
for i=1:100000
F(:,i)=b';
end
F1=interp2(t,f,F,T1,Y1);
figure;
surf(T1,Y1,F);
shading interp;
colormap(jet(111));
axis([0 0.05 0 1.6 0 1]);
view(-60,10);

回复此楼

» 猜你喜欢

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

德国外导原则同意联陪，但有几个问题该怎么回答，大家帮忙看看已经有5人回复
大家帮忙看看VMD安装过程中出现的问题！~ 已经有8人回复
大家帮忙看看这个氢谱存在的问题已经有13人回复
【求助】求UDF高手帮忙看看我的问题已经有18人回复
【求助】fluent的传热问题，计算结果有没有问题，请大家帮忙看看。【已解决】已经有10人回复
小硕面上，请大家帮忙看看专家意见是几个A,几个B,几个C的问题已经有68人回复
一个proof的问题，请帮忙看看已经有9人回复
【求助】欧式距离在很高维空间是否有维数灾难问题(太重要了) 已经有16人回复

1楼 2011-11-17 10:40:39

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

bill_shark

银虫 (小有名气)

应助: 10 (幼儿园)
金币: 1194.7
红花: 4
帖子: 100
在线: 63.3小时
虫号: 1437428
注册: 2011-10-11
专业: 光学

【答案】应助回帖

LZ你行数和列数太多了，点取太密了，我的小电脑可运算不了啊。。。
pt第一行，第三列的0后面的逗号不要用中文的逗号，用英文的；第一行第四列的1/2*(cos(411/10*t(x+1))i*sin(411/10*t(x+1)))*exp(i*f(x+1))*conj(cos(1193/10*t(x+1))+i*sin(1193/10*t(x+1)))，其中1/2*(cos(411/10*t(x+1))与i之间少了个运算符，或者是其他什么？
还有LZ你连要算什么都不讲，用的什么公式，我也不知道你的公式是否写对，要得到什么图形？这些都要说清楚嘛，随便扔个程序上来，没人看得懂，也懒得去猜啊。。。

赞一下