| 1 | 1/1 | 返回列表 |
| 查看: 2732 | 回復(fù): 0 | ||
[求助]
用MATLAB寫了一個解擴(kuò)散方程的程序,需要修改一下,求大神指導(dǎo)一下!
|
|
題是:假設(shè)擴(kuò)散方程du/dt=a du^2/dx^2,(這是偏微分方程,應(yīng)為打不出那個符號,用dt,dx代替),取dx=pi/15,u(x,0)=sin(x),u(0,t)=0,u(pi,t)=0 現(xiàn)在想修改成du/dt=a du^2/dx^2+u^3+u^5 其他條件不變,想問一下大神下面的M文件和運(yùn)行命令,該怎么改變? 函數(shù)文件: function y=myfun(~,u) n=15; x=zeros(1,n); x(1)=-2*u(1)+u(2); x(n-1)=u(n-2)-2*u(n-1); for i=2:n-2 x(i)=u(i-1)-2*u(i)+u(i+1); end y=x'; End 運(yùn)行命令: clc;clear all;format long; n=15;dpi=pi/n; x=dpi:dpi:pi; size(x) u=sin(x); t=0.:0.4:40.; [t u]=ode45(@myfun,t,u); uu(:,2:n)=u(:,1:n-1); uu(:,1)=0;uu(:,n+1)=0; x=[0. x];[xx yy]=meshgrid(x,t); surf(xx,yy,uu); xlabel('x');ylabel('t');zlabel('u(x,t)'); |
找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研]
|
.6lL 2026-03-18 | 8/400 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計(jì)算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|