p47.(实习题-李荣华)用线性元求下列边值问题的数值解 tic; % this method is transform from Galerkin method %also call it as finit method %is used for solving two point BVP which is the first and second term. % %MATLAB 7.0 clear; clc; N=; h=/N; X=:h:; f=inline('(0.5*pi^2)*sin…

MATLAB常微分方程的数值解法 作者:凯鲁嘎吉 - 博客园http://www.cnblogs.com/kailugaji/ 一.实验目的 科学技术中常常要求解常微分方程的定解问题,所谓数值解法就是求未知函数在一系列离散点处的近似值. 二.实验原理 三.实验程序 1. 尤拉公式程序 四.实验内容 选一可求解的常微分方程的定解问题,分别用以上1, 4两种方法求出未知函数在 节点处的近似值,并对所求结果与分析解的(数值或图形)结果进行比较. 五.解答 1. 程序 求解初值问题 取n=10 源程序:…

原文出处http://wenku.baidu.com/view/df412e115f0e7cd184253653.html 因为不太喜欢百度文库的格式,所以写到个人博客里面方便使用 <iframe width='738' height='523' class='preview-iframe' scrolling='no' frameborder='0' src='http://download.csdn.net/source/preview/8126015/065f7104461e354c252…

罗列一下以前自己学习C语言与MATLAB混编的笔记,顺便复习一遍. <C语言与MATLAB接口 编程与实例 李传军编着>(未看完,目前看到P106) 目录P4-8 ****************************************************** C-MEX函数:从MATLAB调用C语言或Fortran语言编写的函数.   MEX文件:是一种按一定格式,使用C语言或者FORTRAN语言编写的,由MATLAB解释器自动调用并执行的动态链接函数.在Windows系统中,M…

欧拉法的来源 在数学和计算机科学中,欧拉方法(Euler method)命名自它的发明者莱昂哈德·欧拉,是一种一阶数值方法,用以对给定初值的常微分方程(即初值问题)求解.它是一种解决常微分方程数值积分的最基本的一类显型方法(Explicit method). [编辑] 什么是欧拉法 欧拉法是以流体质点流经流场中各空间点的运动即以流场作为描述对象研究流动的方法.——流场法 它不直接追究质点的运动过程,而是以充满运动液体质点的空间——流场为对象.研究各时刻质点在流场中的变化规律.将个别流体质点运动过…

机器人运动学逆解的问题经常出现在动画仿真和工业机器人的轨迹规划中:We want to know how the upper joints of the hierarchy would rotate if we want the end effector to reach some goal. IK Solutions: Analytical solutions are desirable because of their speed and exactness of solution. For…

什么是对象 JavaScript 中的所有事物都是对象,如:字符串.数值.数组.函数等,每个对象带有属性和方法. 对象的属性:反映该对象某些特定的性质的,如:字符串的长度.图像的长宽等: 对象的方法:能够在对象上执行的动作.例如,表单的"提交"(Submit),时间的"获取"(getYear)等: JavaScript 提供多个内建对象,比如 String.Date.Array 等等,使用对象前先定义,如下使用数组对象: var objectName =new Ar…

函数文件1: function b=F(f,x0,h,N) % b(1,1)=x0(1)-h*x0(2)-u(1); % b(2,1)=x0(2)+h*x0(1)^2-u(2)-h*f; b=zeros(N,1); b(1,1)=4*x0(1)-x0(2); b(2,1)=h^2*x0(1)^2-2*x0(1)+x0(2)-h^2*f(1) for i=3:N b(i,1)=x0(i-2)+h^2*x0(i-1)^2-2*x0(i-1)+x0(i)-h^2*f(i-1); end 函数文件2:…

函数文件1: function b=F(f,x0,u,h) b(1,1)=x0(1)-h*x0(2)-u(1); b(2,1)=x0(2)+h*x0(1)^2-u(2)-h*f; 函数文件2: function g=Jacobian(x0,u,h) g(1,1)=1; g(1,2)=-h; g(2,1)=2*h*x0(1); g(2,2)=1; 函数文件3: function x=newton_Iterative_method(f,u,h) % u:上一节点的数值解或者初值 % x0 每次迭代的…

tic; % this method is transform from Ritz method %is used for solving two point BVP %this code was writen by HU.D.dong in February 11th 2017 %MATLAB 7.0 clear clc N=50; h=1/N; X=0:h:1; f=inline('pi^2/2*sin(pi/2*x)'); %以下是右端向量: for i=2:N fun1=@(x) f(X…