| 8 | 1/1 | 返回列表 |
| 查看: 1313 | 回復(fù): 7 | ||
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有關(guān)(K1和W)的二維圖形,怎么處理?請牛人指教! |
木蟲 (著名寫手)
金蟲 (正式寫手)
木有頭屑

木蟲 (著名寫手)
|
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) |
木蟲 (著名寫手)
木蟲 (著名寫手)
金蟲 (正式寫手)
木有頭屑
|
1)你的循環(huán)次數(shù)太多,循環(huán)里面又是矩陣運(yùn)算,要是將for語句也改成矩陣形式很難,建議減少向量長度或增大間距。 2)i是復(fù)數(shù)吧,matlab里面最好用1i表示,防止與變量i混淆。 3)程序我暫時(shí)改成如下所示,是不是你期待的結(jié)果我也不知道! clc;clear K = -50:1:50;W=-50:1:50; P = []; % <---------添加此句分配P為向量 for n=1:length(K); K1 = K(n); for 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 -1i*x 1i*x 0 0 0 0; 0 -r2 0 0 0 0 1i*T -1i*T; -2*1i*x -1i*x -r1-1i*K1*0.5 0 -1i*T 0 0 0; 2*1i*x 1i*x 0 -r1+1i*K1*0.5 0 1i*T 0 0; 0 0 -1i*T 0 -r2-1i*(K1+K2)*0.5 0 1i*x 0; 0 0 0 1i*T 0 -r2+1i*(K1+K2)*0.5 0 -1i*x; 1i*T 2*1i*T 0 0 1i*x 0 -r2-r1+1i*(K1-K2)*0.5 0; -1i*T -2*1i*T 0 0 0 -1i*x 0 -r2-r1-1i*(K1-K2)*0.5]; B = [-r1;0;-1i*x;1i*x;0;0;1i*T;-1i*T]; Y = A\B; P(m,n) =imag(Y(3)); % <---------P向量中的元素賦值 end end [X,Y]=meshgrid(K,W); plot3(X,Y,P) xlabel('K') ylabel('W') zlabel('P') |

金蟲 (正式寫手)
木有頭屑

| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +4 | 梁梁校校 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[論文投稿]
申請回稿延期一個月,編輯同意了。但系統(tǒng)上的時(shí)間沒變,給編輯又寫郵件了,沒回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 生物學(xué)調(diào)劑招人。! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +17 | 不知道叫什么! 2026-03-15 | 19/950 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 268求調(diào)劑 +6 | 簡單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[基金申請] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|