| 查看: 2706 | 回复: 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)'); |
回复此楼» 猜你喜欢
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有260人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
急！！！用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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人