| 查看: 1110 | 回复: 1 | ||
[求助]
用matlab求解一偏微分方程，求助！ 已有1人参与
|
|
在单位正方形上有偏微分方程D^2U/Dx^2+D^2U/Dy^2=2(x^2+y^2)，定解条件是u（0，y）=u(x,0)=0;u(1,y)=y^2,u(x,1)=x^2. 分别取N=5,10,15把单位正方形划分成N*N个小格子，并把格子线的交点编号（i,j）， 记u(x(i),x(j))=u(i,j).用熟知微分公式建立关于u（i,j）的线性代数方程组，并求解。 该方程的精确解是u（x，y）=（xy）^2，将精确解和近似解画在同一张三维图片上 原题如图，实在想不出了，14号就要交了，求大神们帮我看看，最好能帮我打好代码，或者指点一二也是很好的，在下感激不尽！  IMG_20150712_124032_1436676074514.jpg |
回复此楼» 猜你喜欢
高分子化学与物理调剂
已经有8人回复
-
290求调剂
已经有7人回复
-
272求调剂
已经有4人回复
-
291分工科求调剂
已经有5人回复
-
0856求调剂285
已经有5人回复
-
求助coordination chemistry reviews 的写作模板
已经有4人回复
-
Optics letters投稿被拒求助
已经有4人回复
-
284求调剂
已经有5人回复
-
材料类求调剂
已经有7人回复
-
298求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急！！！用matlab求解微分方程的问题
已经有10人回复
- 求助二阶偏微分方程求解
已经有13人回复
- 用Matlab求一个二阶微分方程，出错了，求指导（求助）
已经有13人回复
- 偏微分方程Matlab求解
已经有7人回复
- 求助：用Matlab解微分方程，命令怎么编辑？有例子
已经有13人回复
- 变系数偏微分方程matlab求解
已经有6人回复
- 如何用Matlab求时变微分方程的解？
已经有5人回复
- matlab解偏微分方程求助
已经有12人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 偏微分方程数值计算书籍推荐
已经有10人回复
- 分数阶微分方程的求解（用差分法，matlab编程求解）
已经有9人回复
- 求助：如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？
已经有5人回复
- matlab求解偏微分方程组
已经有3人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- Matlab求解偏微分方程组
已经有22人回复
- 请各位前辈指点一下，如何用MATLAB求解如下的微分方程，谢谢！
已经有8人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】求助PDEPE函数求解高阶偏微分方程 ：D
已经有5人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】求助用matlab解偏微分方程
已经有15人回复
zhangzhiguang
至尊木虫 (文坛精英)
- 应助: 121 (高中生)
- 金币: 53088.8
- 红花: 24
- 沙发: 222
- 帖子: 26519
- 在线: 881.4小时
- 虫号: 492270
- 注册: 2008-01-08
- 性别: GG
- 专业: 基础物理学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
faker_: 金币+50, ★★★★★最佳答案, 十分感谢，这么辛苦的帮我打好了，晚上我好好看看，再次感谢！！ 2015-07-13 19:36:12
感谢参与，应助指数 +1
faker_: 金币+50, ★★★★★最佳答案, 十分感谢，这么辛苦的帮我打好了，晚上我好好看看，再次感谢！！ 2015-07-13 19:36:12
|
你这个应该是学习有限差分法的练习题吧。敲字太麻烦，我给你手写了一个思路，你自己做一下吧。求解线性方程组有很多方法，像左除，迭代等，我采用迭代处理了十个网格的，你可以借鉴。列线性方程如图所示。迭代求解及画图代码如下： 本程序为有限差分法解偏微分方程 clear;clc; h=1/10; u=zeros(11); for n=1:10 for i=2:10 for j=2:10 u(11,j)=(h*j)^2; u(i,11)=(h*i)^2; u(11,11)=(h*(i+1))^2; u(i,j)=1/4*(u(i+1,j)+u(i-1,j)+u(i,j+1)+u(i,j-1)+2*h^2*((h*i)^2+(h*j)^2)); end end end x=0:0.1:1; y=0:0.1:1; [X,Y]=meshgrid(x,y); U=X.^2.*Y.^2; mesh(X,Y,u) hold on mesh(X,Y,U)  IMG_20150713_190353.jpg |
2楼2015-07-13 19:11:10
30