| 查看: 2786 | 回复: 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)'); |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有300人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有86人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
30