版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

wzp0098

木蟲 (著名寫手)

應(yīng)助: 2 (幼兒園)
金幣: 2524.4
散金: 39
紅花: 13
帖子: 1218
在線: 226.6小時
蟲號: 853822
注冊: 2009-09-22
性別: GG
專業(yè): 光學(xué)

[求助] 利用for 語句畫圖

K=-10:1:10;
for m=1:K;
K1=Length(K)
w=10;
K2=K1+2*w;x=0.2;
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;
y=imag(Y(3));
end
plot(K,y)
我想利用for循環(huán)語句畫出圖形，可是無法運行，求大蝦賜教！

回復(fù)此樓

» 猜你喜歡

291 求調(diào)劑已經(jīng)有5人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有11人回復(fù)
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑已經(jīng)有7人回復(fù)
316求調(diào)劑已經(jīng)有7人回復(fù)
材料與化工考研調(diào)劑已經(jīng)有4人回復(fù)
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑已經(jīng)有7人回復(fù)
一志愿211 初試270分求調(diào)劑已經(jīng)有4人回復(fù)
08工學(xué)調(diào)劑已經(jīng)有11人回復(fù)
350求調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）已經(jīng)有4人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

這個用matlab怎么畫已經(jīng)有10人回復(fù)
急：求用Fortran或者其他語言編寫一個程序，處理下面的數(shù)據(jù)�。。。。。。。。。� 已經(jīng)有37人回復(fù)
如何用matlab將這個方程代表的圖畫出來，最好能把畫圖的代碼貼出來，謝謝�。� 已經(jīng)有11人回復(fù)
做研究需要學(xué)編程嗎。學(xué)什么語言好已經(jīng)有20人回復(fù)
前處理軟件如何生成網(wǎng)格已經(jīng)有8人回復(fù)
如何利用Matlab畫圖呢？已經(jīng)有7人回復(fù)
C語言編程，以便畫圖。已經(jīng)有5人回復(fù)
擬合方程并作圖已經(jīng)有6人回復(fù)
向各位老師求助三維畫圖坐標(biāo)軸修改及隱去部分圖的問題已經(jīng)有7人回復(fù)
matlab畫圖問題: 怎樣用小圖標(biāo)或自定義圖形代替+ x . o等做圖已經(jīng)有10人回復(fù)
fortran程序運行結(jié)果輸出及繪圖已經(jīng)有9人回復(fù)
如何用MATLAB求解脈沖微分方程并畫圖已經(jīng)有5人回復(fù)
求助Matlab批處理for語句循環(huán) 已經(jīng)有5人回復(fù)
請教mathematic問題？關(guān)于解微分方程組的和畫圖的一些問題。已經(jīng)有19人回復(fù)
matlab中遇到的畫圖小問題，程序沒問題。標(biāo)題要長@@@@@@@@@@@@@@@@@@@@@@ 已經(jīng)有8人回復(fù)
matlab中如何通過for語句得到已有矩陣的部分?jǐn)?shù)據(jù) 已經(jīng)有3人回復(fù)
想用C++畫一個正方形，該怎么做已經(jīng)有20人回復(fù)
matlab畫圖求助已經(jīng)有11人回復(fù)
MCNP5 畫圖的問題已經(jīng)有7人回復(fù)
VC的畫圖框架中，誰用的比較多？已經(jīng)有17人回復(fù)
【原創(chuàng)】在WORD中使用Mathematica畫的圖形的體會已經(jīng)有13人回復(fù)
【求助】Mathematica畫圖問題已經(jīng)有7人回復(fù)
【求助】再次關(guān)于Mathematica畫圖顏色問題（急用）已經(jīng)有3人回復(fù)
【分享】將Mathematica繪圖數(shù)據(jù)轉(zhuǎn)換至Origin軟件已經(jīng)有134人回復(fù)
【求助】gnuplot畫圖軟件【已完結(jié)】已經(jīng)有15人回復(fù)

1樓 2012-07-07 16:15:25

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fjtony163

版主 (文壇精英)

米米

應(yīng)助: 361 (碩士)
貴賓: 6.876
金幣: 70610.6
散金: 36050
紅花: 222
沙發(fā): 616
帖子: 21912
在線: 3793.7小時
蟲號: 1397876
注冊: 2011-09-11
專業(yè): 遺傳學(xué)研究新技術(shù)與方法
管轄: 論文翻譯

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝參與 2012-07-08 09:54:21
wzp0098: 金幣+5, ★有幫助 2012-07-08 11:11:48

K 是一個行向量，m定義為1到一個行向量肯定運行不了啊。m指的是循環(huán)次數(shù)，應(yīng)該就是m=1：10這樣的。能稍微具體說下是個怎樣的式子么。

贊一下

回復(fù)此樓

2樓2012-07-08 09:19:50

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wzp0098

木蟲 (著名寫手)

應(yīng)助: 2 (幼兒園)
金幣: 2524.4
散金: 39
紅花: 13
帖子: 1218
在線: 226.6小時
蟲號: 853822
注冊: 2009-09-22
性別: GG
專業(yè): 光學(xué)

引用回帖:

2樓: Originally posted by fjtony163 at 2012-07-08 09:19:50
K 是一個行向量，m定義為1到一個行向量肯定運行不了啊。m指的是循環(huán)次數(shù)，應(yīng)該就是m=1：10這樣的。能稍微具體說下是個怎樣的式子么。

謝謝，我就是想畫出y隨K1變化的圖像！K1取一個值，就可以得到一個y值，利用for語句循環(huán)。做了修改，還是不行！
K1=-10:1:10;
for m=1:Length(K);
K1=m
w=10;
K2=K1+2*w;x=0.2;
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;
y=imag(Y(3));
end
plot(K1,y)

贊一下

回復(fù)此樓

3樓2012-07-08 11:17:39

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fjtony163

版主 (文壇精英)

米米

應(yīng)助: 361 (碩士)
貴賓: 6.876
金幣: 70610.6
散金: 36050
紅花: 222
沙發(fā): 616
帖子: 21912
在線: 3793.7小時
蟲號: 1397876
注冊: 2011-09-11
專業(yè): 遺傳學(xué)研究新技術(shù)與方法
管轄: 論文翻譯

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★
臭水溝: 金幣+2, 謝謝交流~~ 2012-07-08 18:02:36
wzp0098: 金幣+5, ★有幫助, 謝謝 2012-07-09 10:30:27

引用回帖:

3樓: Originally posted by wzp0098 at 2012-07-08 11:17:39
謝謝，我就是想畫出y隨K1變化的圖像！K1取一個值，就可以得到一個y值，利用for語句循環(huán)。做了修改，還是不行！
K1=-10:1:10;
for m=1:Length(K);
K1=m
w=10;
K2=K1+2*w;x=0.2;
r1=8;r2=0.001*r1;T=3 ...

從中間的賦值到end之前的部分我不知道是啥，但只要你認(rèn)為輸對了就沒事。

關(guān)鍵是賦值之前的部分，就我的理解應(yīng)該是
K=-10：1：10
for m=1:length(K)
K1=K(m)
w=10;
………………

贊一下

回復(fù)此樓

4樓2012-07-08 11:30:10

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wzp0098

木蟲 (著名寫手)

應(yīng)助: 2 (幼兒園)
金幣: 2524.4
散金: 39
紅花: 13
帖子: 1218
在線: 226.6小時
蟲號: 853822
注冊: 2009-09-22
性別: GG
專業(yè): 光學(xué)

引用回帖:

4樓: Originally posted by fjtony163 at 2012-07-08 11:30:10
從中間的賦值到end之前的部分我不知道是啥，但只要你認(rèn)為輸對了就沒事。

關(guān)鍵是賦值之前的部分，就我的理解應(yīng)該是
K=-10：1：10
for m=1:length(K)
K1=K(m)
w=10;
………………...

根據(jù)你的意見，我做了些修改，具體如下，
K=-10:0.1:10;
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=imag(Y(3));
end
plot(K,P)
但是得不到我想要的圖形，望賜教！
從中間的賦值到end之前的部分是矩陣求解過程。

贊一下

回復(fù)此樓

5樓2012-07-09 10:29:59

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fjtony163

版主 (文壇精英)

米米

應(yīng)助: 361 (碩士)
貴賓: 6.876
金幣: 70610.6
散金: 36050
紅花: 222
沙發(fā): 616
帖子: 21912
在線: 3793.7小時
蟲號: 1397876
注冊: 2011-09-11
專業(yè): 遺傳學(xué)研究新技術(shù)與方法
管轄: 論文翻譯

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
wzp0098: 金幣+5, 謝謝 2012-07-09 11:16:55

引用回帖:

5樓: Originally posted by wzp0098 at 2012-07-09 10:29:59
根據(jù)你的意見，我做了些修改，具體如下，
K=-10:0.1:10;
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=;
B=;
Y=A\B;
P=imag(Y(3));
end
plot(K,P)
但是得不 ...

中間求解過程你確定對了么。。。我是弄不太清楚求解的東西。

贊一下

回復(fù)此樓

6樓2012-07-09 10:58:10

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗: +4
仿真EPI: 6
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
wzp0098: 金幣+30, ★★★很有幫助, 牛人��！ 2012-07-09 11:16:29
臭水溝: 金幣+3, 謝謝交流~~ 2012-07-09 21:24:09

function parafit
K=-10:0.1:10;
P=0*K;
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));
end
plot(K,P)

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

7樓2012-07-09 11:05:11

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

duotojh

金蟲 (小有名氣)

應(yīng)助: 9 (幼兒園)
金幣: 1066.6
紅花: 5
帖子: 225
在線: 137.8小時
蟲號: 1290636
注冊: 2011-05-09
專業(yè): 生物物理、生物化學(xué)與分子

★ ★
臭水溝: 金幣+2, 謝謝交流~~ 2012-07-09 21:24:26

你的意圖應(yīng)該是每輪循環(huán)算出的P連接成曲線，但程序中保留的是最后一輪P值（單值），
因此畫出的圖形是單點當(dāng)然看不到（除非將單點作標(biāo)記才可以看到，例如三角符標(biāo)記等等）
只需添加P向量化，并修改循環(huán)語句中的P賦值這一句即可，修改程序如下：
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)
用freemat畫出的圖形如下

贊一下

回復(fù)此樓

8樓2012-07-09 13:06:49

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wzp0098

木蟲 (著名寫手)

應(yīng)助: 2 (幼兒園)
金幣: 2524.4
散金: 39
紅花: 13
帖子: 1218
在線: 226.6小時
蟲號: 853822
注冊: 2009-09-22
性別: GG
專業(yè): 光學(xué)

引用回帖:

7樓: Originally posted by dbb627 at 2012-07-09 11:05:11
function parafit
K=-10:0.1:10;
P=0*K;
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=;
B=;
Y=A\B;
P(n)=imag(Y(3));
end
plot(K,P)...

我如果想畫出P有關(guān)（K1，w）的二維圖像，怎么處理？請牛人指教！謝謝！

贊一下

回復(fù)此樓

9樓2012-07-09 18:24:16

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

duotojh

金蟲 (小有名氣)

應(yīng)助: 9 (幼兒園)
金幣: 1066.6
紅花: 5
帖子: 225
在線: 137.8小時
蟲號: 1290636
注冊: 2011-05-09
專業(yè): 生物物理、生物化學(xué)與分子

【答案】應(yīng)助回帖

★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝應(yīng)助 2012-07-10 11:13:40
wzp0098: 金幣+2, 請給出程序 2012-07-11 09:57:41

引用回帖:

9樓: Originally posted by wzp0098 at 2012-07-09 18:24:16
我如果想畫出P有關(guān)（K1，w）的二維圖像，怎么處理？請牛人指教！謝謝！...

將 plot 用 scatter 替換即可

贊一下

回復(fù)此樓

10樓2012-07-10 10:51:42

已閱關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 wzp0098 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 材料與化工考研調(diào)劑 +4	孅華 2026-03-22	4/200	2026-03-23 16:13 by 一休哥FU
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +6	墨墨漠 2026-03-23	7/350	2026-03-23 16:09 by 墨墨漠
[考研] 工科0856求調(diào)劑 +4	沐析汀汀 2026-03-21	4/200	2026-03-23 13:08 by 醉在風(fēng)里
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3	13756423260 2026-03-22	3/150	2026-03-22 20:11 by edmund7
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 346求調(diào)劑[0856] +4	WayneLim327 2026-03-16	7/350	2026-03-21 04:02 by JourneyLucky
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 化學(xué)求調(diào)劑 +4	臨澤境llllll 2026-03-17	5/250	2026-03-21 02:23 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 一志愿武理材料305分求調(diào)劑 +6	想上岸的鯉魚 2026-03-18	7/350	2026-03-21 01:03 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計算機(jī)專業(yè)調(diào)劑） +3	123阿標(biāo) 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷