galerkin有限元法matlab实现,PDE的Galerkin和有限元的MATLAB程序

%Galerkin方法 clear all clc syms x;

A=zeros(4,4); for i=1:4

for j=1:4

phy1=x^i; phy2=x^j;

dphy1=diff(phy1,1); dphy2=diff(phy2,1); phy=pi^2*phy1*phy2; dphy=dphy1*dphy2;

A(i,j)=int(phy+dphy,x,0,1); end end D=[]; for k=1:4

f1=2*pi^2*sin(pi*x)+pi^3*x; f2=x^k; f=f1*f2;

D(k)=int(f,x,0,1); end D=D'; C=A\\D; C=C';

X=linspace(0,1,6); F=0; for i=1:4

F=F+C(i)*x^i; end for j=1:6

Y(j)=subs(F,X(j)); end

Y=Y-pi.*X; Y1=sin(pi.*X);

err=norm(abs(Y-Y1)); disp('数值解') disp(Y)

disp('整体误差') disp(err)

%%%%%%%%%追赶法 function x=chase(a,b,c,d) n=length(b); u(1)=c(1)/b(1); q(1)=d(1)/b(1); for i=2:n-1

h(i)=b(i)-u(i-1)*a(i-1); u(i)=c(i)/h(i);

q(i)=(d(i)-q(i-1)*a(i-1))/h(i); end

q(n)=(d(n)-q(n-1)*a(n-1))/(b(n)-u(n-1)*a(n-1)); x(n)=q(n); for i=n-1:-1:1

x(i)=q(i)-u(i)*x(i+1); end end

%有限元法 clear all clc m=5;

xspan=[0 1];

a=xspan(1);b=xspan(2);h=(b-a)/m; x=[a:h:b];

f=zeros(m-1,1);

F1=@(t)(-1/h+(pi^2)*t.*(1-t)*h); c1=quadl(F1,0,1);

F2=@(t)(2/h+(pi^2)*(t.^2+(1-t).^2)*h); c2=quadl(F2,0,1);

A=diag(c2*ones(m-1,1))+diag(c1*ones(m-2,1),1)+... diag(c1*ones(m-2,1),-1); for i=2:m

F=@(t)h*t.*(2*pi^2*sin(pi*(x(i-1)+h*t)))+... h*(1-t).*(2*pi^2*sin(pi*(x(i)+h*t))); f(i-1)=quadl(F,0,1); end

v=chase(c1*ones(m-2,1),c2*ones(m-1,1),c1*ones(m-2,1),f); u=[0,v,0]

galerkin有限元法matlab实现,PDE的Galerkin和有限元的MATLAB程序相关推荐

matlab计算应力位移,2012年-有限元作业-matlab编程实现有限元求解简单结构位移及应力.doc...
<2012年-有限元作业-matlab编程实现有限元求解简单结构位移及应力.doc>由会员分享,可在线阅读,更多相关<2012年-有限元作业-matlab编程实现有限元求解简单结构位 ...
二维有限元方程matlab,有限元法求解二维Poisson方程的MATLAB实现
有限元法求解二维 Poisson 方程的 MATLAB 实现陈莲a ,郭元辉b ,邹叶童a ( 西华师范大学 a．数学与信息学院; b．教育信息技术中心,四川南充 6437009) 摘要: ...
matlab的PDE工具箱的简单使用
matlab的PDE工具箱的简单使用_LSEC小陆的博客-CSDN博客_matlab pde pdetool工具箱的使用浅析_DavidEnterpriseStar的博客-CSDN博客_pdetool
matlab模拟线圈电磁场,利用MATLAB的PDE工具箱对电场和磁场进行模拟
中学物理Vol. 32 No.巧疑的态度,认为弊大于利,笔者认为这些想法不无道理,只是大家的关注点应该在于如何更好发挥它的优势,而避免它的不足．在实践的过程中笔者有以下心得和体会,与大家分享: 3, ...
matlab求偏微分方法解析解,偏微分数值解(2,MATLAB求解方法)学案.ppt
这部分主要讨论如何用MATLAB实现对偏微分方程的数值仿真求解．MATLAB的偏微分方程工具箱(PDE Toolbox)的出现,为偏微分方程的求解以及定性研究提供了捷径．主要步骤为: 2.1 用偏微分 ...
李萍matlab实验报告,李萍, 张磊, 王垚廷. 基于Matlab的偏微分方程数值计算[J]. 齐鲁工业大学学报, 2017, 31(4): 39-43....
ABSTRACT: 半导体是近代发现的一种新型材料,其导电性能介于导体与绝缘体之间,并且导电性也会随温度变化,是制作电子器件的重要材料,因此半导体器件问题的数值模拟也越来越重要.本文主要研究用Matl ...
MATLAB计算杨氏模量,四阶弹性模量Cijkl如何在matlab里表示啊？ - 计算模拟 - 小木虫 - 学术科研互动社区...
matlab 四元数运算计算包就可以了吧 Matlab 四元数操作函数 2012-06-03 21:02:55| 分类: MATLAB&Mathemati | 标签:四元数 quater ...
在matlab中有几种获得帮助的途径,matlab经典习题及解答
精品文档 . 1欢迎下载第1章 MATLAB 概论 1.1 与其他计算机语言相比较,MATLAB 语言突出的特点是什么? MATLAB 具有功能强大.使用方便.输入简捷.库函数丰富.开放性强等特点. ...
【数字信号处理】卷积编程实现 ( Matlab 卷积和多项式乘法 conv 函数 | 使用 matlab 代码求卷积并绘图 )
文章目录一.Matlab 卷积和多项式乘法 conv 函数二.使用 matlab 代码求卷积并绘图一.Matlab 卷积和多项式乘法 conv 函数 Matlab 文档地址 : https:// ...

galerkin有限元法matlab实现,PDE的Galerkin和有限元的MATLAB程序

galerkin有限元法matlab实现,PDE的Galerkin和有限元的MATLAB程序相关推荐

最新文章

热门文章