| 查看: 318 | 回复: 1 | ||
[求助]
维数问题 ，帮忙看看错在哪里
|
|
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); |
回复此楼» 猜你喜欢
26应届毕业生考博求助
已经有3人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有7人回复
-
教学课件你会给同学吗
已经有6人回复
-
求助2，4-二氯-5-嘧啶甲醛的合成方法
已经有5人回复
-
本人最近太闲了，谁有问题可以提，每天会统一回复
已经有19人回复
-
申博自荐
已经有6人回复
-
这年头没有找到涵评专家，还有中面上的可能吗
已经有8人回复
-
材料类只有一篇综述能申博么
已经有4人回复
-
风把牡丹吹跑了
已经有3人回复
-
请问大佬b0816评完了吗
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 德国外导原则同意联陪，但有几个问题该怎么回答，大家帮忙看看
已经有5人回复
- 大家帮忙看看VMD安装过程中出现的问题！~
已经有8人回复
- 大家帮忙看看这个氢谱存在的问题
已经有13人回复
- 【求助】求UDF高手帮忙看看我的问题
已经有18人回复
- 【求助】fluent的传热问题，计算结果有没有问题，请大家帮忙看看。【已解决】
已经有10人回复
- 小硕面上，请大家帮忙看看专家意见是几个A,几个B,几个C的问题
已经有68人回复
- 一个proof的问题，请帮忙看看
已经有9人回复
- 【求助】欧式距离在很高维空间是否有维数灾难问题(太重要了)
已经有16人回复
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你连要算什么都不讲，用的什么公式，我也不知道你的公式是否写对，要得到什么图形？这些都要说清楚嘛，随便扔个程序上来，没人看得懂，也懒得去猜啊。。。 |
2楼2011-11-18 22:40:50
20
精华III:
