| 5 | 1/1 | 返回列表 |
| 查看: 1956 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
wzp0098木蟲(chóng) (著名寫(xiě)手)
|
[求助]
利用for 語(yǔ)句畫(huà)圖
|
||
|
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ǔ)句畫(huà)出圖形,可是無(wú)法運(yùn)行,求大蝦賜教! |
木蟲(chóng) (著名寫(xiě)手)
|
謝謝,我就是想畫(huà)出y隨K1變化的圖像!K1取一個(gè)值,就可以得到一個(gè)y值,利用for語(yǔ)句循環(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) |
版主 (文壇精英)
米米
版主 (文壇精英)
米米
木蟲(chóng) (著名寫(xiě)手)
|
根據(jù)你的意見(jiàn),我做了些修改,具體如下, 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之前的部分是矩陣求解過(guò)程。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑推薦 材料 304 +17 | 荷包蛋hyj 2026-03-26 | 17/850 |
|
|---|---|---|---|---|
|
[考研] 土木304求調(diào)劑 +4 | 兔突突突, 2026-03-31 | 5/250 |
|
|
[考研] 266分,求材料相關(guān)專(zhuān)業(yè)調(diào)劑 +10 | 哇呼哼呼哼 2026-03-30 | 12/600 |
|
|
[考研] 0703化學(xué) +19 | 妮妮ninicgb 2026-03-27 | 19/950 |
|
|
[考研] 材料專(zhuān)碩調(diào)劑 +12 | 椰椰。 2026-03-29 | 12/600 |
|
|
[考研] 本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分 +5 | 里子木yy 2026-03-29 | 5/250 |
|
|
[論文投稿]
chinese chemical letters英文版投稿求助
130+3
|
Yishengeryi 2026-03-30 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9 | hanamiko 2026-03-27 | 9/450 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚(yú) 2026-03-27 | 12/600 |
|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
[考研] 283求調(diào)劑(080500) +14 | A child 2026-03-27 | 14/700 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專(zhuān)業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|