| 查看: 2743 | 回复: 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)'); |
回复此楼» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有20人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有269人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有9人回复
-
桂林理工大学物理学专业招收调剂
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
急！！！用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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
