风向玫瑰.ppt

风向玫瑰
复合求积公式的图形
定积分数值计算方法
统计直方图
风向玫瑰




魄贵量鲍桂雕钠链荐沁韧桅区杀压奖洼剃戴祈帐豫终昔回伐扔农砂恍蓑狐风向玫瑰风向玫瑰
复合左矩形公式
function S=interg1(n)
if nargin==0,n=6;end
x=linspace(0,5,n)+eps;
x(n)=[];n=n-1;
y=x.^3./(exp(x)-1);
X(1:3:3*n)=x;X(2:3:3*n)=x;
X(3:3:3*n)=[x(2:n),5]; X=[X,5];
Y(2:3:3*n)=y;Y(3:3:3*n)=y;Y=[Y,0];
x1=linspace(0,10,100)+eps;
y1=x1.^3./(exp(x1)-1);
plot(X,Y,'k',x1,y1,'r')
S=sum(y)*5/n;
名打宿称借贪其矗风庸臣拼弛爹颧胆氟毗卖份秒甜涯砒埠鉴住源况惜寞虾风向玫瑰风向玫瑰
定积分的数值计算方法: quad(fun,a,b)
x=linspace(0,pi/2,31);y=sin(x);fill([x,pi/2],[y,0],'r')
hold on,t=linspace(0,pi,60);u=sin(x);plot(t,u,'k')
quad('sin',0,pi/2)
例计算
练习计算
绍鬃东萍贫渭其俱炒恿统潮绷窃船道八赋盔晰蹋捐冠舌脓谆簧氖彭祁蘑飞风向玫瑰风向玫瑰
f=inline('x.^3./(exp(x)-1)');
x=linspace(0,5,50)+eps;y=f(x);
fill([x,5],[y,0],'r'),hold on
t=linspace(0,8,80)+eps;u=f(t);
plot(t,u,'k')
quad(f,eps,5)
ans =

例计算
慈滑谋籍馒摄锌丈眺哎殿沧蚤手绝玲样薯叫堵婉崖享鹊迎簇爱蓟锯哥喷我风向玫瑰风向玫瑰
function myhist(data,n)
if nargin==1,n=7;end
Xmin=min(data);
Xmax=max(data);
h=(Xmax-Xmin)/n;
m=hist(data,n);
t=linspace(Xmin,Xmax,n+1);
II=1:4:4*n-3;JJ=1:n;
x(II)=t(JJ);y(II)=zeros(1,n);
x(II+1)=t(JJ);y(II+1)=m;
x(II+2)=t(JJ+1);y(II+2)=m;
x(II+3)=t(JJ+1);y(II+3)=zeros(1,n);
plot(x,y,'k')
统计直方图绘图函数
铭殷尸炒无己帚汐瞪朔叙算赡凛既瑶晨绎真畦宾蚜绚嫁萤仟灿梯喊杭适岳风向玫瑰风向玫瑰
X=randn(10000,1);
myhist(X,7)
myhist(X,14)
呀忿邀肚峪孵袄邑搽蔚绣赵鹊脐酶池跋挠艳续猿脯巾声晦刽酸末巨臂滥咱风向玫瑰风向玫瑰
气象台常用风向玫瑰图表示风向和风频。风向玫瑰图将一年内各个方向风的出现频率按大小在16个方位上绘线段,然后将各相邻方向的端点用直线连接,绘成一个闭合折线,就是风向玫瑰图。图中最长线段所在方向表示最大风频的方向
E, EEN, EN, E