有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,一维有限元法解常微分方程相关推荐

  1. 中点法 matlab,中点法解常微分方程(组)

    &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp 预备知识 常微分方程(组)的数值解 我们先来尝试用欧拉法解一 ...

  2. 一维热传导方程 matlab隐式解,一维热传导偏微分方程的数值解的matlab程序问题出在哪儿?...

    我现在编写了一个求解一维热传导的偏微分方程,调程序都调了好多天了 不知道问题在哪儿,求各位高手帮忙看一下好么? 我在此表示万分感谢 需求解的方程看图片,我的程序如下 %---------------- ...

  3. 有限差分法的一维扩散MATLAB,一维扩散方程的有限差分法matlab

    用matlab编程实现一维扩散方程的有限差分法 1 一维扩散方程的有限差分法 --计算物理实验作业七 陈万 物理学2013级 130******** ● 题目: 编程求解一维扩散方程的解 ⎪⎪⎪⎪⎩ ...

  4. galerkin有限元法matlab实现,有限元法求解二维Poisson方程的MATLAB实现

    有限元法求解二维Poisson方程的MATLAB实现 陈莲a,郭元辉b,邹叶童a [摘要]文章讨论了圆形区域上的三角形单元剖分.有限元空间,通过变分形式离散得到有限元方程. 用MATLAB编程求得数值 ...

  5. matlab寻访一维数组,MATLAB 一维数组的创建和寻访

    3.2 一维数组的创建和寻访 3.2.2 一维数组的子数组寻访和赋值 [ * 例 3.2.2 -1 ]子数组的寻访( Address ). rand('state',0) %把均匀分布伪随机发生器置为 ...

  6. [Matlab科学计算] 四阶Runge-Kutta法解常微分方程

    四阶Runge-Kutta法格式的详细推导请查找相关数值分析书籍,这里直接给出四阶Runge-Kutta法的经典格式和Matlab代码 Matlab代码如下:自行修改常微分方程即可 %% 四阶Rung ...

  7. 【Matlab】一、解常微分方程ODE

    文章目录 求解常微分方程 ODE (1)求解解析解 (2)求解数值解 求解常微分方程 ODE ​ 在matlab中,我们可以求解常微分方程的解析解,和数值解,一般使用dsolve来求解常微分方程的解析 ...

  8. matlab求微分方程的初值问题,Matlab 解常微分方程的初值问题.doc

    Matlab 解常微分方程的初值问题 题目:Matlab 解常微分方程的初值问题 设计目的: 1.熟练掌握Matlab的基本编程方法,及其编程风格. 2.熟练掌握Matlab常用函数的使用. 3.与本 ...

  9. matlab:一维数组

    1:如何生成一维数组 方法一:用:生成等差向量 >> a=2:3:20 a =2 5 8 11 14 17 20 方法二:linspace() 函数 >> a=linspace ...

最新文章

  1. JS操作JSON总结
  2. 机房布线的最高境界……
  3. Centos7 搭建FastDFS文件管理系统及简单的迁移方法
  4. Docker批量操作容器
  5. gevent.joinall()开启协程
  6. Faas 和无服务器架构
  7. Lesson 1#03-Python安装与Hello Python World
  8. MiniGUI+qvfb+Eclipse完成嵌入式系统图形开发(开发环境的搭建)
  9. 政府采购就是支持民族软件吗?
  10. 小爱同学100个奇葩回复_小米小爱音箱Pro开箱评测,看到的不仅是全面升级更是小米loT的高速发展...
  11. 启动mysql失败:net start mysql服务名无效
  12. futuretask java 并发请求_Java并发机制(9)--Callable、Future、FutureTask的使用
  13. 用户、组和权限命令练习
  14. Java的常用集合框架
  15. ADAMS2016启动证书错误解决
  16. 计算机分屏显示知识,如何使用Windows的分屏显示功能
  17. backtracking line Serach
  18. 抓取lol全英雄图(不含皮肤)
  19. 用C语言求并集和交集
  20. POJ 1436 Horizontally Visible Segments(线段树区间修改)

热门文章

  1. css类选择器优先级 大于 标签选择器
  2. SAP C4C,CRM和S4HANA的Saved Query使用介绍
  3. Angular component的职责
  4. javax.ws.rs.Path注解@Path的工作原理解析
  5. 如何使用SAP UI5 web Component的React框架的柱状图和折线图
  6. Event handling in Angular
  7. Fiori GM4取metadata 从http自动重定向到https原因
  8. Message popover
  9. SalesOrder Text determination
  10. How to find all event handler list for a given event