 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 仿真模擬 » MATLAB » 【求助】畫三維圖【已解決】
| 5 | 1/1 | 返回列表 |
| 查看: 552 | 回復: 4 | |||
| 當前主題已經(jīng)存檔。 | |||
[交流]
【求助】畫三維圖【已解決】
|
|||
|
function equ1=funb(t,y,flag,A,B,C,D,E,F,K,M) equ1=[y(2);(-E*y(2)-F*y(1)+((A*t+B)./M)*(cosh(C*D*t.^2./2)-cos(C*D*t.^2./2)-K*(sinh(C*D*t.^2./2)-sin(C*D*t.^2./2))))]; function ctn A=1;B=1;C=1;D=1;E=0.1;F=1;K=1;M=2; y0=[0;0]; tspan=[0:0.1:40]; [t,y]=ode45('funb',tspan,y0,[],A,B,C,D,E,F,K,M); 以上我是求的一個微分方程。 接下來，我要畫一個三維圖，Z(x,t)=Y(t)*W(x) W(x)表達式已知，現(xiàn)在我的問題是我不知道該如何把解微分方程得到的數(shù)值解Y(t)加入到畫圖程序中，我試著用過surf and plot3但總是出錯，，，，望高手不吝賜教！ W(x)=(cosh(D*x)-cos(D.*x)-K.*(sinh(D.*x)-sin(D.*x))) [ Last edited by nono2009 on 2010-3-30 at 09:21 ] |
回復此樓» 猜你喜歡
332求調劑
已經(jīng)有5人回復
-
085701求調劑初試286分
已經(jīng)有3人回復
-
材料工程294求調劑
已經(jīng)有4人回復
-
本科新能源科學與工程，一志愿華理能動285求調劑
已經(jīng)有10人回復
-
070305高分子化學與物理 304分求調劑
已經(jīng)有5人回復
-
083000學碩274求調劑
已經(jīng)有8人回復
-
一志愿北京化工大學材料與化工（085600）296求調劑
已經(jīng)有10人回復
-
求調劑一志愿海大，0703化學學碩304分，有大創(chuàng)項目，四級已過
已經(jīng)有11人回復
-
071000生物學求調劑，初試成績343
已經(jīng)有7人回復
-
0703化學調劑，求導師收
已經(jīng)有8人回復
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應助: 44 (小學生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時
- 蟲號: 496517
- 注冊: 2008-01-19
- 專業(yè): 化學反應工程
★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
nono2009(金幣+2):鼓勵應助！ 2010-03-26 18:25
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
nono2009(金幣+2):鼓勵應助！ 2010-03-26 18:25
|
function ctn global A B C D E F K M A=1; B=1; C=1; D=1; E=0.1; F=1; K=1; M=2; y0=[0;0]; tspan=0:0.1:40; [t,y]=ode45(@funb,tspan,y0); x=0:0.1:40; Wx=(cosh(D*x)-cos(D.*x)-K.*(sinh(D.*x)-sin(D.*x))); Z1=y(:,1)*Wx; Z2=y(:,2)*Wx; [T,X]=meshgrid(0:0.1:40,0:0.1:40); figure(1) mesh(T,X,Z1) title('y1') figure(2) mesh(T,X,Z2) title('y2') function equ1=funb(t,y) global A B C D E F K M equ1=[y(2); (-E*y(2)-F*y(1)+((A*t+B)./M)*(cosh(C*D*t.^2./2)-cos(C*D*t.^2./2)-K*(sinh(C*D*t.^2./2)-sin(C*D*t.^2./2))))]; |
2樓2010-03-26 16:58:06
|
謝謝，我后來也編出來了，用的是個循環(huán)。但是圖不如你的好看。。。 function ctn A=1;B=1;C=1;D=1;E=0.1;F=1;K=1;M=2; y0=[0;0]; tspan=[0:1:40]; [t,y]=ode45('funb',tspan,y0,[],A,B,C,D,E,F,K,M); x=0:0.025:1; a=length(t); b=length(x); k=1; for i=1:a for j=1:b z(j,i)=(cosh(D*x(j))-cos(D.*x(j))-K.*(sinh(D.*x(j))-sin(D.*x(j)))).*y(i,1); end end mesh(t,x,z); xlabel('t'),ylabel('x'),zlabel('z') colorbar function equ1=funb(t,y,flag,A,B,C,D,E,F,K,M) equ1=[y(2);(-E*y(2)-F*y(1)+((A*t+B)./M)*(cosh(C*D*t.^2./2)-cos(C*D*t.^2./2)-K*(sinh(C*D*t.^2./2)-sin(C*D*t.^2./2))))]; 還有一個問題，我在保存圖片的時候，總是保存不上，打開之后就是一片黑的，沒有圖。有的說是內存不夠了，但是我覺得這個圖沒有多大。PS：我換了好幾種格式還是保存不上。 |
3樓2010-03-29 09:22:44
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應助: 44 (小學生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時
- 蟲號: 496517
- 注冊: 2008-01-19
- 專業(yè): 化學反應工程
4樓2010-03-29 12:10:18
5樓2010-03-30 09:06:48
相關版塊跳轉
我要訂閱樓主 岸芷汀蘭- 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 0703化學 +10 | 妮妮ninicgb 2026-03-27 | 10/500 |
|
|
|
[考研] 321求調劑 +6 | 材料cailiao 2026-03-21 | 6/300 |
|
|
|
[考研] 299求調劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 291求調劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
|
[考研]
 材料與化工考研調劑
+17
|
孅華 2026-03-22 | 17/850 |
|
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學-環(huán)境工程專碩-298 +4 | 手機用戶 2026-03-26 | 4/200 |
|
|
|
[有機交流]
高溫高壓反應求助
 10+4
|
chibby 2026-03-25 | 4/200 |
|
|
|
[考研] 08開頭275求調劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 322求調劑 +4 | 我真的很想學習 2026-03-23 | 4/200 |
|
|
|
[考研] 材料專碩 335 分求調劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
|
[考研] 26考研-291分-廈門大學（085601）-柔性電子學院材料工程專業(yè)求調劑 +3 | min3 2026-03-24 | 4/200 |
|
|
|
[考研] 347求調劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 284求調劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
|
[考研] 生物學學碩求調劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
|
[考研] 一志愿武理085500機械專業(yè)總分300求調劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
|
[考研] 求調劑一志愿武漢理工大學材料工程（085601） +5 | WW.' 2026-03-23 | 7/350 |
|
|
|
[考研] 335求調劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
|
[考研] 280分求調劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|