函数文件1: function b=F(x0,h,u,tau) b(,)=x0()-u(); b(,)=x0()-u()+*h*1e8*cos(tau)*x0(); 函数文件2: function g=Jacobian(x0,h,tau) g(,)=; g(,)=; g(,)=; g(,)=+*h*1e8*cos(tau); 函数文件3: function x=Euler(h,x0,u,tau) % x0 表示迭代初值 % u 表示上一节点数值 tol=1e-; x1=x0-Jacobian(x…

% forward Euler demo % take two steps in the solution of % dy/dt = y, y(0) = 1 % exact solution is y(t) = exp(t) clear all close all % the exact solution t = 0:0.01:1; yexact = exp(t); % initial conditions y0 = 1 h = 0.5 tt = [0]; yy = [y0]; % one st…

% matlab script to test efficiency of % Euler's method, classical Runge-Kutta, and ode45 % on Arenstorf orbit problem close all clear all % these are variables we would like the right-hand function to "see" % without actually passing them as arg…

一.显示Euler 函数文件:Euler.m function f=Euler(h,Y) f(1,1)=Y(1)+h*(0.01-(1+(Y(1)+1000)*(Y(1)+1))*(0.01+Y(1)+Y(2))); f(2,1)=Y(2)+h*(0.01-(1+Y(2)^2)*(0.01+Y(1)+Y(2))); 脚本文件: tic; clear clc %% %显示Euler方法求刚性微分方程,要求用Richardson外推法估计近似误差从而控制步长 y(1:2,1)=[0;0];%初值 e…

函数文件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 每次迭代的…

MATLAB地图工具箱学习总结(三)地图工具箱的基本知识 今天想要介绍的是一些比较基础的函数.了解了这些函数,地图投影的基本概念才能真正明白.而要想继续研究MATLAB中有关地图投影的函数,尤其是未来我要提到的投影文件源代码,知晓这些函数的功能必不可少.本篇文章将会罗列三个案例,并在后面一一进行讲解. 1                    作业案例:地图投影作业1 这次的案例从作业1开始.作业1是要求计算出地球椭球体的一些基本参数,包括子午圈曲率半径.卯酉圈曲率半径.平均曲率半径和纬圈半径等…

举例:分别用欧拉法和龙哥库塔法求解下面的微分方程 我们知道的欧拉法(Euler)"思想是用先前的差商近似代替倒数",直白一些的编程说法即:f(i+1)=f(i)+h*f(x,y)其中h是设定的迭代步长,若精度要求不高,一般可取0.01.在定义区间内迭代求解即可.龙哥库塔法一般用于高精度的求解,即高阶精度的改进欧拉法,常用的是四阶龙哥库塔,编程语言如下:y(i+1)=y(i)+h*(k1+2*K2+2*k3+k4)/6;k1=f(xi,yi)k2=f(xi+h/2,yi+h*k1/2);…

生成简单多边形后,有时还需要对多边形各顶点的凹凸性做判断. 先计算待处理点与相邻点的两个向量,再计算两向量的叉乘,根据求得结果的正负可以判断凹凸性. 结果为负则为凹顶点,为正则为凸顶点. 凹顶点用o表示,凸顶点用*表示. 结果如下: matlab代码如下: clear all;close all;clc; n=; p=rand(n,); p=createSimplyPoly(p); %创建简单多边形 hold on; :n %处理第一个点 v1=p(n,:)-p(,:); %当前点到前一点向量…

题目链接:poj2284 That Nice Euler Circuit 欧拉公式:如果G是一个阶为n,边数为m且含有r个区域的连通平面图,则有恒等式:n-m+r=2. 欧拉公式的推广: 对于具有k(k≥2)个连通分支的平面图G,有:n-m+r=k+1. 题意:给出连通平面图的各顶点,求这个欧拉回路将平面分成多少区域. 题解:根据平面图的欧拉定理“n-m+r=2”来求解区域数r. 顶点个数n:两两线段求交点,每个交点都是图中的顶点. 边数m:在求交点时判断每个交点落在第几条边上,如果一个交点落在…

K-means是一种经典的聚类算法,是十大经典数据挖掘算法之一.K-means算法的基本思想是:以空间中k个点为中心进行聚类,对最靠近他们的对象归类.通过迭代的方法,逐次更新各聚类中心的值,直至得到最好的聚类结果. 假设要把样本集分为c个类别,算法描述如下: (1)适当选择c个类的初始中心: (2)在第k次迭代中,对任意一个样本,求其到c个中心的距离,将该样本归到距离最短的中心所在的类: (3)利用均值等方法更新该类的中心值: (4)对于所有的c个聚类中心,如果利用(2)(3)的迭代法更新后,值…