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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 利用for语句画二维图
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1330  |  回复: 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的回帖

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的回帖
查看全部 8 个回答

wzp0098

木虫 (著名写手)
画P有关（K1和w）的二维图形，请大虾指教！
赞一下 回复此楼
2楼2012-07-15 09:43:25
已阅   回复此楼   关注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的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +6 202451007219 2026-04-02 6/300 2026-04-03 01:16 by BruceLiu320
[考研] 环境285分，过六级，求调剂 +7 xhr12 2026-04-02 7/350 2026-04-02 22:45 by chran16
[考研] 土木水利328分求调剂 +4 疾风知劲草666 2026-04-02 4/200 2026-04-02 21:59 by barlinike
[考研] 349求调剂 +10 zwjjjjjj 2026-03-31 10/500 2026-04-02 20:13 by dongzh2009
[考研] 求生物学调剂 +10 15172915737 2026-04-01 10/500 2026-04-02 18:53 by 哦哦嗯哈
[考研] 一志愿厦门大学材料工程专硕354找调剂！！！ +7 贝呗钡钡 2026-03-30 7/350 2026-04-02 12:19 by 1753564080
[考研] 085601一志愿中山大学深圳材料工程330求调剂 +8 pipiver 2026-03-30 8/400 2026-04-02 12:01 by ms629
[考研] 292求调剂 +14 是妍子也是研子 2026-03-30 14/700 2026-04-02 10:17 by sanrepian
[考研] 354求调剂 +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[考研] 土木304求调剂 +6 兔突突突， 2026-03-31 7/350 2026-04-02 09:06 by coolminer
[考研] 070300化学279求调剂 +15 哈哈哈^_^ 2026-03-31 17/850 2026-04-01 21:37 by 给你你注意休息
[考研] 320分，材料与化工专业，求调剂 +14 一定上岸aaa 2026-03-27 18/900 2026-04-01 20:10 by 积极调剂的小学�
[考研] 08工科275求调剂，可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 化学0703 调剂 306分 一志愿211 +12 26要上岸 2026-03-28 12/600 2026-04-01 11:10 by chemdavid
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
[考研] 326求调剂 +4 崽崽仔 2026-03-31 4/200 2026-04-01 09:58 by 我的船我的海
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-04-01 09:04 by oooqiao
[考研] 求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
[考研] 346求调剂 一志愿070303有机化学 +11 萝卜炖青菜 2026-03-28 12/600 2026-03-31 17:54 by xhai2011
[考研] 数一英一271专硕（085401）求调剂，可跨 +7 前行必有光 2026-03-28 8/400 2026-03-28 23:22 by 小木虫tim
信息提示
关闭
请填处理意见
关闭