24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wzp0098

木虫 (著名写手)

应助: 2 (幼儿园)
金币: 2524.4
散金: 39
红花: 13
帖子: 1218
在线: 226.6小时
虫号: 853822
注册: 2009-09-22
性别: GG
专业: 光学

[求助] 利用for语句画二维图

K = -10:0.1:10;
P = []; % <-------------------添加此句分配P为向量
for n=1:length(K);
      K1 = K(n);
      w = 10;
      K2 = K1+2*w; x = 0.5;
      r1 = 8; r2 = 0.001*r1; T = 3;
      A=[  -r1 r2-r1  -i*x    i*x    0          0    0             0;
               0 -r2    0    0    0          0    i*T          -i*T;
         -2*i*x  -i*x -r1-i*K1*0.5 0    -i*T       0    0             0;
            2*i*x  i*x    0  -r1+i*K1*0.5 0       i*T    0             0;
               0    0    -i*T    0 -r2-i*(K1+K2)*0.5 0    i*x          0;
               0    0       0    i*T    0 -r2+i*(K1+K2)*0.5 0          -i*x;
               i*T 2*i*T    0    0    i*x       0 -r2-r1+i*(K1-K2)*0.5 0;
            -i*T -2*i*T    0    0    0       -i*x    0 -r2-r1-i*(K1-K2)*0.5];
      B = [-r1;0;-i*x;i*x;0;0;i*T;-i*T];
      Y = A\B;
P(n) = imag(Y(3)); % <---------------------------P向量中的元素赋值
end
plot(K,P)
该程序可以画一维图形，如果我想画P有关（K1和W）的二维图形，怎么处理？请牛人指教！

回复此楼

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有11人回复
救命帖已经有5人回复
限项规定已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
招博士已经有3人回复
存款400万可以在学校里躺平吗已经有35人回复
最失望的一年已经有18人回复
求推荐博导已经有4人回复
求推荐英文EI期刊已经有5人回复
疑惑？已经有5人回复

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

gambit建模想在一个二维的长方形内画一段横线表示一个木板，该怎么建模呢？已经有8人回复
在matlab画二维图时，图形的纵坐标达到了十的十几次方，太大了已经有8人回复
matlab 画三维怎么画呢，在二维上有颜色渐进表示已经有3人回复
Matlab二维图形给图形加上颜色已经有4人回复
matlab里怎么画两个自变量的二维图已经有3人回复
请教一下，在pcolor画出的二维图形后，怎么样画一条等值（即取一个特定颜色值）线已经有6人回复
求直接模拟蒙特卡罗方法的程序，包括二维和三维的用fortran编写的，谢谢已经有9人回复
AutoCAD画二维图？已经有10人回复
Catia画的二维图怎么导入ICEM呢？已经有13人回复
如何用matlab画二维网格已经有5人回复
新手求助：matlab 画二维图已经有7人回复
求助MATLAB画图，二维曲线，总是出错，求助高手已经有3人回复
画H(theta,r)函数的二维图，其中theta为x轴，H为y轴，r用颜色表示已经有6人回复
这个图怎么画啊，是二维的。。。已经有7人回复
GAMBIT 二维画图中的一个问题已经有5人回复
CCA二维关系排序图怎么画？已经有5人回复
请教Fortran下如何给二维或三维数组赋初值，谢谢已经有3人回复
有什么方法可以知道自己做的晶体是否具有一维、二维和三维结构？用Diamond画图？已经有14人回复
大家看看二维图用什么软件画既好用也比较容易上手？已经有14人回复
【求助】二维快速傅里叶（Fourier）变换FORTRAN程序已经有5人回复
【请教】在Diamond里画一维、二维、三维图【请大家帮忙】已经有4人回复
【求助】matlab 画出类似的二维图已经有3人回复

1楼 2012-07-15 09:41:12

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

wzp0098

木虫 (著名写手)

应助: 2 (幼儿园)
金币: 2524.4
散金: 39
红花: 13
帖子: 1218
在线: 226.6小时
虫号: 853822
注册: 2009-09-22
性别: GG
专业: 光学

引用回帖:

3楼: Originally posted by mei3014 at 2012-07-15 12:39:44
你的w又不是变量，怎么画二维图，除非w=10*ones（1，length（K））。
这样就要再用一个for循环，最后输出的P应该是201*201的矩阵，因为1*201的向量是画不了二维图的！

K = -50:0.1:50;W=-50:0.1:50;
P = [] % <---------添加此句分配P为向量
for n=1:length(K);
      K1 = K(n);
      m=1:length(W);
      w=W(m);
      x = 0.1;
      K2 = K1+2*w;
      r1 = 6.6; r2 = 0.001*r1; T = 3;
      A=[  -r1 r2-r1  -i*x    i*x    0          0    0             0;
               0 -r2    0    0    0          0    i*T          -i*T;
         -2*i*x  -i*x -r1-i*K1*0.5 0    -i*T       0    0             0;
            2*i*x  i*x    0  -r1+i*K1*0.5 0       i*T    0             0;
               0    0    -i*T    0 -r2-i*(K1+K2)*0.5 0    i*x          0;
               0    0       0    i*T    0 -r2+i*(K1+K2)*0.5 0          -i*x;
               i*T 2*i*T    0    0    i*x       0 -r2-r1+i*(K1-K2)*0.5 0;
            -i*T -2*i*T    0    0    0       -i*x    0 -r2-r1-i*(K1-K2)*0.5];
      B = [-r1;0;-i*x;i*x;0;0;i*T;-i*T];
      Y = A\B;
P(n,m) =imag(Y(3)); % <---------P向量中的元素赋值
end
plot3(K,W,P)