第六章 MATLAB IN ENGINEERING

Polynomial Differentiation多项式微分

 %幂级数
 f(x) = x^3-2x-5;
 p = [1 0 -2 -5]  %自变量前系数
 ​
 polyval()
 eg:
 a = [9, -5, 3, 7];
 x = -2:0.01:5;
 f = polyval(a,x);  %表示a系数下的X处的微商数值
 plot(x,f,'LineWidth',2);
 xlabel('x');
 ylabel('f(x)');
 set(gca,'FontSize',14)
 ​
 polyder()  %给的是微商式子
 eg:
 p = [5 0 -2 0 1];
 polyder(p)    %微分符号
 eg:
 p = [5 0 -2 0 1];
 polyder(p);
 polyval(polyder(p),7)

Polynomial Integration多项式积分

 polyint()  %不定积分
 eg:
 p = [5 0 -2 0 1];
 polyint(p,3);    %积分后的常数项是3
 polyval(polyint(p,3),7)    %不定积分7处的积分值

Numerical Differentiation

 diff()    %计算前后差异  
 eg1:
 x = [1 2 5 2 1];
 diff(x)   % ans = [1 3 -3 -1];
 eg2(两点斜率(1,5),(2,7)):
 x = [1 2];
 y = [5 7];
 slope = diff(y)./diff(x)
 ​
 eg3:
 x0 = pi/2;
 h = 0.1;
 x = [x0 x0+h];
 y = [sin(x0) sin(x0+h)];
 m = diff(y)./diff(x)
 eg4:
 h = 0.5;
 x = 0:h:2*pi;
 y = sin(x);
 m = diff(y)./diff(x);
 //h的大小对于导函数的影响:
 g = colormap(lines);
 hold on;
 for i = 1:4
  x = 0:power(10,-i):pi;
  y = sin(x);
  m = diff(y)./diff(x);
  plot(x(1:end-1),m,'Color',g(i,:));
 end
 hold off;
 set(gca,'XLim',[0,pi/2]);
 set(gca,'YLim',[0,1.2]);
 set(gca,'FontSize',18);
 set(gca,'FontName','symbol');
 set(gca,'XTick',0:pi/4:pi/2);
 set(gca,'XTickLabel',{'0','p/4','p/2'});
 h = legend('h = 0.1','h = 0.01','h = 0.001','h = 0.0001');
 set(h,'FontName','Times New Roman');
 box on;

Second and Third Derivatives

 x = -2:0.005:2;
 y = x.^3;
 m = diff(y)./diff(x);
 m2 = diff(m)./diff(x(1:end-1));
 ​
 plot(x,y,x(1:end-1),m,x(1:end-2),m2);
 xlabel('x','FontSize',18);
 ylabel('y','FontSize',18);
 legend('f(x) = x^3','f''(x)','f''''(x)');
 set(gca, 'FontSize',18);

Numerical Integration

Midpoint Rule用矩形去逼近

 h = 0.05;
 x = 0:h:2;
 midpoint = (x(1:end-1)+x(2:end))./2;
 y = 4*midpoint.^3;   %f(x) = 4*x^3;
 s = sum(h*y)

Trapezoid Rule用梯形去逼近

 trapz();%梯形数值积分
 eg:
 h = 0.05;
 x = 0:h:2;
 y = 4*x.^3;
 s = h*trapz(y);  
 OR:
 h = 0.05;
 x = 0:h:2;
 y = 4*x.^3;
 trapezoid = (y(1:end-1)+y(2:end))/2;
 s = h*sum(trapezoid);

Simpson's Rule

 h = 0.05;
 x = 0:h:2;
 y = 4*x.^3;
 s = h/3*(y(1)+2*sum(y(3:2:end-2))+4*sum(y(2:2:end))+y(end))

Funtion Handles(@)

function [y] = xy_plot(input,x)   %丢入类似@sin ; 一个函数不能作为input
% xy_plot receives the handle of a function and plots that
%function of x
y = input(x);
plot(x,y,'r--');
xlabel('x');
ylabel('function(x)');
end integral();
eg:
y = @(x) 1/(x.^3-2*x-5);
integral(y,0,2) %0,2是上下限

Double and Triple Integrals

eg1:
f = @(x,y) y.*sin(x)+x.*cos(y);
integral2(f,pi,2*pi,0,pi) %二重积分,顺序跟书写的一样 eg2:
f = @(x,y,z) y.*sin(x)+z.*cos(y);
integral3(f,0,pi,0,1,-1,1);

Matlab %陆的更多相关文章

  1. Matlab 绘制三维立体图(以地质异常体为例)

    前言:在地球物理勘探,流体空间分布等多种场景中,定位空间点P(x,y,x)的物理属性值Q,并绘制三维空间分布图,对我们洞察空间场景有十分重要的意义. 1. 三维立体图的基本要件: 全空间网格化 网格节 ...

  2. Matlab slice方法和包络法绘制三维立体图

    前言:在地球物理勘探,流体空间分布等多种场景中,定位空间点P(x,y,x)的物理属性值Q,并绘制三维空间分布图,对我们洞察空间场景有十分重要的意义. 1. 三维立体图的基本要件: 全空间网格化 网格节 ...

  3. Matlab 高斯_拉普拉斯滤波器处理医学图像

    前言:本程序是我去年实现论文算法时所做.主要功能为标记切割肝脏区域.时间有点久,很多细节已经模糊加上代码做了很多注释,因此在博客中不再详述. NOTE: 程序分几大段功能模块,仔细阅读,对解决医学图像 ...

  4. MATLAB中绘制质点轨迹动图并保存成GIF

    工作需要在MATLAB中绘制质点轨迹并保存成GIF以便展示. 绘制质点轨迹动图可用comet和comet3命令,使用例子如下: t = 0:.01:2*pi;x = cos(2*t).*(cos(t) ...

  5. linux下配置matlab运行环境(MCR)

    在安装好的matlab下有MCR(MatlabCompilerRuntime)在matlab2011/toolbox/compiler/deploy/glnxa64下找到MCRInstaller.zi ...

  6. EMD分析 Matlab 精华总结 附开源工具箱(全)

    前言: 本贴写于2016年12与15日,UK.最近在学习EMD(Empirical Mode Decomposition)和HHT(Hilbert-Huang Transform)多分辨信号处理,FQ ...

  7. Atitit MATLAB 图像处理 经典书籍attilax总结

    Atitit MATLAB 图像处理 经典书籍attilax总结 1.1. MATLAB数字图像处理1 1.2. <MATLAB实用教程(第二版)>((美)穆尔 著)[简介_书评_在线阅读 ...

  8. Atitit MATLAB 图像处理attilax总结

    Atitit MATLAB 图像处理attilax总结 1.1. 下载 Matlab7.0官方下载_Matlab2012 v7.0 官方简体中文版-办公软件-系统大全.html1 1.2. Matla ...

  9. Atitit java c# php c++ js跨语言调用matlab实现边缘检测等功能attilax总结

    Atitit java c# php c++ js跨语言调用matlab实现边缘检测等功能attilax总结 1.1. 边缘检测的基本方法Canny最常用了1 1.2. 编写matlab边缘检测代码, ...

  10. 使用MATLAB对图像处理的几种方法(下)

     试验报告 一.试验原理: 图像点处理是图像处理系列的基础,主要用于让我们熟悉Matlab图像处理的编程环境.灰度线性变换和灰度拉伸是对像素灰度值的变换操作,直方图是对像素灰度值的统计,直方图均衡是对 ...

随机推荐

  1. Shell-通过expect工具脚本的自动交互及实例

    安装expect工具 expect是建立在tcl基础上的一个自动化交互套件, 在一些需要交互输入指令的场景下, 可通过脚本设置自动进行交互通信. 其交互流程是: spawn启动指定进程 -> e ...

  2. sap 付费支持 fico付费求助

    从事SAP多年,SAP付费求助,fico有偿服务,月结问题 有偿处理,物料分类账异常处理 每次每个问题最少 500 CNY , 有需要的联系 wx :erpworld sap  fico  有偿服务

  3. TS 之 reduce

    一.函数介绍 Array.reduce()方法是对数组进行遍历,返回一个计算后的值 使用方法: Array.reduce((acc, cur, idx, src) => { }, initial ...

  4. JDBC之Statement

    Statement 目录 Statement Statement Statement概述 要执行的SQL分为两类 查询 增删改 Statement继承体系 SQL注入问题 SQL注入问题解决 获取得到 ...

  5. java8 检查

    // 检查放款日期 boolean allMatch = cdiscountMonthBillDetails.stream().map(CdiscountMonthBillDetail::getTra ...

  6. 计数 dp 部分例题(十一~十五部分)

    十一.矩阵的利用(行列を用いたテクニック) 1. 快速幂(二分累乗) (1) 推导转移矩阵(行列の導出) 例题:Placing Squares 题解 (2) BM 优化递推(?)(コンパニオン行列の累 ...

  7. 37.Spring注解相关面试题

    1.@Resource和@Autowired 区别 2.@Repository.@Component.@Service.@Controller 区别

  8. Python Cli 编写指南

    Python Cli 编写指南 python实现cli 环境: python 3.8 库 python自带argparse 指南 简单示例 : cli.py import argparse def c ...

  9. 【ADB命令】安装app

    在电脑上安装以下指令 adb install app的文件位置

  10. windows U盘启动装机工具(Ventoy)

    Windows优盘启动装机工具 一.WePE 纯净 WePE 启动盘/ U盘重装系统工具 1.U盘重装系统工具下载地址:https://www.wepe.com.cn/download.html?td ...