| 查看: 387 | 回复: 2 | ||
[求助]
画出来 的图只有空格子没有图形，帮我看看问题在哪里
|
|
clc;clear; F = zeros(1,100); for t = 0:0.05:0.05*99; po=[1/2,0,0,1/2; 0,0,0,0; 0,0,0,0; 1/2,0,0,1/2] pt=[ (1/2*cos(1781/10*t)-1/2*i*sin(1781/10*t))*conj(cos(1781/10*t)-i*sin(1781/10*t)),0,0,(1/2*cos(1781/10*t)-1/2*i*sin(1781/10*t))*conj(cos(2563/10*t)+i*sin(2563/10*t)); 0,0,0,0; 0,0,0,0; (1/2*cos(2563/10*t)+1/2*i*sin(2563/10*t))*conj(cos(1781/10*t)-i*sin(1781/10*t)),0,0,(1/2*cos(2563/10*t)+1/2*i*sin(2563/10*t))*conj(cos(2563/10*t)+i*sin(2563/10*t))] F=[Trace(sqrt((sqrt(po)*pt*(sqrt(po)))))]^2; end plot(t,F) |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 大家帮我看看图片中的黄色的支撑哪有卖的，万谢
已经有5人回复
- 大家帮我看看图片中的黄色的支撑哪有卖的，万谢
已经有4人回复
- 麻烦塔菲尔曲线高手帮我看看这几个图，说明什么问题，或者能从中得出什么信息？
已经有12人回复
- 请各位专家，帮我看看图片中是什么蕨类植物
已经有4人回复
- 求高手帮我看看DGGE图片，为什么跑不开，非常感谢
已经有32人回复
- 哪位高手帮我看下这种三维图是用什么软件画的？
已经有11人回复
- 谁能帮我看下这张图是用什么软件画出来的？悬赏金币3枚
已经有5人回复
- 哪位大侠帮我看看这个能级图是如何画出来的，不是能带图
已经有3人回复
- 【求助】求助大家帮我看看这个核磁图
已经有9人回复
- 【求助】【求助】请大家帮我看看，图片里溶胶的颜色是银溶胶的颜色吗
已经有4人回复
- 大家帮我看看这个图是怎么画出来的
已经有12人回复
- 【求助/交流】CTAB法提取了dna，帮我看看图片，帮忙提点改进方法
已经有20人回复
- 大家帮我看看cDNA图片对吗？
已经有5人回复
- 【求助】请大家帮我看看，我的循环伏安图说明了什么问题？
已经有7人回复
- 【求助】大家帮我看看怎么画这个图
已经有6人回复
- 【求助】请大家帮我看看下面这个图用fluent模拟的时候是怎么画出来的？
已经有7人回复
- 大家帮我看看这样的图是怎么画出来的？
已经有6人回复
ygEric
金虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 4801.4
- 散金: 128
- 红花: 15
- 帖子: 1421
- 在线: 891.7小时
- 虫号: 1233602
- 注册: 2011-03-15
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
|
clc;clear; F = zeros(1,100); k=1; for t = 0:0.05:0.05*99; po=[1/2,0,0,1/2; 0,0,0,0; 0,0,0,0; 1/2,0,0,1/2] pt=[ (1/2*cos(1781/10*t)-1/2*i*sin(1781/10*t))*conj(cos(1781/10*t)-i*sin(1781/10*t)),0,0,(1/2*cos(1781/10*t)-1/2*i*sin(1781/10*t))*conj(cos(2563/10*t)+i*sin(2563/10*t)); 0,0,0,0; 0,0,0,0; (1/2*cos(2563/10*t)+1/2*i*sin(2563/10*t))*conj(cos(1781/10*t)-i*sin(1781/10*t)),0,0,(1/2*cos(2563/10*t)+1/2*i*sin(2563/10.*t))*conj(cos(2563/10.*t)+i*sin(2563/10*t))] F(k)=[Trace(sqrt((sqrt(po)*pt*(sqrt(po)))))]^2; k=k+1; end plot(t,F) －－－－－－－－－－－－－ F应该是向量，所以要用向量赋值的形式F(k)=? |
2楼2011-11-11 22:41:54
ygEric
金虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 4801.4
- 散金: 128
- 红花: 15
- 帖子: 1421
- 在线: 891.7小时
- 虫号: 1233602
- 注册: 2011-03-15
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
lhl851128(金币+1): 万分感谢 2011-11-14 11:16:07
|
你的程序写错了，上面的也得不出正确结果，下面的为正解： －－－－－－－－－－－－－－－－－－－－ clc;clear; for k = 0:99; t=k*0.05;% t=0:0.05:0.05*99 po=[1/2,0,0,1/2; 0,0,0,0; 0,0,0,0; 1/2,0,0,1/2] pt=[ (1/2*cos(1781/10*t)-1/2*i*sin(1781/10*t))*conj(cos(1781/10*t)-i*sin(1781/10*t)),0,0,(1/2*cos(1781/10*t)-1/2*i*sin(1781/10*t))*conj(cos(2563/10*t)+i*sin(2563/10*t)); 0,0,0,0; 0,0,0,0; (1/2*cos(2563/10*t)+1/2*i*sin(2563/10*t))*conj(cos(1781/10*t)-i*sin(1781/10*t)),0,0,(1/2*cos(2563/10*t)+1/2*i*sin(2563/10.*t))*conj(cos(2563/10.*t)+i*sin(2563/10*t))] F(k+1)=[Trace(sqrt((sqrt(po)*pt*(sqrt(po)))))]^2; % F(1),F(2),F(3),... end x=0:0.05:0.05*99; plot(x,F) axis([0 5.5 -1 6]) 你的错误原因：for中的t是单个变量，不是向量，因此只有一个值。plot时需要x向量和y向量 |
3楼2011-11-11 22:53:58
20