matlab求抛物问题,高手进，用matlab求解一维抛物方程的初边值问题 - 数学 - 小木虫

下面这个程序我运行过，可以了

%Finite different method for 1D parabolic equation

%Backward Euler for Time, Center Different for Space

function u=parabolicFD(m,n)

h=1/m;%space size

tao=1/n;%time step

r=tao/h^2;

u0=zeros(m+1,1);

u1=zeros(m+1,1);

A=zeros(m+1,m+1);

F=zeros(m+1,1);

for j=1:m+1

u0(j)=cos(pi*(j-1)*h);

end

%assembling the Matrix A

for i=1:m+1

for j=1:m+1

if (abs(i-j)==1)

A(i,j)=-1*r;

end

if (i==j)

A(i,j)=1+2*r;

end

A(1,2)=-2*r;

A(m+1,m)=-2*r;

for i=1:n

for j=1:m+1

F(j)=u0(j)+tao*sin(tao*i);

end

u1=A\F;

for j=1:m+1

u0(j)=u1(j);

end

for j=1:m+1

u(j)=u1(j);

end

最后在命令窗口分别运行parabolicFD(40,1600)，parabolicFD(80,3200)，

matlab求抛物问题,高手进，用matlab求解一维抛物方程的初边值问题 - 数学 - 小木虫 - 学术科研互动社区...相关推荐

matlab文件序号超出511,求教一段matlab的代码 - 数学 - 小木虫 - 学术科研互动社区...
不知楼主说的没法求出优化值,是什么意思.是报错还是怎么. 我也是刚接触优化,觉得楼主的错误可能是在:主函数里 A=[-1;1];的括号中间应该用逗号,而不是分号.即改为A=[-1,1]; 我对你的程序 ...
学习 lt MATLAB gt 心得,lt;lt;MATLAB可视化大学物理学gt;gt;使大学物理更具体，更有趣。 - 物理 - 小木虫 - 学术科研互动社区...
现代大学生除了会运用高等数学解决大学物理中的问题之外,还应该学会应用计算机语言--特别是MATLAB--解决问题. 通过MATLAB可帮助学生学习物理知识.例如,麦克斯韦速率分布率看起来很复杂,但是, ...
matlab 比较两个函数,Matlab同时拟合两个函数 - 数学 - 小木虫 - 学术科研互动社区...
大家好!我想对两组数据(数据和图见下)同时进行拟合,得到一组参数.不知哪位高手能编写出程序,最后生成:原始数据,拟合数据,参数值,每个参数的误差以及总方差. 当x较小时,y变化较快,需要对x求一下对数 ...
matlab换热器优化设计,matlab用于求解一个简单的优化设计问题 - 数学 - 小木虫 - 学术科研互动社区...
最近在学习工程优化设计与MATLAB实现,在模拟课本中的例子的时候,程序提示不正确. Error using ==> mrdivide Matrix dimensions must agree. ...
matlab级数求和阶乘,(阶乘)数列求和 - 数学 - 小木虫 - 学术科研互动社区
翻看了一下坛子里的相关内容,觉得应该是发在这里最合适了.要是这里不合适,请大家告知我一下. 我在对一化学过程数学建模时,得到一数学关系,是一数列,其中各项表达如下: [latex]{x_1} = \f ...
matlab 打开软件报错,matlab程序错误，提示如下【弄了半夜还是不行】 - 数学 - 小木虫 - 学术科研互动社区...
function [x,y,Y]=levydistrintheorytest(t0,lambda,T,a,b); % 计算理论的二维levy分布的图象画出来的3维曲面图 % t0 :lev ...
用MATLAB做结构灵敏度分析怎么做,如何做灵敏度分析 - 数学 - 小木虫 - 学术科研互动社区...
我有以下数据年份平均气温最高气温最低气温相对湿度日照时数实际风速日均ET 1953 ...
hermite矩阵matlab,Hermite矩阵的特征值计算问题 - 数学 - 小木虫 - 学术科研互动社区...
楼主,首先要知道Hermite矩阵就是实对称矩阵的复数版本.我更愿意用:A^H=A来定义Hermite矩阵,这里A是任何一个n阶的复数矩阵,H表示共轭转置.你的这个问题有三个步骤,但是我不清楚你的意图 ...
matlab中lower,【求助】matlab,这个错误究竟是什么？ - 数学 - 小木虫 - 学术科研互动社区...
谢谢以上两位,程序改动了一下,以前的问题暂时没有出现,但是出现了下边的问题: 主程序中是这样写的,调用curvefunzscanAbs作为拟合方程,其中NonAbs是拟合参数,transIntUp,p ...
matlab如何进行四维拟合,如何用matlab对四维数据进行非线性拟合？ - 数学 - 小木虫 - 学术科研互动社区...
这个很简单啊,用lsqcurvefit或者lsqnonlin都可以.比如: function ydata = myfun(a,xdata) ydata=xdata(:,1).*xdata(:,2).* ...

matlab求抛物问题,高手进，用matlab求解一维抛物方程的初边值问题 - 数学 - 小木虫 - 学术科研互动社区...

matlab求抛物问题,高手进，用matlab求解一维抛物方程的初边值问题 - 数学 - 小木虫 - 学术科研互动社区...相关推荐

最新文章

热门文章

matlab求抛物问题,高手进，用matlab求解一维抛物方程的初边值问题 - 数学 - 小木虫 - 学术 科研 互动社区...

matlab求抛物问题,高手进，用matlab求解一维抛物方程的初边值问题 - 数学 - 小木虫 - 学术 科研 互动社区...相关推荐

最新文章

热门文章

matlab求抛物问题,高手进，用matlab求解一维抛物方程的初边值问题 - 数学 - 小木虫 - 学术科研互动社区...

matlab求抛物问题,高手进，用matlab求解一维抛物方程的初边值问题 - 数学 - 小木虫 - 学术科研互动社区...相关推荐