matlab函数作图指令

MATLAB 分段函数 绘图

1.一元分段函数绘图

例如：

把下面的函数保存为Piecewise_x.m文件

function F=Piecewise_x(x)

F=x.^2.*(x>=0 & x<1)+cos(pi*(x-1)).*(x>=1 & x<2)+(-x.^2./(x+2)).*(x>=2 & x<=4); end

运行：

x=linspace(0,4);

F=Piecewise_x(x);%计算相应函数值

plot(x,F);%绘制曲线

hold on;

plot(1*ones(1,2),ylim,'r:');%画区间间隔线

plot(2*ones(1,2),ylim,'r:');%画区间间隔线

一元分段函数图

2.二元分段函数绘图

例如：

把下面的函数保存为Piecewise_xy.m文件

function Pxy=Piecewise_xy(x,y)

Pxy=0.5457*exp(-0.75*y^2 - 3.75*x^2 - 1.5*x).*(x+y>1)+... 0.7575*exp(-y^2 - 6*x^2).*(x+y>-1)+...

0.5457*exp(-0.75*y^2 -

第 6 章 数据可视及探索

视觉是人们感受世界、认识自然的最重要依靠。数据可视化的目的：借助几何、色彩媒质表现一堆貌似杂乱的离散数据集合的形态，暴露数据内在关系和总体趋势，进而揭示出数据所传递的内在本质。随着计算机软硬件的发展，能力越来越强的图形表现，反过来对科学计算方法本身也产生了越来越大的影响。 针对符号计算和数值计算，MATLAB配置了两套形式上不同的绘图指令：“图形易绘指令”和“数值绘图指令”。前者配合符号计算，已安排在第5.8节；而本章内容将完全围绕数值绘图展开。 本章将系统地阐述：离散数据表示成图形的基本机理；曲线、曲面绘制的基本技法和指令；特殊图形的生成和使用示例；如何使用线型、色彩、数据点标记凸现不同数据的特征；如何利用着色、灯光照明、反射效果、材质体现和透明度处理渲染、烘托表现高维函数的性状；如何生成和运用标识，画龙点睛般地注释图形；如何表现变址、灰度、真彩图象；如何制作动画等。 随着MATLAB版本的升级，现今的MATLAB 图形窗不再是“单向性的图形显示工具”，而已成为进行“双向性探索的图形交互界面”。为此，占本章四分之一篇幅的第6.8节用于专门叙述全交互式绘图、图形对象属性的交互式设置、绘图用M函数文件的自

matlab命令

一、常用对象操作：除了一般windows窗口的常用功能键外。 1、!dir 可以查看当前工作目录的文件。 !dir& 可以在dos状态下查看。 2、who 可以查看当前工作空间变量名， whos 可以查看变量名细节。 3、功能键： 功能键 快捷键 说明 方向上键 Ctrl+P 返回前一行输入 方向下键 Ctrl+N 返回下一行输入 方向左键 Ctrl+B 光标向后移一个字符 方向右键 Ctrl+F 光标向前移一个字符 Ctrl+方向右键 Ctrl+R 光标向右移一个字符 Ctrl+方向左键 Ctrl+L 光标向左移一个字符 Home

Ctrl+A 光标移到行首 End Ctrl+E 光标移到行尾 Esc Ctrl+U 清除一行

Del Ctrl+D 清除光标所在的字符 Backspace Ctrl+H 删除光标前一个字符 Ctrl+K 删除到行尾 Ctrl+C 中断正在执行的命令

4、clc可以命令窗口显示的内容，但并不清除工作空间。

二、函数及运算 1、运算符：

＋：加， －：减， *：乘， /： 除， \\：左除 ^：幂， ‘：复数的共轭转置， （）：制定运算顺序

实验作业2 MATLAB作图

1、用plot,fplot 绘制函数y=cos(tan(?x))的图形。

解： 代码如下:

x=linspace(-0.5,0.5,30); y=cos(tan(pi*x)); plot(x,y) 结果：

fplot('[cos(tan(pi*x))]',[-0.5,0.5])

结果：

注意事项： fplot('[cos(tan(pi*x))]',[-0.5,0.5])不能写成fplot('[cos(tan(pi*x))]',（-0.5,0.5）)否则是错误程序。

2. 用ezplot绘制函数 exy-sin(x+y)=0 在[-3, 3]上的图形。 解： 代码如下：

ezplot('exp(x*y)-sin(x+y)=0',[-3,3])

3. 用matlab以动画的方式绘制出摆线 x=a(t-sint) y=a(1-cost) （a自己赋值）的渐屈线

clear; clc; close; a=1; syms t

x=a*(t-sin(t)); y=a*(1-cos(t));

ezplot(x,y,[0,2*pi]),grid on;hold on; dy=diff(y)/diff(x); dyy=d

%函数

function dxdt=differential(t,x)

dxdt=[x(1)*(1-x(1))*(5-6*x(2));x(2)*(1-x(2))*(1-6*x(1))]; end %主程序 clear %y-x

for i=0:0.1:1 for j=0:0.1:1

[T,Y]=ode45('differential',[0 5],[i j]); figure(1) grid on

plot(Y(:,1),Y(:,2)); hold on end end %dx/dt-t for i=0:0.1:1 for j=0:0.1:1

[T,Y]=ode45('differential',[0 5],[i j]); figure(2) grid on

plot(T,Y(:,1)); hold on end end %dy/dt-t for i=0:0.1:1 for j=0:0.1:1

[T,Y]=ode45('di

赚点财富值用一下 良心下载券

第二讲MATLAB作图

赚点财富值用一下 良心下载券

二维图形

三维图形

图形处理实例

特殊二、三维图形

作 业

赚点财富值用一下 良心下载券

1.曲线图 Matlab作图是通过描点、连线来实现的，故在 画一个曲线图形之前，必须先取得该图形上的 一系列的点的坐标(即横坐标和纵坐标),然 后将该点集的坐标传给Matlab函数画图.命令为： plot(X,Y,S)

线型

y 黄色 . 点 - 连线 m 洋红 o 圈 : 短虚线 c 蓝绿色 x x-符号 -. 长短线 r 红色 + 加号 -- 长虚线

X,Y是向量,分别表示点集的横坐标和纵坐标

plot(X,Y)--画实线 plot(X,Y1,S1,X,Y2,S2,……,X,Yn,Sn) --将多条线画在一起

赚点财富值用一下 良心下载券

图形元素设定颜色b 蓝(默认) g r y m c w k 绿 红 黄 洋红 青 白 黑

线形

标记

—— 实线(默认) 无标记(默认) * -划线 . -: 虚线 --. 点划线 。 x + s d 等

赚点财富值用一下 良心下载券

例 在[0,2*pi]用红线画sin(x),用绿圈画cos(x).解 x=linspace(0,2*pi,30)

%函数

function dxdt=differential(t,x)

dxdt=[x(1)*(1-x(1))*(5-6*x(2));x(2)*(1-x(2))*(1-6*x(1))]; end %主程序 clear %y-x

for i=0:0.1:1 for j=0:0.1:1

[T,Y]=ode45('differential',[0 5],[i j]); figure(1) grid on

plot(Y(:,1),Y(:,2)); hold on end end %dx/dt-t for i=0:0.1:1 for j=0:0.1:1

[T,Y]=ode45('differential',[0 5],[i j]); figure(2) grid on

plot(T,Y(:,1)); hold on end end %dy/dt-t for i=0:0.1:1 for j=0:0.1:1

[T,Y]=ode45('di

实验二 MATLAB编程与作图

1、设x为一个长度为n的数组，编程求下列均值和标准差

2?1n1?n2x??xi,s?x?nx?i?,n?1ni?1n?1??i?1?

并以x?[81,70,65,51,76,66,90,87,77]为例计算。 解： function [m,s]=myf_1(x) n=length(x); m=sum(x)/n;

s=sqrt((sum(x.^2)-n*m^2)/(n-1)); >>x=[81,70,65,51,76,66,90,87,61,77]; >> [m,s]=myf_1(x) m = 72.4000 s =

12.1124

2、求满足n?0?ln(1?n)?100m的最小m值。

解：>> s=log(1);n=0; >> while s<=100 n=n+1;

s=s+log(1+n); end >> m=n m =

37

3、用循环语句形成Fibonacci数列

F1?F2?1,Fk?Fk?1?Fk?2,k?3,4,L，并验证极限

limFk1?5?k??F2。 k?1解：>> clear;

>> F(1)=1;F(2)=1;k=2;x=0; >> e

实验二 MATLAB编程与作图

1、设x为一个长度为n的数组，编程求下列均值和标准差

2?1n1?n2x??xi,s?x?nx?i?,n?1ni?1n?1??i?1?

并以x?[81,70,65,51,76,66,90,87,77]为例计算。 解： function [m,s]=myf_1(x) n=length(x); m=sum(x)/n;

s=sqrt((sum(x.^2)-n*m^2)/(n-1)); >>x=[81,70,65,51,76,66,90,87,61,77]; >> [m,s]=myf_1(x) m = 72.4000 s =

12.1124

2、求满足n?0?ln(1?n)?100m的最小m值。

解：>> s=log(1);n=0; >> while s<=100 n=n+1;

s=s+log(1+n); end >> m=n m =

37

3、用循环语句形成Fibonacci数列

F1?F2?1,Fk?Fk?1?Fk?2,k?3,4,L，并验证极限

limFk1?5?k??F2。 k?1解：>> clear;

>> F(1)=1;F(2)=1;k=2;x=0; >> e

信源函数

randerr 产生比特误差样本

randint 产生均匀分布的随机整数矩阵 randsrc 根据给定的数字表产生随机矩阵 wgn 产生高斯白噪声 信号分析函数

biterr 计算比特误差数和比特误差率 eyediagram 绘制眼图 scatterplot 绘制分布图

symerr 计算符号误差数和符号误差率 信源编码

compand mu律/A律 压缩/扩张

dpcmdeco DPCM（差分脉冲编码调制）解码dpcmenco DPCM编码 dpcmopt 优化DPCM参数

lloyds Lloyd法则优化量化器参数 quantiz 给出量化后的级和输出值 误差控制编码

bchpoly 给出二进制BCH码的性能参数和产生多项式

convenc 产生卷积码

cyclgen 产生循环码的奇偶校验阵和生成矩阵

cyclpoly 产生循环码的生成多项式 decode 分组码解码器 encode 分组码编码器

gen2par 将奇偶校验阵和生成矩阵互相转换

gfweight 计算线性分组码的最小距离

hammgen 产生汉明码的奇偶校验阵和生成矩阵

rsdecof 对Reed-Solomon编码的ASCII文件解码

rsencof