MATLAB绘图和数据可视化课件.pptxVIP

专题1：数学物理建模思想与 数据可视化方法 数学物理建模与计算机辅助设计 (一) 数学物理方法问题建模的基本原理 . 主要解决的问题： . (1) 三维可视化方法使得“数学物理方法”中复杂、 抽象、烦琐的理论模型变得直观、具体、明了。 例1: 绘制复变函数 的图像，并验证其泰勒展开和罗朗展开的函数。 (一) 数学物理方法问题建模的基本原理 图3 在1|z|2内罗朗级数展开图形 图4 |z|2内罗朗级数展开图形 Page (一) 数学物理方法问题建模的基本原理 例2: 绘制复变量Γ(z)函数的图形 Page (一) 数学物理方法问题建模的基本原理 例3： 绘制球函数 的图形，取l=3，m=0 图6 球函数 的图形 Page (一) 数学物理方法问题建模的基本原理 例4：绘制四周固定的矩形膜本征振动的图形 图7 矩形膜振动的前四个本征函数三维可视化结果 Page (一) 数学物理方法问题建模的基本原理 主要解决的问题： . (2) 通过仿真设计与验证，解决前沿性的科研中缺乏实 图9 布拉格光子晶体光纤场振幅Hz三维分布 图8 布拉格光纤横截面结构示意图 验条件的问题。 Page (二) 利用MATLAB进行数学物理建模基础 例1：已知函数，如何求导及高阶导数？ 思路：① 问题：求导过程很繁杂，容易出错 sin(x)/(x^2+4*x+3)+4*cos(x)/(x^2+4*x+3)^2*(2*x+4)-12*sin(x)/(x^2 +4*x+3)^3*(2*x+4)^2+12*sin(x)/(x^2+4*x+3)^2-24*cos(x)/(x^2+ 4*x+3)^4*(2*x+4)^3+48*cos(x)/(x^2+4*x+3)^3*(2*x+4)+24*sin(x)/(x^2+4*x+3)^5*(2 *x+4)^4-72*sin(x)/(x^2+4*x +3)^4*(2*x+4)^2+ 24*sin(x)/(x^2+4*x+3)^3 使用Matlab的符号运算功能 syms x diff(sin(x)/(x^2+4*x+3),4) ans= 由分式求导公式，得出 ②逐次求导则可以得出 Page MATLAB的常量和变量 • Matlab预定义一些常量 Page 1 直接输入向量 a=[1 2 3 4 5 6 7 8 9] b=[1;2;3;4;5;6;7;8;9] – 用空格或逗号生成行向量,用分号生成列向量 2 利用冒号表达式基本形式:x=x1:step:xn a=1: 2: 12 a=1 3 5 7 9 11 3 线性等分向量生成 基本格式y=linspace(x1,x2,n) MATLAB的向量的生成 • 向量的生成的三种方式 Page MATLAB的向量的运算 A + B A - PBge 加减与数加减运算 • • 生 成 矩 阵 的 四 种 方 式 – 直 接 输 入 小 矩 Page MATLAB的矩阵的生成 阵 • 矩阵的四则运算 – 矩阵加减运算 + - • 规则：矩阵同阶(维数相同，各维度长度对应相同)， 对应做加减 A=[1 2 3;2 3 4;3 4 5]; B=[1 1 1;2 2 2;3 3 3]; C=A+B; MATLAB的矩阵的运算 Page MATLAB的矩阵的运算 矩阵的乘法运算 矩阵的四则运算 – Page • MATLAB的矩阵的运算 • 矩 阵 的 四 则 运 算 – 矩 阵 的 除 法 运 算 ， 分 Page 为 左 MATLAB的矩阵的运算 矩阵的行列式运算 d e t ( A P)age • • 矩阵的幂运算 A^n 和 A.^n – 规则：(1)A^n中A必须为方阵，相当于n个方阵A连乘 (2)A.^n中A不必为方阵，结果为A中各元素取n次幂 D=[1,2,3;4,5,6;7,8,9]; E=D^3; F=D.^3; MATLAB的矩阵的运算 Page ans = 5.0939 -0.0470 + 2.5448i -0.0470 - 2.5448i • 求多项式的乘除 conv(乘法) deconv(除法) p= [1 2 3]; q=[4 5 6 7]; pq=conv(p,q); deconv(pq,q); pq=[4 13 28 34 32 21] ans=[1 2 3] • 求多项式的