| 查看: 2719 | 回复: 0 | |||
[求助]
用MATLAB写了一个解扩散方程的程序,需要修改一下，求大神指导一下！
|
|
题是：假设扩散方程du/dt=a du^2/dx^2,(这是偏微分方程，应为打不出那个符号，用dt,dx代替），取dx=pi/15,u(x,0)=sin(x),u(0,t)=0,u(pi,t)=0 现在想修改成du/dt=a du^2/dx^2+u^3+u^5 其他条件不变，想问一下大神下面的M文件和运行命令，该怎么改变？ 函数文件： 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 运行命令: 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)'); |
回复此楼» 猜你喜欢
投稿chemical physical letters不送审？
已经有2人回复
-
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有269人回复
-
杂志要求提供official language-editing certificate
已经有4人回复
-
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有16人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
-
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律
已经有1人回复
-
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
找到一些相关的精华帖子，希望有用哦~
急！！！用matlab求解微分方程的问题
已经有10人回复
- 求大神解下这个微分方程，求出y关于x的函数关系，，谢谢了！
已经有9人回复
- 刚刚接触matlab，需要编写三阶矩阵，但其中有方程，不知道该怎么处理，求大神指导
已经有4人回复
- 求助一个matlab解半监督SVM的程序实例，或者正文如图的matlab求解代码
已经有7人回复
- 用Matlab求一个二阶微分方程，出错了，求指导（求助）
已经有13人回复
- 病态矩阵的精确解求解
已经有17人回复
- 用MATLAB回归方程的参数，是个隐函数方程
已经有10人回复
- 急！求大神用Matlab四阶龙格库塔解个方程！
已经有5人回复
- 求大神用matlab给我处理一组数据，有人没有啊？
已经有15人回复
- 关于matlab进行傅里叶变换和逆变换的程序，求大神指导！
已经有9人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- Matlab求解偏微分方程组
已经有22人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人