matlab计算国土面积程序

MATLAB--计算国土面积 数值分析 课程设计 历时两周 完美通过!!

国土面积的计算

摘要：数学建模方法是处理科学理论的一种经典方法，也是解决各类实际问题的常用方法。本文采用插值、复化梯形公式的方法，并利用数学软件MATLAB对国土面积进行计算，计算结果与实际记录基本吻合，最终将误差控制在1%内。 关键词：建模 分段 插值 复化梯形公式(trapz)方法 MATLAB

1．问题分析

在生产与科学试验中，反映自然规律的函数关系，往往是通过实验、观测得到的。从题目我们知道本题的本质是就是求不规则平面图形的面积。题目给出地图测量数据，我们就可以使用MATLAB软件的plot函数画出该不规则平面图形的模拟图。题中已经把图形分成上下边疆线,所以可以运用插值与复化梯形公式分别求出上下边疆的面积，最后求出图形面积。

2．算法设计

2．1作图

因为表中数据点较多，并且横轴有重复的点，无法直接进行插值，所以需要进行分段。分别把上下边疆线按数字连续截取，遇到重复的值，或小的值进行分段。

记为xn,其中n=1,2,3,……

以下边疆线为例：

把下边疆线分为6段，分别是：

x1=[17 18 20 31 41 58 66 72 ];

x2=[57 60 71 104 1

可以大致了解我们的国家

中国历代国土面积

中国历代国土面积

秦朝的面积为：960万-608.61万+3.30万平方公里=354.69万平方公里。

西汉的面积为：354.69万+311.68万=666.37万平方公里

东汉的面积为：666.37万+33.53万-45.28万=654.62万平方公里

西晋的面积为：654.62万-39.12万=615.5万平方公里

唐（唐高宗时期）的面积：470.57万+809.47万-28.85万=1251.19万平方公里

唐（唐玄宗时期）的面积：1251.19万+184.13万-545.73万=889.59万平方公里

吐蕃：453.88万平方公里

辽：448.54万平方公里

北宋的面积：889.59万-606.03万=283.56万平方公里

元朝的面积：960万+1207.82万-45.08万=2122.74万平方公里

明朝的面积：960万+449.19万-175.81万=1233.38万平方公里

清朝的面积：960万+324.65万=1284.65万平方公里

中国各朝国土的平方公里数如下：

1、夏（210万，包括长江以北的湖北，河南，安徽，山东，河北，山西，京津和江苏，陕西的一部分）

2，商（320万，除上述土地外，还包括陕西江苏的剩余土地

可以大致了解我们的国家

中国历代国土面积

中国历代国土面积

秦朝的面积为：960万-608.61万+3.30万平方公里=354.69万平方公里。

西汉的面积为：354.69万+311.68万=666.37万平方公里

东汉的面积为：666.37万+33.53万-45.28万=654.62万平方公里

西晋的面积为：654.62万-39.12万=615.5万平方公里

唐（唐高宗时期）的面积：470.57万+809.47万-28.85万=1251.19万平方公里

唐（唐玄宗时期）的面积：1251.19万+184.13万-545.73万=889.59万平方公里

吐蕃：453.88万平方公里

辽：448.54万平方公里

北宋的面积：889.59万-606.03万=283.56万平方公里

元朝的面积：960万+1207.82万-45.08万=2122.74万平方公里

明朝的面积：960万+449.19万-175.81万=1233.38万平方公里

清朝的面积：960万+324.65万=1284.65万平方公里

中国各朝国土的平方公里数如下：

1、夏（210万，包括长江以北的湖北，河南，安徽，山东，河北，山西，京津和江苏，陕西的一部分）

2，商（320万，除上述土地外，还包括陕西江苏的剩余土地

世界各国国土面积排名以及语言&中国各省面积排名

一,世界各国国土面积排名(万平方公里)

面积排名国家国土面积

1 俄罗斯1,707.5

2 加拿大997.1

3 中国960.1

4 美国936.4

5 巴西854.7

6 澳大利亚774.1

7 印度328.8

8 阿根廷278.0

9 哈萨克斯坦271.7

10 苏丹250.6

11 阿尔及利亚238.2

12 刚果(金) 234.5

13 沙特阿拉伯215.0

14 墨西哥195.8

15 印度尼西亚190.5

16 利比亚176.0

17 伊朗163.3

18 蒙古156.7

19 秘鲁128.5

20 乍得128.4

21 尼日尔126.7

22 安哥拉124.7

23 马里124.0

24 南非122.1

25 哥伦比亚113.9

26 埃塞俄比亚110.4

27 玻利维亚109.9

28 毛里塔尼亚102.6

29 埃及100.1

30 坦桑尼亚94.5

31 尼日利亚92.4

32 委内瑞拉91.2

33 纳米比亚82.4

34 莫桑比克80.2

35 巴基斯坦79.6

36 土耳其77.5

37 智利75.7

38 赞比亚75.3

39 缅甸67.7

40 阿富汗65.2

41 索马里63.8

42 中非62.3

43 乌克兰60.4

44 马达加斯加58.7

45 博

北航强度与振动大作业

%叶片强度，离心拉伸应力，气动弯矩，离心力弯矩，合成弯矩，ABC三点的弯曲应力和总应力

midu=8.2*10^3; %密度

zhuansu=4700*pi/30; %转速rad/s

w=midu*zhuansu*zhuansu;

X=[0.0053;0.0041;0.0041;0.0040;0.0024;0.0012];

Y=[-0.0041;-0.0038;-0.0030;-0.0019;-0.0011;-0.0002];

Z=[0.628;0.591;0.56;0.53;0.494;0.458];

A=[0.00018;0.000232;0.000312;0.00041;0.000548;0.000705];

XX=X; YY=Y; ZZ=Z; AA=A;

X(1)=[]; Y(1)=[]; Z(1)=[]; A(1)=[]; %后五行

XX(6)=[]; YY(6)=[]; ZZ(6)=[]; AA(6)=[]; %前五行

Xm=(XX+X)/2; %中间插值

Ym=(YY+Y)/2; %中间插值

Zm=(ZZ+Z)/2; %中间插值

Am=(AA+A)/2; %中间插值

dZ=ZZ-Z; %逐差

%离心拉伸应力

dPc=w.*A

中国历代疆域地图&中国历代国土面积一览表

夏朝:

大约在公元前2070年，禹的儿子启废除禅让制度，在河南禹县的钧台称王，建国号夏。夏朝的领土在陕西、山西、河北、山东、河南之间，面积大约有50万平方公里。夏朝的末代君王叫桀，扰民太甚，天下人反。

商朝:

大约在公元前1600年，夏所属部落首领成汤灭夏，定都河南商丘，建国号商。商君武丁扩土四方，占领湖北、安徽、江苏、内蒙古等地区,面积扩大到100万平方公里。公元前12世纪，商王盘庚迁都

到殷（今河南安阳），所以商朝又叫殷朝或殷商。商朝的末代君王是纣王。纣王扰民太甚，天下人反。

西周：

大约在公元前1046年，商所属部落首领姬发灭商，定都陕西西安，建国号周，史称西周。姬发大封诸侯。周朝初期的面积与商朝差不多，大约有100平方公里，拥有河南、河北、山西、陕西、甘肃、山东等地区。后来又扩充，东北边到达黑龙江地区，东边临海，东南方最远到达江苏无锡地区，南边到了安徽巢湖地区，西边到了甘肃东部，北部占有内蒙古的部分地区，面积大约有250万平方公里。

春秋战国：

公元前770年，西周遗臣姬宜臼称王，定都河南洛阳，建国号周，史称东周。东周各代君王无能，所属各诸侯国纷纷自强，史称“春秋战

电工专业必看

%简单潮流计算的牛顿拉夫逊程序，相关的原始数据数据数据输入格式如下： %B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写 %对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点 %编号，将变压器的串联阻抗置于低压侧处理。

%第三列为支路的串列阻抗参数。

%第四列为支路的对地导纳参数。

%第五列为含变压器支路的变压器的变比

%第六列为变压器是否含有变压器的参数，其中“1”为含有变压器，

%“0”为不含有变压器。

%B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。

%X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地 %参数。

n=input('请输入节点数:n=');

n1=input('请输入支路数:n1=');

isb=input('请输入平衡节点号:isb=');

pr=input('请输入误差精度:pr=');

B1=input('请输入支路参数:B1=');

B2=input('请输入节点参数:B2=');

X=input('节点号和对地参数:X=');

Y=zeros(n

%% 汽车动力性计算（自己编的动力性计算程序，供大家计算动力性时参考，具体参数大家根据所给程序对应输入，并对坐标轴数值按需要进行修改） clc; clear; close all;

%%根据所给发动机数据拟合外特性曲线（发动机数据按照你所得到的数据进行输入）

n_test=[500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 2100 2200];

T_test=[975 1108 1298 1496 1546 1620 1670 1785 1974 1974 1970 1889 1829 1748 1669 1700 1524 1105]; figure(1)

plot(n_test,T_test,'g'); hold on grid on

%p=polyfit(n_test,T_test,7); p=polyfit(n_test,T_test,2); n=[450:1:2200]; Ttq=polyval(p,n); plot(n,Ttq,'k');

xlabel('发动机转速n(r/min)'); ylabel('

matlab编写的Lyapunov指数计算程序

已有 2406 次阅读 2009-12-29 08:37 |个人分类:其它|系统分类:科普集锦|关键词:李雅普诺夫指数

一、计算连续方程Lyapunov指数的程序 其中给出了两个例子：

计算Rossler吸引子的Lyapunov指数

计算超混沌Rossler吸引子的Lyapunov指数

http://www.pudn.com/downloads39/sourcecode/math/detail132231.html

二、recnstitution重构相空间，在非线性系统分析中有重要的作用 function [Texp,Lexp]=lyapunov(n,tstart,stept,tend,ystart,ioutp); global DS; global P;

global calculation_progress first_call; global driver_window;

global TRJ_bufer Time_bufer bufer_i; %

% Lyapunov exponent calcullation for ODE-system. %

% The al

disp(' ')

disp(' srs.m ver 2.0 July 3, 2006')

disp(' by Tom Irvine Email: tomirvine@aol.com') disp(' ')

disp(' This program calculates the shock response spectrum') disp(' of an acceleration time history, which is pre-loaded into Matlab.') disp(' The time history must have two columns: time(sec) & acceleration') disp(' ') %

clear t; clear y; clear yy; clear n; clear fn; clear a1; clear a2 clear b1; clear b2; clear jnum; clear THM; clear resp; clear x_pos; clear x_neg; %

iunit=input(' Enter acceleration unit: 1= G