| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1955 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +12 | Y-cap 2026-03-29 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +8 | 墨墨漠 2026-03-25 | 8/400 |
|
|
[考研] 0817化工學(xué)碩調(diào)劑 +4 | 努力上岸中! 2026-03-31 | 4/200 |
|
|
[考研] 347求調(diào)劑 +6 | 山頂見(jiàn)α 2026-03-25 | 6/300 |
|
|
[考研] 環(huán)境工程 085701,267求調(diào)劑 +11 | minht 2026-03-29 | 11/550 |
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
[考研] 0703一志愿9,初試成績(jī):338,四六級(jí)已過(guò),有科研經(jīng)歷,求調(diào)劑! +7 | Zuhui0306 2026-03-25 | 7/350 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +10 | 景隨風(fēng) 2026-03-29 | 15/750 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 266求調(diào)劑 +11 | 陽(yáng)陽(yáng)哇塞 2026-03-27 | 12/600 |
|
|
[考研] 08開(kāi)頭275求調(diào)劑 +4 | 拉誰(shuí)不重要 2026-03-26 | 4/200 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 296求調(diào)劑 +4 | 汪!?! 2026-03-25 | 7/350 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|