matlab实现数值积分【一】（trapz函数）

总述

数值积分问题是传统数值分析课程中的重要内容。如果被积函数的数学表达式未知，则需要由实测数据通过梯形算法求出积分的近似值。本文将介绍被积函数的数学表达式未知时数值积分问题的求解方法，即已知数据点求积分。

函数调用格式

S = trapz(x, y);

应用举例

例1：梯形法求积分

x = [0:pi/30:pi]';
y = [sin(x) cos(x) sin(x/2)];
S = trapz(x,y)

结果为：S = [1.9982 0.0000 1.9995]

由于选择的步距较大，为 h＝π／30=0.1h＝\pi／30 = 0.1h＝π／30=0.1 , 故得出的结果有较大的误差。其实可以将积分问题与样条插值技术相结合，给出一个能精确计算积分的MATLAB函数。（待补充）

例2：不同步长对积分结果的影响

题目： 用定步长法求解积分∫03π2cos15xdx\int_0^{\frac{3\pi}{2}}{cos15xdx}∫023πcos15xdx，并讨论不同步长对积分值的影响。

首先，绘制被积函数的图像：

x=[0:0.01:3*pi/2, 3*pi/2];  % //这样赋值能确保 3*pi/2 点被包含在内
y=cos(15*x); plot(x,y)

由图像观察出在求解区域内被积函数有很强的振荡。

对不同的步距 h=0.1,0.01,0.001,0.0001,0.00001,0.000001h= 0.1, 0.01, 0.001, 0.0001, 0.00001, 0.000001h=0.1,0.01,0.001,0.0001,0.00001,0.000001 ,可以用下面的语句求出采用不同步长的积分近似结果。

syms x, A=int(cos(15*x),0,3*pi/2) % //求取理论值为1/15
h0=[0.1,0.01,0.001,0.0001,0.00001,0.000001]; v=[];
for h=h0x = [0:h:3*pi/2, 3*pi/2];  y=cos(15*x); I = trapz(x,y); v = [v; h,I,A-I];
end

得出结果如下：

可见，随着步距hhh的减小，计算精度逐渐增加。

matlab实现数值积分【一】（trapz函数）相关推荐

matlab进行数值积分的主要函数使用方法
MATLAB进行数值积分的主要函数: 1.trapz 梯形法求解积分 x=0:pi/10:pi; y=sin(x); trapz(x,y) 2.quad 基于变步长simpso法求积分 q = qua ...
matlab trapz 离散点,MATLAB求解方程之trapz函数
EDA365欢迎您登录! 您需要登录才可以下载或查看,没有帐号?注册 x ( T( B3 I- e% Q& H3 m trapz 是基于梯形法则的离散点积分函数. 调用形式:6 H* C! ...
matlab中 trapz函数,trapz(matlab中trapz函数用法)
采用梯形法求积分下面两种用法等效 h=pi/200; t=0:h:pi/2; y=sin(t)+cos(2*t); z1=trapz(t,y) z2=trapz(y)*h z1 = 1.0000 z ...
matlab trapz无穷,trapz(matlab中trapz函数用法)
采用梯形法求积分下面两种用法等效 h=pi/200; t=0:h:pi/2; y=sin(t)+cos(2*t); z1=trapz(t,y) z2=trapz(y)*h z1 = 1.0000 z ...
matlab中trapz的用,matlab使用trapz函数对积分的运算
Matlab中find函数使用输出是对的.因为他找到的是B(2:n,2)中的第一个元素.你可以好好看看find函数的作用.你如果写下面代码应该就对了:A=B(2:n,2);A(q=find(B(2: ...
matlab trapz x y,请问matlab trapz函数用法？？？？？？？
大野瘦子的回答: 使用方法如下: I=trapz(x,y) 其中x和y分别是自变量和对应的值,例如有函数y=x^3-2x-3,为了计算在[0,1]上的积分,可以这么做: >> format ...
trapz 函数在matlab 和 python中的区别
trapz 函数! 1.matlab 中,贝塞尔函数的表达形式他山之石--官方文档 // matlab中的函数 Q = trapz(Y) /* 如果 Y 为向量,则 trapz(Y) 是 Y 的近似 ...
matlab实现数值积分【二】（integral函数）
如果被积函数的数学表达式已知,但解析解不易求,可使用数值积分的方法求解积分. 目录函数调用格式应用举例例1:求解数值解并检验其精度例2:分段函数积分例3:与梯形法比较例4:大范围积分例5 ...
matlab多重数值积分,Matlab数值积分和微分(一重、多重都有)，简单实用。
第8章 MATLAB数值积分与微分 8.1 数值积分 8.2 数值微分 8.1 数值积分 8.1.1 数值积分基本原理求解定积分的数值方法多种多样,如简单的梯形法.辛普生(Simpson)法.牛顿- ...
2021-01-07 matlab数值分析数值积分与数值微分复合梯形公式复合Simpson公式
matlab数值分析数值积分与数值微分 1 复合梯形公式 function I=ftrapz(f,a,b,n) format long %显示15位双精度 h=(b-a)/n; x=linspace ...

matlab实现数值积分【一】（trapz函数）

目录

总述

函数调用格式

应用举例

例1：梯形法求积分

例2：不同步长对积分结果的影响

matlab实现数值积分【一】（trapz函数）相关推荐

最新文章

热门文章

matlab实现数值积分 【一】（trapz函数）

目录

总述

函数调用格式

应用举例

例1：梯形法求积分

例2：不同步长对积分结果的影响

matlab实现数值积分 【一】（trapz函数）相关推荐

最新文章

热门文章

matlab实现数值积分【一】（trapz函数）

matlab实现数值积分【一】（trapz函数）相关推荐