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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 利用for语句画二维图
北京石油化工学院2026年研究生招生接收调剂公告
81/1返回列表
查看: 1326  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wzp0098

木虫 (著名写手)

[求助] 利用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）的二维图形，怎么处理？请牛人指教！
回复此楼

» 猜你喜欢

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

wzp0098

木虫 (著名写手)
画P有关（K1和w）的二维图形，请大虾指教！
赞一下 回复此楼
2楼2012-07-15 09:43:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mei3014

金虫 (正式写手)

木有头屑

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
臭水沟: 金币+2, 谢谢交流~~ 2012-07-15 19:24:27
你的w又不是变量，怎么画二维图，除非w=10*ones（1，length（K））。
这样就要再用一个for循环，最后输出的P应该是201*201的矩阵，因为1*201的向量是画不了二维图的！
赞一下 回复此楼
一勤天下无难事，处世无难忍最高。
3楼2012-07-15 12:39:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzp0098

木虫 (著名写手)
引用回帖:
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)
赞一下 回复此楼
4楼2012-07-15 14:56:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzp0098

木虫 (著名写手)
引用回帖:
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的向量是画不了二维图的！

在这里  
w ;      
K2 ；
K1 ;
x ;
T;
都可以是变量的，但是我没办法画出二维图，下的程序运行不起来。求指教！谢谢！
赞一下 回复此楼
5楼2012-07-15 15:00:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzp0098

木虫 (著名写手)
如果有大虾能够解决这个问题，将翻倍酬谢！
赞一下 回复此楼
6楼2012-07-15 15:02:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mei3014

金虫 (正式写手)

木有头屑

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
wzp0098: 金币+100, ★★★很有帮助, 非常谢谢！受教了！高人 2012-07-16 12:47:54
xiegangmai: 金币+3, 谢谢应助！ 2012-07-16 23:48:59
引用回帖:
4楼: Originally posted by wzp0098 at 2012-07-15 14:56:28
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;
   ...

1）你的循环次数太多，循环里面又是矩阵运算，要是将for语句也改成矩阵形式很难，建议减少向量长度或增大间距。
2）i是复数吧，matlab里面最好用1i表示，防止与变量i混淆。
3）程序我暂时改成如下所示，是不是你期待的结果我也不知道！
clc;clear
K = -50:1:50;W=-50:1:50;
P = []; % <---------添加此句分配P为向量
for n=1:length(K);
        K1 = K(n);
     for   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  -1i*x      1i*x      0           0       0              0;
                  0    -r2       0       0       0           0      1i*T           -1i*T;
           -2*1i*x  -1i*x -r1-1i*K1*0.5 0     -1i*T          0       0              0;
                2*1i*x  1i*x       0  -r1+1i*K1*0.5 0          1i*T      0              0;
                  0     0      -1i*T      0 -r2-1i*(K1+K2)*0.5 0      1i*x             0;
                  0     0        0      1i*T      0 -r2+1i*(K1+K2)*0.5 0             -1i*x;
                 1i*T 2*1i*T       0       0      1i*x          0 -r2-r1+1i*(K1-K2)*0.5 0;
                -1i*T -2*1i*T      0       0       0         -1i*x      0 -r2-r1-1i*(K1-K2)*0.5];
        B = [-r1;0;-1i*x;1i*x;0;0;1i*T;-1i*T];
        Y = A\B;
    P(m,n) =imag(Y(3)); % <---------P向量中的元素赋值
     end
end
[X,Y]=meshgrid(K,W);
plot3(X,Y,P)
xlabel('K')
ylabel('W')
zlabel('P')
赞一下(1人) 回复此楼
一勤天下无难事，处世无难忍最高。
7楼2012-07-16 11:28:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mei3014

金虫 (正式写手)

木有头屑
★
xiegangmai: 金币+1, 辛苦了！ 2012-07-16 23:49:08
[X,Y]=meshgrid(K,W);
plot3(X,Y,P)
以上两条还是改成
[X1,Y1]=meshgrid(K,W);
plot3(X1,Y1,P)
赞一下 回复此楼
一勤天下无难事，处世无难忍最高。
8楼2012-07-16 11:30:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wzp0098 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 295求调剂 +6 愿旅途永远坦然 2026-04-02 6/300 2026-04-02 21:57 by macy2011
[基金申请] esi高被引论文是不是能对中标有所加分和帮助呢 +4 redcom 2026-04-01 5/250 2026-04-02 20:48 by redcom
[考研] 296材料专硕求调剂 +15 202451007219 2026-04-02 16/800 2026-04-02 20:35 by dongzh2009
[考研] 环境工程297分求调剂一志愿杭高院 +15 GENJIOW 2026-03-31 16/800 2026-04-02 17:56 by cyh—315
[考研] 272求调剂，接受跨专业调剂！ +4 闲鱼卢 2026-03-31 4/200 2026-04-02 11:18 by guyan1000
[考研] 385分 生物学（071000）求调剂 +6 qf626 2026-04-01 6/300 2026-04-02 10:39 by guoweigw
[考研] 070300化学求调剂 +14 小黄鸭宝 2026-03-30 14/700 2026-04-02 10:07 by oooqiao
[考研] 304求调剂 +12 素年祭语 2026-03-31 15/750 2026-04-01 22:41 by peike
[考研] 一志愿南昌大学324求调剂 +12 hanamiko 2026-03-27 12/600 2026-04-01 22:34 by 运气yunqi
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 320分，材料与化工专业，求调剂 +14 一定上岸aaa 2026-03-27 18/900 2026-04-01 20:10 by 积极调剂的小学�
[考研] 调剂 +3 好好读书。 2026-04-01 3/150 2026-04-01 17:06 by zhouyuwinner
[考研] 材料专业调剂 +5 啦啦啦哭 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[考研] 291求调剂 +3 迷蒙木木 2026-04-01 4/200 2026-04-01 11:07 by 逆水乘风
[考研] 求调剂，一志愿北林食品与营养095500，301分，已过六级，有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[考研] 生物考研337分求调剂 +4 cgxin 2026-03-30 6/300 2026-03-31 14:18 by 记事本2026
[考研] 吉大生物学326分求调剂 +3 sunnyupup 2026-03-31 3/150 2026-03-31 09:28 by longlotian
[考研] 105500药学求调剂，一志愿山东大学药学，348分 +3 gr哈哈哈 2026-03-28 3/150 2026-03-30 18:56 by 源_2020
[考研] 材料专硕 085600求调剂 +7 BBQ233 2026-03-30 7/350 2026-03-30 17:44 by oooqiao
[考研] 332求调剂 +6 @MZB382400 2026-03-28 6/300 2026-03-30 16:57 by 无际的草原
信息提示
关闭
请填处理意见
关闭