- 博客(3)
- 收藏
- 关注
RSS订阅原创 二维边值问题（II）波动方程、动图处理
波动方程①源代码：cleart=linspace(0,100,5000);x=linspace(0,10,500);alpha=1.2;x0=0;c=1; %初始化dx=x(2)-x(1);dt=t(2)-t(1);r=c*(dt/dx);for j=1:2 for i=1:length(x) y(i,j)=sin(i*dx-x0)*exp(-1.2*i*dx); endendfor k=3:length(t) y(1...
2021-10-28 18:12:33
703
原创 二维边值问题的数值方法（I）（例子为泊松方程）雅克比算法和GS算法以及SOR算法优劣的探究
1.微分方程为。（1）分析：对于此类微分方程，我们可以找到一定的物理意义，即有源（或者汇）的静电场满足的微分方程，又叫泊松方程。因此V应该是有物理意义的，即电势。这个方程我们很容易就可以解出它的解之一为，所以我们可以根据这个解给定它的边界条件，其中给定的边界条件为①通过代码来说明：maxgrid=50;V=zeros(maxgrid,maxgrid)；%我们对研究的对象V建立一个50×50的矩阵。y=linspace(-10,10,50);%其中我们让这个矩阵对-10到10进行等间距划分。
2021-10-28 18:10:37
1324
原创 求解第一类边值问题 源：matlab
matlab爱好者。求解第一类边值问题的四种方法（完整），包括相互的验证：任何第一类边值问题都可以用这四种方法快速解决。矩阵法这里只能适用于5个离散点，想得到通解请自己推一遍矩阵法。我们以一道例题开始：题目：解法：①初值解法解边值：通过给dy来求y（end），直到拟合。首先做一个函数，给一个dy,返回一个y（end）。函数名为findf，如下:function s=findf(dy) t=linspace(0,4,11);dt=t(2)-t(1);y(1)=0; v(1)=dy;fo
2021-10-24 23:44:33
1181
1
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人