| 查看: 2712 | 回复: 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)'); |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有136人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
找到一些相关的精华帖子，希望有用哦~
急！！！用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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人