 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 仿真模擬 » MATLAB » 【求助】畫三維圖【已解決】
| 5 | 1/1 | 返回列表 |
| 查看: 551 | 回復: 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 ] |
回復此樓» 猜你喜歡
311（085601）求調劑
已經(jīng)有4人回復
-
305求調劑
已經(jīng)有4人回復
-
一志愿北化085600材料專碩275|有文章專利｜求調劑
已經(jīng)有7人回復
-
一志愿華理，數(shù)一英一285求A區(qū)調劑
已經(jīng)有12人回復
-
289求調劑
已經(jīng)有11人回復
-
食品工程專碩一志愿中海洋309求調劑
已經(jīng)有7人回復
-
生物學學碩，一志愿湖南大學，初試成績338
已經(jīng)有6人回復
-
321求調劑
已經(jīng)有7人回復
-
343求調劑
已經(jīng)有4人回復
-
本科新能源科學與工程，一志愿華理能動285求調劑
已經(jīng)有6人回復
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ā)表 | |
|---|---|---|---|---|
|
[考研] 311（085601）求調劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 本科新能源科學與工程，一志愿華理能動285求調劑 +3 | AZMK 2026-03-28 | 6/300 |
|
|
[考研] 070305高分子化學與物理 304分求調劑 +4 | c297914 2026-03-28 | 4/200 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
|
[考研] 一志愿中南大學化學0703總分337求調劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
|
[考研] 085602 307分 求調劑 +7 | 不知道叫什么！ 2026-03-26 | 7/350 |
|
|
|
[考研] 291求調劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
|
[考研] 086000調劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調劑 +10 | 幸運的醬醬 2026-03-22 | 12/600 |
|
|
|
[考研] 279 分 求調劑 +4 | 睡個好覺_16 2026-03-24 | 4/200 |
|
|
|
[考研] 315調劑 +4 | 0860求調劑 2026-03-26 | 5/250 |
|
|
|
[考研] 292求調劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 321求調劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研] 340求調劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
|
[考研] 材料與化工304求B區(qū)調劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
|
[考研] 293求調劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
|
[考研]
 一志愿070300浙大化學358分，求調劑！
+4
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
|
[考研] 293求調劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
|
[考研] 石河子大學（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 306求調劑 +5 | 來好運來來來 2026-03-22 | 5/250 |
|