| 10 | 1/1 | 返回列表 |
| 查看: 1947 | 回復: 9 | ||||
wzp0098木蟲 (著名寫手)
|
[求助]
利用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)語句畫出圖形,可是無法運行,求大蝦賜教! |
版主 (文壇精英)
米米
木蟲 (著名寫手)
|
謝謝,我就是想畫出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) |
版主 (文壇精英)
米米
木蟲 (著名寫手)
|
根據(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之前的部分是矩陣求解過程。 |
版主 (文壇精英)
米米
榮譽版主 (著名寫手)
![]() |
專家經(jīng)驗: +4 |
|
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) |

金蟲 (小有名氣)
木蟲 (著名寫手)
金蟲 (小有名氣)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
求調(diào)劑材料學碩080500,總分289分
5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|---|---|---|---|---|
|
[考研] 070300化學求調(diào)劑 +5 | 苑豆豆 2026-03-20 | 5/250 |
|
|
[考研] 284求調(diào)劑 +6 | Zhao anqi 2026-03-22 | 6/300 |
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +5 |  ̄^ ̄゜汗 2026-03-19 | 7/350 |
|
|
[考研] 293求調(diào)劑 +12 | zjl的號 2026-03-16 | 17/850 |
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 南昌大學材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
[考研]
|
.6lL 2026-03-18 | 8/400 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
[考研] 南京大學化學376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
[考研] 0703化學調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|