加载中…MATLAB小技巧之积分求面积
(2013-04-16 21:25:45) | 分类： Matlab |
clc;
clear;
% 假设有一组数据,求取拟合曲线
x = [1 2 3 4 5 6 7 8]
y=[1.2 1.5 1.3 2.4 2.3 2.6 3.2 3.4];
a=polyfit(x,y,2); %首先拟合曲线函数，选择合理的拟合曲线参数
yi=polyval(a,x);
plot(x,y,'o',x,yi,'r');
% trapz
梯形法沿列方向求函数Y关于自变量X的积分(向量形式，数值方法)
mm = max(x);
nn = 100; % 网格划分密度
m = 0:mm/nn:mm % 拟合曲线的x轴范围
for i = 1:length(m)
end
S = trapz(m,n); % 求解面积
hold on;
title(['the area is ',num2str(S)]);
喜欢
0
赠金笔