一维有限元法matlab,一维有限元法解常微分方程
有matlab程序
一维问题的有限元法一.算法构思
考虑下面的两点边值问题
''2
(0,1) Fu=-u+u=sin x 1+ =f,x
y(0) 0,y(1) 0
设F是一个微分算子,则
1
Fu=0 Fu vdx=0, v C0,即v(0)=v(1)=0,且v'(x)连续
01
则把问题中的微分方程化为积分方程,得 u'v'+uvdx= fvdx
1
1
令 a u,v u'v'+uvdx, (f,v)= fvdx
1则问题就是求u C10,使得a(u,v)=(f,v), v C0
1对于一般的u C0,其范围很广泛,但样条函数理论给我们提供了解决问题的有力工具。
11
对[0,1]进行等分:
10=x0
n
1
x-xk-1),x [xk-1,xk) h
取 k x = ,k 1,2,...,n 1.
1 - x-x ,x [x,x)
k+1kk+1
h则u(x)= ui i x ,其中ui是u(x)再xi处得函数值.
i=1n-1
a(u,v)=(f,v) a(u,w)=(f,w),w k x ,k 1,2,...,n 1.即 a i, j ui=(f, j),j=1,2,...,n-1.
i=1n-1
, 1 1
AU=b,其中A=
,
1n-1
xk+1
2k
(f, 1) u 1
,b= . ,U= (f, n-1) n-1, n-1 n-1 u
n-1, 1
22h22
k, k = xk-1 x dx+=+
h3h
xk+11h21
k, k+1 = xk k x k+1 x dx-=-h6h
xk+11+ 22
f, k = 1+ xk-1 k x sin x dx=2 2sin xk -sin xk 1 sin xk 1 h
二,程序实现
clear; clc;
x=zeros(1,101); X=0.01:0.01:0.99;
一维有限元法matlab,一维有限元法解常微分方程相关推荐
- 中点法 matlab,中点法解常微分方程(组)
              预备知识 常微分方程(组)的数值解 我们先来尝试用欧拉法解一 ...
- 一维热传导方程 matlab隐式解,一维热传导偏微分方程的数值解的matlab程序问题出在哪儿?...
我现在编写了一个求解一维热传导的偏微分方程,调程序都调了好多天了 不知道问题在哪儿,求各位高手帮忙看一下好么? 我在此表示万分感谢 需求解的方程看图片,我的程序如下 %---------------- ...
- 有限差分法的一维扩散MATLAB,一维扩散方程的有限差分法matlab
用matlab编程实现一维扩散方程的有限差分法 1 一维扩散方程的有限差分法 --计算物理实验作业七 陈万 物理学2013级 130******** ● 题目: 编程求解一维扩散方程的解 ⎪⎪⎪⎪⎩ ...
- galerkin有限元法matlab实现,有限元法求解二维Poisson方程的MATLAB实现
有限元法求解二维Poisson方程的MATLAB实现 陈莲a,郭元辉b,邹叶童a [摘要]文章讨论了圆形区域上的三角形单元剖分.有限元空间,通过变分形式离散得到有限元方程. 用MATLAB编程求得数值 ...
- matlab寻访一维数组,MATLAB 一维数组的创建和寻访
3.2 一维数组的创建和寻访 3.2.2 一维数组的子数组寻访和赋值 [ * 例 3.2.2 -1 ]子数组的寻访( Address ). rand('state',0) %把均匀分布伪随机发生器置为 ...
- [Matlab科学计算] 四阶Runge-Kutta法解常微分方程
四阶Runge-Kutta法格式的详细推导请查找相关数值分析书籍,这里直接给出四阶Runge-Kutta法的经典格式和Matlab代码 Matlab代码如下:自行修改常微分方程即可 %% 四阶Rung ...
- 【Matlab】一、解常微分方程ODE
文章目录 求解常微分方程 ODE (1)求解解析解 (2)求解数值解 求解常微分方程 ODE 在matlab中,我们可以求解常微分方程的解析解,和数值解,一般使用dsolve来求解常微分方程的解析 ...
- matlab求微分方程的初值问题,Matlab 解常微分方程的初值问题.doc
Matlab 解常微分方程的初值问题 题目:Matlab 解常微分方程的初值问题 设计目的: 1.熟练掌握Matlab的基本编程方法,及其编程风格. 2.熟练掌握Matlab常用函数的使用. 3.与本 ...
- matlab:一维数组
1:如何生成一维数组 方法一:用:生成等差向量 >> a=2:3:20 a =2 5 8 11 14 17 20 方法二:linspace() 函数 >> a=linspace ...
最新文章
- JS操作JSON总结
- 机房布线的最高境界……
- Centos7 搭建FastDFS文件管理系统及简单的迁移方法
- Docker批量操作容器
- gevent.joinall()开启协程
- Faas 和无服务器架构
- Lesson 1#03-Python安装与Hello Python World
- MiniGUI+qvfb+Eclipse完成嵌入式系统图形开发(开发环境的搭建)
- 政府采购就是支持民族软件吗?
- 小爱同学100个奇葩回复_小米小爱音箱Pro开箱评测,看到的不仅是全面升级更是小米loT的高速发展...
- 启动mysql失败:net start mysql服务名无效
- futuretask java 并发请求_Java并发机制(9)--Callable、Future、FutureTask的使用
- 用户、组和权限命令练习
- Java的常用集合框架
- ADAMS2016启动证书错误解决
- 计算机分屏显示知识,如何使用Windows的分屏显示功能
- backtracking line Serach
- 抓取lol全英雄图(不含皮肤)
- 用C语言求并集和交集
- POJ 1436 Horizontally Visible Segments(线段树区间修改)
热门文章
- css类选择器优先级 大于 标签选择器
- SAP C4C,CRM和S4HANA的Saved Query使用介绍
- Angular component的职责
- javax.ws.rs.Path注解@Path的工作原理解析
- 如何使用SAP UI5 web Component的React框架的柱状图和折线图
- Event handling in Angular
- Fiori GM4取metadata 从http自动重定向到https原因
- Message popover
- SalesOrder Text determination
- How to find all event handler list for a given event