matlab中trapz,MATLAB Trapz
MATLAB函数trapz(x, y, n), 其中y是x的积分, 使用梯形法则逼近函数y = f(x)的积分, 并且n(可选)与维度n进行积分。
句法
Z = trapz(Y)
Z = trapz (X, Y)
Z = trapz (..., dim)
例1
使用MATLAB函数trapz(x, y)估算积分成本
并与精确值进行比较, 评估n = 5和n = 10时的百分比误差。
解
确切的值从
对于使用梯形法则的逼近, 我们让x5代表n = 5的行向量, x10代表n = 10的向量, 即∆x = 1/5和∆x = 1/10。相应的值表示为y5和y10, 曲线下的面积分别表示为area5和area10。
创建以下脚本
x5=linspace (1, 2, 5);
x10=linspace (1, 2, 10);
y5=1. /x5; y10=1. /x10;
area5=trapz (x5, y5), area10=trapz (x10, y10)
MATLAB显示以下结果:
area5 =
0.6970
area10 =
0.6939
使用∆x = 1/5时的百分比误差为
使用∆x = 1/10时的百分比误差为
例2
积分
其中τ是积分的虚拟变量, 称为误差函数, 在通讯理论中得到广泛使用。当积分的上限为t = 2时, 使用MATLAB trapz(x, y)函数在n = 10的情况下找到该积分下的面积。
解
创建以下脚本
t=linspace (0, 2, 10);
y=exp(-t.^2);
area=trapz (t, y)
MATLAB显示以下结果:
area =
0.8818
范例3
非线性电机的i-v(电流-电压)关系由下式给出
其中v(t)= sin3t。
无论如何, 找到
瞬时功率为
p(t)= v(t)i(t)= 0.1 sin3t(e0.2sin3t-1)
能量是瞬时能量的整数
使用零件积分可以对最后一个积分进行解析, 但这并不容易。我们可以尝试MATLAB int(f, a, b)函数, 其中f是符号表达式, 而a和b分别是积分的下限和上限。
当MATLAB无法找到解决方案时, 它将返回警告。例如, 当使用等式的符号表达式尝试积分时, MATLAB返回以下消息
t=sym('t');
s=int (0.1*sin(3*t) *(exp (0.2*sin(3*t))-1), 0, 10)
执行此脚本后, MATLAB将显示以下消息。
Warning: Explicit integral could not be found.
接下来, 我们将使用MATLAB trapz(x, y)函数通过梯形法则找到并绘制功率和能量。在此示例中, 我们选择n = 100, 因此∆x = 1/100。下面的MATLAB脚本将计算和绘制功率。
t=linspace (0, 10, 100);
v=sin (3. *t); i=0.1. *(exp (0.2.*v)-1); p=v.*i;
plot(t, p); grid; title('Power vs Time'); xlabel('seconds'); ylabel('watts')
功率以统一的方式变化, 如图所示
下面的MATLAB脚本计算并绘制能量。
energy=trapz (t, p), plot (t, energy, '+'); grid; title ('Energy vs Time'); ...
xlabel('seconds'); ylabel('joules')
energy =
0.1013
因此, 能量的值为0.1013焦耳。能量如图所示:
matlab中trapz,MATLAB Trapz相关推荐
- matlab中数值积分函数trapz的使用注意
在matlab中,使用数值积分的时候,应当注意,积分边界是否会使得被积函数产生无意义的值(无穷大). 用函数trapz举例说明. 此处,我们是求积分 ∫0111−x2dx\int_0^1{1 \ove ...
- matlab中dblquad,MATLAB积分.ppt
MATLAB积分 (2) 用int命令求解: 先将重积分化成累次积分: MATLAB代码为 clear; syms x y z; int(int(int('y*sin(x)+z*cos(x)',z,- ...
- matlab中normfit,MATLAB中如何得到一组统计数据的分布特征
MATLAB中如何得到一组统计数据的分布特征 我想通过MATLAB知道一组统计数据的分布特征,如属于哪种分布,均值方差等,请问如何实现?谢谢! function f=p_judge(A,alpha) ...
- matlab中fdyn,Matlab的用法总结
1. 对序列进行洗牌 randperm() randperm()产生随机的序列 %if filepaths 是一个5*1的结构体,then cshuffle = randperm(length(fil ...
- matlab中rsa,matlab - Matlab中的RSA代码 - 堆栈内存溢出
我想加密一条消息,例如" HELO1234" ,然后解密以获取原始消息.我已经在matlab中编写了无法正常工作的RSA代码. 参数计算 temp=1; range=1:10; k ...
- 在matlab中function,Matlab中function函数使用操作方法
亲们或许不知道Matlab中function函数如何使用,那么今天小编就讲解Matlab中function函数使用操作方法哦,希望能够帮助到大家呢. Matlab中function函数使用操作方法 打 ...
- matlab中投影,MATLAB在极射赤平投影中的应用
文章编号: 100926825 (2010) 360357202 MATLAB在极射赤平投影中的应用 收稿日期: 20100822 作者简介:潘冀川 (1988) ,男 ,石家庄经济学院本科生 ,河北 ...
- matlab中lambertw,MATLAB解常微分方程
在MATLAB中,由函数dsolve()解决常微分方程(组)的求解问题,其具体格式如下: X=dsolve('eqn1','eqn2',-) 函数dsolve用来解符号常微分方程.方程组,如果没有初始 ...
- tyvector在matlab中代表,MATLAB曲线绘制
信号源产生的方法 来源:http://www.2cto.com/kf/201401/270494.html matlab的checkerboard说明,GOOD! 来源:http://www.chi ...
- matlab中facedetector,Matlab人脸检测方法(Face Parts Detection)详解
今天同学让我帮忙制作一个人脸表情识别的样本库,其中主要是对人脸进行裁剪,这里用到了一个相对较新的Matlab人脸检测方法Face Parts Detection,网上百度了一下发现关于Matlab人脸 ...
最新文章
- ARCGIS中某字段递增赋值
- cordova使用cordova-plugin-baidumaplocation插件获取定位
- 风控评分模型全流程的开发及应用
- Django(part42)--限制IP请求次数
- mysql 常用命令与备份恢复 整理
- Kubernetes 1.20 版本开始将弃用 Docker,是时候拥抱 Containerd 和 Podman 了!
- php隐藏IP最后位,替换手机号中间数字为*号
- CraftAR入门教程
- 20190829:(leetcode习题)环形链表
- mysql批量修改http为https,搜索和将数据库中的“ http”替换为“ https”
- Atitit 单项功能开发 最佳实践规范 标准化流程attilax总结.docx
- 史上最全wireshark使用教程,8万字整理总结,建议先收藏再耐心研读
- Centos 防火墙开启/关闭端口
- 通达信接口的开发工具?
- C++ Primer 中文版(第 5 版)练习解答合集
- Golang 解析XML
- 最新最全的阿里云产品手册出炉
- 一个简单的出货补货网页php
- 【OpenCV】 ⚠️实战⚠️ 银行卡卡号读取 ☢️建议手收藏☢️
- java 取一个数的各个位数