第5章 数值计算

  • (5.1) 一重积分
  • (5.2) 多重积分
  • (5.3) 函数零点
  • (5.4) 函数极值
  • (5.5) 一维插值
  • (5.6) 多维插值

(5.1) 一重积分

clc;clear;
x = linspace(0,2*pi,100);
y = sin(x);
I = trapz(x,y) %9.1073e-17,梯形积分法
clc;clear;
f = @(x) sin(x).*exp(x);
I_Simpson = quad(f,0,2*pi); %辛普森(Simpson)法
I_Guass_Lobatto = quadl(f,0,2*pi); %高斯-洛巴托(Gauss-Lobatto)法
I_Guass_Cronrod = quadgk(f,0,2*pi); %高斯-克朗罗德(Gauss-Cronrod)法
I_Integral = integral(f,0,2*pi); %全局自适应积分法

(5.2) 多重积分

I = integral2(f,a,b,c,d);
I = quad2d(f,a,b,c,d);
I = dblquad(f,a,b,c,d);

I = integral3(f,a,b,c,d,e,f);
I = triplequad(f,a,b,c,d,e,f);

(5.3) 函数零点

clc;clear;
f = @(x) cos(x);
zero_point1 = fsolve(f,1) %求f在3附近的零点
zero_point2 = fsolve(f,[1,2]) %求f在区间[1 2]内的零点,区间必须符合零点存在定理
clc;clear;
f = @(x) [ sin(x(1)) + x(2) + x(3).^2 .*exp(x(1)),x(1) + x(2) + x(3),x(1) .*x(2) .* x(3) ];
x = fsolve(f,[1,1,1]) %求方程组在(1,1,1)附近的解

(5.4) 函数极值

x = fminsearch(f,x0); %在点x0附近寻找f的极小值点x
x = fminbnd(f,x1,x2); %在区间[x1,x2]内寻找f的极小值点x
[x,fmin] = fminsearch(f,x0); %在点x0附近寻找f的极小值点x和极小值fmin
[x,fmin] = fminbnd(f,x1,x2); %在区间[x1,x2]内寻找f的极小值点x和极小值fmin

(5.5) 一维插值

clc;clear;
x = 0:0.2:2;
y = (x.^2 - 3*x + 5).*exp(-3*x).*sin(x);
xi = 0:0.003:2;
yi_linear = interp1(x,y,xi); %线性插值
yi_nearest = interp1(x,y,xi,'nearest'); %临近点插值
yi_spline = interp1(x,y,xi,'spline'); %三次样条插值
yi_pchip = interp1(x,y,xi,'pchip'); %分段三次Hermite插值
yi_v5cubic = interp1(x,y,xi,'v5cubic'); %MATLAB5中的三次多项式插值
figure;
subplot(231);
plot(x,y,'rp');
title('数据点');
subplot(232);
plot(x,y,'rp',xi,yi_linear,'g-');
title('线性插值');
subplot(233);
plot(x,y,'rp',xi,yi_nearest,'g-');
title('临近点插值');
subplot(234);
plot(x,y,'rp',xi,yi_spline,'g-');
title('三次样条插值');
subplot(235);
plot(x,y,'rp',xi,yi_pchip,'g-');
title('分段三次Hermite插值');
subplot(236);
plot(x,y,'rp',xi,yi_v5cubic,'g-');
title('MATLAB5中的三次多项式插值');

(5.6) 多维插值

clc;clear;
x = linspace(-3,3,10);
y = linspace(-3,3,10);
[X,Y] = meshgrid(x,y);
Z = peaks(X,Y);
xi = linspace(-3,3,20);
yi = linspace(-3,3,20);
[Xi,Yi] = meshgrid(xi,yi);
Zi_linear = interp2(X,Y,Z,Xi,Yi); %线性插值
Zi_nearest = interp2(X,Y,Z,Xi,Yi,'nearest'); %临近点插值
Zi_spline = interp2(X,Y,Z,Xi,Yi,'spline'); %三次样条插值
Zi_cubic = interp2(X,Y,Z,Xi,Yi,'cubic'); %三次多项式插值
figure;
hold on;
colormap(jet);
subplot(3,2,[1 2]);
surf(X,Y,Z);
title('原始图');
subplot(323);
surf(Xi,Yi,Zi_linear);
title('线性插值');
subplot(324);
surf(Xi,Yi,Zi_nearest);
title('临近点插值');
subplot(325);
surf(Xi,Yi,Zi_spline);
title('三次样条插值');
subplot(326);
surf(Xi,Yi,Zi_cubic);
title('三次多项式插值');


《MATLABSYNTAX》系列博客创作参考资料来源《MATLAB\ SYNTAX》系列博客创作参考资料来源《MATLAB SYNTAX》系列博客创作参考资料来源

  1. 《自动控制原理实验教程》.巨林仓.西安交通大学出版社.
  2. 《MATLAB工程与科学绘图》.周博.薛世峰.清华大学出版社.
  3. 《MATLAB R2018a完全自学一本通》.刘浩.韩晶.电子工业出版社.
  4. 《科学计算与MATLAB语言》.刘卫国.蔡旭晖.吕格莉.何小贤.中国大学MOOC.
  5. 《MATLAB软件与基础数学实验》.李换琴.朱旭.王勇茂.籍万新.西安交通大学出版社.
  6. 《Matlab教程 - 图像处理》@正月点灯笼.https://www.bilibili.com.
  7. 《MATLAB从入门到秃头》@古德谓尔.https://www.bilibili.com.

博客创作:AidenLee博客创作:Aiden\ Lee博客创作:Aiden Lee
特别声明:文章仅供学习参考,转载请注明出处,严禁盗用!

《MATLAB SYNTAX》第5章 数值计算相关推荐

  1. 利用matlab来计算双曲函数的值,关于MATLAB第2章数值计算与数据分析的介绍

    关于MATLAB第2章数值计算与数据分析的介绍 MATLAB6.0 数学手册62第 2 章 数值计算与数据分析2.1 基本数学函数2.1.1 三角函数与双曲函数函数 sin.sinh功能 正弦函数与双 ...

  2. 《MATLAB SYNTAX》第4章 多项式

    第4章 多项式 (4.1) 多项式的建立 (4.2) 多项式的求值 (4.3) 多项式的求根 (4.4) 多项式的乘法 (4.5) 多项式的导数 (4.6) 多项式的积分 (4.7) 多项式的展开 ( ...

  3. 《MATLAB SYNTAX》第1章 数据类型

    第1章 数据类型 (1.1) 常用函数 (1.2) 数值型 (1.3) 字符(串)型 (1.4) 复数型 (1.5) 逻辑型 (1.6) 函数句柄 (1.7) 单元数组 (1.8) 结构体 (1.1) ...

  4. 《MATLAB SYNTAX》第9章 三维图形可视化

    第9章 三维图形可视化 (9.1) plot3 (9.2) fplot3 (9.3) scatter3 (9.4) pie3 (9.5) bar3 (9.6) bar3h (9.7) stem3 (9 ...

  5. Matlab抽样和,《统计学原理》与MATLAB编程-第三章 抽样和抽样分布

    <统计学原理>与MATLAB编程-第三章 抽样和抽样分布 第一节 排列与组合 排列: perms(x) x为向量,求x的全排列.如: a=perms([2 3 7 ]) a= 7 3 2 ...

  6. 自适应滤波器在matlab仿真的程序_电气信息类专业课程之matlab系统仿真 第五章 BPSK通信系统(3)...

    继续讲解! 上一篇文章提出了那么多问题, 不知道大家是否能回答上来啊? 如果回答不了, 有没有去问问度娘呢? 程序写完了,回想了一下辅导2017届学生毕业设计的过程,那是真累.我要求他们在现有程序基础 ...

  7. matlab高级数据结构,Matlab讲义-第四章Matlab语言数据结构1-10.pdf

    Matlab讲义-第四章Matlab语言数据结构1-10.pdf 第四章第四章 Matlab语言基础语言基础 Matlab语言程序简洁.可语言程序简洁.可 读性很强而且调试十分容易.读性很强而且调试十 ...

  8. 将数学模型用matlab,matlab课件第3章_控制系统的数学模型及其转换.ppt

    matlab课件第3章_控制系统的数学模型及其转换 3.1.1 传递函数 传递函数输入举例 3.1.2 零极点增益形式 3)零极点增益模型: 3.1.3 部分分式形式 2.1.4 状态空间表达式 3. ...

  9. 控制系统数学模型的matlab仿真,第7章 控制系统的MATLAB仿真

    <第7章 控制系统的MATLAB仿真>由会员分享,可在线阅读,更多相关<第7章 控制系统的MATLAB仿真(101页珍藏版)>请在人人文库网上搜索. 1.1,本章主要教学内容在 ...

最新文章

  1. Centos7上安装docker 详细教程
  2. 超越英伟达的,不会是另一款GPU!中国公司发布首款数据流AI芯片
  3. C#9引入的自定义类型record
  4. 【图论】清理牛棚/Cleaning Shifts S(luogu 4644)
  5. 提升用户体验,你不得不知道的事儿——三种提醒框的微技巧
  6. 【pytorch】torch.meshgrid()==>常用于生成二维网格,比如图像的坐标点
  7. android视频压缩框架,GitHub - tangpeng/VideoCompressor: Android 使用自带的MediaCodec 框架进行本地视频压缩,速度嗖嗖的,亲测有效!!!...
  8. C#发现之旅第八讲 ASP.NET图形开发带超链接的饼图
  9. python征程3.1(列表,迭代,函数,dic,set,的简单应用)
  10. google 翻译 api
  11. UltraEdit 25注册机及免费破解注册教程(附带工具)
  12. html做自我介绍模板,应聘自我介绍模板简短
  13. 【问题解决方案】电脑版微信如何使用外接摄像头
  14. mariadb mysql.h_MariaDB(MySQL)的常用命令1 【检索数据】
  15. 趣味js之斗牛算法详解
  16. Harpoon: an OSINT / Threat Intelligence tool
  17. python自定义函数实例 lite_python 函数小实例
  18. 考研英语近义词与反义词·十一
  19. Java执行Excel公式
  20. 什么是创建对象java_java创建对象的几种方式分别是什么

热门文章

  1. 一文看懂三极管在电路中的应用原理
  2. html中display有哪些属性值,css中display属性是什么及其各个值的作用和用法
  3. JAVA快速获取网络图片或者URL图片并保存到本地
  4. 【课件制作软件】Focusky教程 | 对齐工具快刀斩乱麻
  5. edit中网格线 ultra_【常用】Ultra Edit常用设置及快捷键
  6. R 评分卡 scorecard包
  7. ABOUT ME/OI回忆录
  8. iPhone 4.3.3 重刷 4.3.3 固件教程 及 一键越狱
  9. 微信红包php算法,基于PHP微信红包的算法探讨
  10. python3封装pymysql,方便快速调用