f=@(z1,z2)exp(abs(z1-z2));%对该函数进行积分,积分范围0

quad2d(@(z1,z2)f(z1,z2),0,3,1,3)

ans =

19.7563

>> %去掉绝对值号计算,对f函数的积分转化成对f1,f2两个函数的积分

f1=@(z1,z2)exp(z1-z2);%积分范围0

f2=@(z2,z1)exp(z2-z1);%积分范围1

c=@(z1)z1;

d=@(z2)z2;

quad2d(@(z1,z2)f1(z1,z2),0,3,1,c)+quad2d(@(z2,z1)f2(z2,z1),1,3,0,d)

ans =

19.3884

末学已经确定第二种方法是错误的,因为对f1的积分quad2d(@(z1,z2)f1(z1,z2),0,3,1,c)来说,并不表示1

但不知道第一种方法是否正确?

如果要去掉绝对值号,应该怎么改写这个函数?

==================================================

还有一个问题,对函数f4积分

f4=@(z1,z2)exp(z1-z2);

quad2d(@(z1,z2)f4(z1,z2),2,5,3,8)

ans =

6.9739

但我用分段求和的方法来积分这个函数,结果如下

sum1=0;

N=1000;       %分段数

d1=(5-2)/N;  %变量z1的步长

d2=(8-3)/N;  %变量z2的步长

for i=1:N

sum2=0;

z1=(2+(i-1)*d1+2+i*d1)/2;   %先固定z1,对z2积分,sum2表示对z2的积分结果,sum1表示对z1的积分结果

for j=1:1000

z2=(3+(i-1)*d2+3+i*d2)/2;

sum2=sum2+f4(z1,z2)*d2;

end

sum1=sum1+sum2*d1;

end

>> sum1

sum1 =

2.3857

这是什么原因,哪里有错误吗,还是精度产生的误差?末学对matlab非常外行,请大贤指点,感激涕零!

matlab 取绝对值最快的犯法,求助被积函数有绝对值号的问题。相关推荐

  1. MATLAB仿真QPSK调制信号通过AWGN信道的误符号率和误比特率分析

    MATLAB仿真QPSK调制信号通过AWGN信道的误符号率和误比特率分析 形式:程序 程序实现功能: 仿真正交相移键控QPSK信号调制的基带数字通信系统通过AWGN信道的误符号率(SER)和误比特率( ...

  2. matlab无限循环怎么停下来,[求助] 关于matlab无限循环的问题

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 各路大神,求助啊啊啊啊~ 以下是我的代码: n=0 for i=1:length(NAXLES) if CODE(i)==0 n=n+1 GVM1(n)= ...

  3. matlab 12位 显示不出来,求助大神,为何不同机器运行MATLAB结果不同

    求助:不同机器运行MATLAB结果不同 我调用MATLAB优化工具箱的库函数fmincon,使用相同的初始解(可行解),对同一个问题进行局部搜索(算法为序列二次规划,即SQP),但在不同机器上得到的结 ...

  4. matlab安装无效距离过远,求助matlab的远程序

    求助matlab的远程序 function varargout = a1(varargin) % A1 M-file for a1.fig % A1, by itself, creates a new ...

  5. matlab求函数的无穷积分,[求助]贝塞尔函数的无穷积分怎么求解?

    回复: [求助]贝塞尔函数的无穷积分怎么求解? 对于这种积分格式,上限是正无穷,下限是zero,我编了个函数文件,回头你自己把函数的表达式代入就可以了. function I = IntGaussLa ...

  6. matlab gpu cpu 快,(Matlab)GPU计算简介,及其与CPU计算性能的比较

    1.GPU与CPU结构上的对比 原文: Multicore machines and hyper-threading technology have enabled scientists, engin ...

  7. matlab 中调用s函数表达式,[求助]S函数中能否调用M函数

    [求助]S函数中能否调用M函数 我现在写了一个S-FUNCTION 里面调用了一些M函数,当把S函数的名字写入S-FUNCTION块的时候 总是出现 ??? SWITCH expression mus ...

  8. MATLAB安装c 错误,matlab编译c文件错误,求助

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 >> mex -setup Please choose your compiler for building external interfa ...

  9. matlab里excel汉字怎么显示,[求助]matlab如何导入excel中的汉字,求助~~~

    第一种方法,你可以使用xlsread函数来读取excel中的数据 第二种方法,就是把字符转化为数字,使用函数str2num xlsread的使用 EXAMPLES: 1. Default operat ...

最新文章

  1. 收藏:SqlServer
  2. taro 在components文件夹中 新建组件时,组件支持自定义命名,但是不能大写开头...
  3. 创新创新再创新(4)
  4. 电子书下载:The C# Programming Language, 4th Edition
  5. activiti候选人的多个场景应用
  6. css如何实现背景透明,文字不透明?
  7. vm 安装jdk1.8_JDK 13:VM.events已添加到jcmd
  8. mysql数据库open函数_C#访问Mysql数据库方法,以及库函数
  9. 极大似然估计、拉普拉斯平滑定理、M-估计详解
  10. 新荣耀员工现金补偿最高可拿 N+5;天津立法禁止采集人脸识别信息;IntelliJ IDEA 新版发布|极客头条...
  11. swift 脚本_Swift脚本和命令行参数
  12. 世界各国国家代码简称 - 备用
  13. 基于51单片机智能车流量检测车量统计系统设计
  14. 2021-5-25有限元从0开始第六天(Lp空间,弱导数)
  15. 高德地图定位、画线 基础功能
  16. Windows 10 修改 Alt+Tab 键 切换 新版 Microsoft Edge 单个标签页窗口
  17. 证件照制作小程序源码下载-超级实用工具
  18. npm 配置项registry修改为淘宝镜像
  19. 文献管理软件 Endnote 简易上手
  20. 规范化理论:如何求属性集X关于F的闭包?

热门文章

  1. 算法高级(10)-如何实现浏览器的前进、后退功能
  2. 数据不平衡问题及解决方案
  3. mysql 编辑距离 搜索_最小编辑距离(Minimum Edit Distance)
  4. ajax then jquery,使用Jquery.ajax()。then()时无法.catch()错误
  5. react使用less预编译语言和本地代理配置
  6. 【Demo】HTML5 拍照上传
  7. easyui汉化啊!
  8. myeclipse8.5打包jar并引入第三方jar包
  9. PHP高效获取远程图片尺寸和大小(转)
  10. Chapter 15 配置服务器存储和群集 第1课