matlab进行数据处理

matlab入门经典范例

num1=[13]; den1=conv([1,1],[1,0]); G1=tf(num1,den1); num2=[5.096,13]; ssys1=conv([1,1],[1,0]); ssys2=conv([0.098,1],[1]); den2=conv( ssys1,ssys2); G2=tf(num2,den2); figure(1) margin(G1); hold on margin(G2); num1=[13]; den1=conv([1,1],[1,0]); G1=tf(num1,den1); num2=[5.096,13]; ssys1=conv([1,1],[1,0]); ssys2=conv([0.098,1],[1]); den2=conv( ssys1,ssys2); G2=tf(num2,den2); figure(1) margin(G1); hold on margin(G2); num=[4.56,10]; ssys1=conv([1,1],[1,0]); ssys2=conv([0.114,1],[1]); den=conv( ssys1,ssys2); G=tf(num,den); figure(1) bode(G) num=[4.56,10]; ssys1=conv([1,1],[1,0]); ssys2=conv([0.114,1],[1]); den=conv( ssys1,ssys2); G=tf(num,den); figure(1) msrgin(G); num1=[13]; den1=conv([1,1],[1,0]);

matlab经典习题及解答

第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 页脚内容1

命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？ 在MATLAB中有多种获得帮助的途径： （1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器； （2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息； （3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数 （4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息。 第2章MATLAB矩阵运算基础 页脚内容2

BP神经网络matlab实例(简单而经典)

p=p1';t=t1'; [pn,minp,maxp,tn,mint,maxt]=premnmx(p,t); %原始数据归一化 net=newff(minmax(pn),[5,1],{'tansig','purelin'},'traingdx'); %设置网络,建立相应的BP网络 net.trainParam.show=2000; % 训练网络 net.trainParam.lr=0.01; net.trainParam.epochs=100000; net.trainParam.goal=1e-5; [net,tr]=train(net ,pn,tn); %调用TRAINGDM算法训练BP网络 pnew=pnew1'; pnewn=tramnmx(pnew,minp,maxp); anewn=sim(net,pnewn); %对BP网络进行仿真 anew=postmnmx(anewn,mint,maxt); %还原数据 y=anew'; 1、BP网络构建 （1）生成BP网络 = (,[1 2...],{ 1 2...},,,) net newff PR S S SNl TF TF TFNl BTF BLF PF R?维矩阵。 PR：由R维的输入样本最小最大值构成的2 S S SNl：各层的神经元个数。 [ 1 2...] TF TF TFNl：各层的神经元传递函数。 { 1 2...} BTF：训练用函数的名称。 （2）网络训练 = net tr Y E Pf Af train net P T Pi Ai VV TV [,,,,,] (,,,,,,) （3）网络仿真 = [,,,,] (,,,,) Y Pf Af E perf sim net P Pi Ai T {'tansig','purelin'},'trainrp' BP网络的训练函数 训练方法训练函数 梯度下降法traingd 有动量的梯度下降法traingdm 自适应lr梯度下降法traingda 自适应lr动量梯度下降法traingdx 弹性梯度下降法trainrp Fletcher-Reeves共轭梯度法traincgf

matlab入门教程

MATLAB入门教程 1．MATLAB的基本知识 1-1、基本运算与函数 在MATLAB下进行基本数学运算，只需将运算式直接打入提示号（>>）之後，并按入Enter键即可。例如： >> (5*2+1.3-0.8)*10/25 ans =4.2000 MATLAB会将运算结果直接存入一变数ans，代表MATLAB运算後的答案（Answer）并显示其数值於萤幕上。 小提示： ">>"是MATLAB的提示符号（Prompt），但在PC中文视窗系统下，由於编码方式不同，此提示符号常会消失不见，但这并不会影响到MATLAB的运算结果。 我们也可将上述运算式的结果设定给另一个变数x： x = (5*2+1.3-0.8)*10^2/25 x = 42 此时MATLAB会直接显示x的值。由上例可知，MATLAB认识所有一般常用到的加（+）、减（-）、乘（*）、除（/）的数学运算符号，以及幂次运算（^）。 小提示：MATLAB将所有变数均存成double的形式，所以不需经过变数宣告（Variable declaration）。MATLAB同时也会自动进行记忆体的使用和回收，而不必像C语言,必须由使用者一一指定.这些功能使的MATLAB易学易用，使用者可专心致力於撰写程式，而不必被软体枝节问题所干扰。 若不想让MATLAB每次都显示运算结果，只需在运算式最後加上分号（；）即可，如下例： y = sin(10)*exp(-0.3*4^2); 若要显示变数y的值，直接键入y即可： >>y y =-0.0045 在上例中，sin是正弦函数，exp是指数函数，这些都是MATLAB常用到的数学函数。 下表即为MATLAB常用的基本数学函数及三角函数： 小整理：MATLAB常用的基本数学函数 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数z的相角(Phase angle)

matlab经典编程例题

以下各题均要求编程实现，并将程序贴在题目下方。 1．从键盘输入任意个正整数，以0结束，输出那些正整数中的素数。 clc;clear; zzs(1)=input('请输入正整数：');k=1; n=0;%素数个数 while zzs(k)~=0 flag=0;%是否是素数，是则为1 for yz=2:sqrt(zzs(k))%因子从2至此数平方根 if mod(zzs(k),yz)==0 flag=1;break;%非素数跳出循环 end end if flag==0&zzs(k)>1%忽略0和1的素数 n=n+1;sus(n)=zzs(k); end k=k+1; zzs(k)=input('请输入正整数：'); end disp(['你共输入了' num2str(k-1) '个正整数。它们是：']) disp(zzs(1:k-1))%不显示最后一个数0 if n==0 disp('这些数中没有素数！')%无素数时显示 else disp('其中的素数是：') disp(sus) end 2．若某数等于其所有因子（不含这个数本身）的和，则称其为完全数。编程求10000以内所有的完全数。 clc;clear;

wq=[];%完全数赋空数组 for ii=2:10000 yz=[];%ii的因子赋空数组 for jj=2:ii/2 %从2到ii/2考察是否为ii的因子 if mod(ii,jj)==0 yz=[yz jj];%因子数组扩展，加上jj end end if ii==sum(yz)+1 wq=[wq ii];%完全数数组扩展，加上ii end end disp(['10000以内的完全数为：' num2str(wq)])%输出 3．下列这组数据是美国1900—2000年人口的近似值（单位：百万）。 （1）若. 2c + = y+ 与试编写程序计算出上式中的a、b、c; 的经验公式为 t at bt y （2）若.bt 的经验公式为 y= 与试编写程序计算出上式中的a、b; y ae t （3）在一个坐标系下，画出数表中的散点图（红色五角星），c + =2中 ax bx y+拟合曲线图(蓝色实心线),以及.bt y=(黑色点划线)。 ae （4）图形标注要求:无网格线,横标注“时间t”，纵标注“人口数(百万)”，图形标题“美国1900—2000年的人口数据”。 （5）程序中要有注释，将你的程序和作好的图粘贴到这里。 clf;clc;clear %清除图形窗、屏幕、工作空间 t=1900:10:2000; y=[76 92 106 123 132 151 179 203 227 250 281]; p1=polyfit(t,y,2);%二次多项式拟合

Matlab基础教程

1-1、基本运算与函数 在MATLAB下进行基本数学运算，只需将运算式直接打入提示号（>>）之後，并按入Enter键即可。例如： >> (5*2+1.3-0.8)*10/25 ans =4.2000 MATLAB会将运算结果直接存入一变数ans，代表MATLAB运算後的答案（Answer）并显示其数值於萤幕上。 小提示： ">>"是MATLAB的提示符号（Prompt），但在PC中文视窗系统下，由於编码方式不同，此提示符号常会消失不见，但这并不会影响到MATLAB的运算结果。 我们也可将上述运算式的结果设定给另一个变数x： x = (5*2+1.3-0.8)*10^2/25 x = 42 此时MATLAB会直接显示x的值。由上例可知，MATLAB认识所有一般常用到的加（+）、减（-）、乘（*）、除（/）的数学运算符号，以及幂次运算（^）。 小提示： MATLAB将所有变数均存成double的形式，所以不需经过变数宣告（Variable declaration）。MATLAB同时也会自动进行记忆体的使用和回收，而不必像C语言,必须由使用者一一指定.这些功能使的MATLAB易学易用，使用者可专心致力於撰写程式，而不必被软体枝节问题所干扰。 若不想让MATLAB每次都显示运算结果，只需在运算式最後加上分号（；）即可，如下例： y = sin(10)*exp(-0.3*4^2);

若要显示变数y的值，直接键入y即可： >>y y =-0.0045 在上例中，sin是正弦函数，exp是指数函数，这些都是MATLAB常用到的数学函数。 下表即为MATLAB常用的基本数学函数及三角函数： 小整理：MATLAB常用的基本数学函数 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数z的相角(Phase angle) sqrt(x)：开平方 real(z)：复数z的实部 imag(z)：复数z的虚部 conj(z)：复数z的共轭复数 round(x)：四舍五入至最近整数 fix(x)：无论正负，舍去小数至最近整数 floor(x)：地板函数，即舍去正小数至最近整数 ceil(x)：天花板函数，即加入正小数至最近整数 rat(x)：将实数x化为分数表示 rats(x)：将实数x化为多项分数展开

实验一数据处理方法MATLAB实现

实验一数据处理方法的MATLAB实现 一、实验目的 学会在MATLAB环境下对已知的数据进行处理。 二、实验方法 1. 求取数据的最大值或最小值。 2. 求取向量的均值、标准方差和中间值。 3．在MATLAB环境下，对已知的数据分别进行曲线拟合和插值。 三、实验设备 1．586以上微机，16M以上内存，400M硬盘空间，2X CD-ROM 2．MATLAB5.3以上含CONTROL SYSTEM TOOLBOX。 四、实验内容 1．在MATLAB环境下，利用MATLAB控制系统工具箱中的函数直接求取数据的最大值或最小值，以及向量的均值、标准方差和中间值。 2．在MATLAB环境下，选择合适的曲线拟合和插值方法，编写程序，对已知的数据分别进行曲线拟合和插值。 五、实验步骤 1. 在MATLAB环境下，将已知的数据存到数据文件mydat.mat中。 双击打开Matlab，在命令窗口（command window）中，输入一组数据：实验一数据处理方法的MATLAB实现 一、实验目的 学会在MATLAB环境下对已知的数据进行处理。 二、实验方法 1. 求取数据的最大值或最小值。 2. 求取向量的均值、标准方差和中间值。 3．在MATLAB环境下，对已知的数据分别进行曲线拟合和插值。 三、实验设备 1．586以上微机，16M以上内存，400M硬盘空间，2X CD-ROM 2．MATLAB5.3以上含CONTROL SYSTEM TOOLBOX。 四、实验内容

1．在MATLAB环境下，利用MATLAB控制系统工具箱中的函数直接求取数据的最大值或最小值，以及向量的均值、标准方差和中间值。 2．在MATLAB环境下，选择合适的曲线拟合和插值方法，编写程序，对已知的数据分别进行曲线拟合和插值。 五、实验步骤 1. 在MATLAB环境下，将已知的数据存到数据文件mydat.mat中。 双击打开Matlab，在命令窗口（command window）中，输入一组数据： x=[1,4,2,81,23,45] x = 1 4 2 81 2 3 45 单击保存按钮，保存在Matlab指定目录（C:\Program Files\MATLAB71）下，文件名为“mydat.mat”。 2. 在MATLAB环境下，利用MATLAB控制系统工具箱中的函数直接求取数据的最大值或最小值，以及向量的均值、标准方差和中间值。 继续在命令窗口中输入命令： （1）求取最大值“max（a）”； >> max(x) ans = 81 （2）求取最小值“min（a）”； >> min(x) ans = 1 （3）求取均值“mean（a）”； >> mean(x) ans =

matlab学习笔记(入门)

数据类：double，unit8，unit16，unit32，int8，int16，int32，single，char，logical！Matlab中所有数值计算都可以用double类来进行！，unit8实际中最常用的图像 图像类型：亮度图像，二值图像，索引图像，RGB图像 亮度图像：是数据矩阵，若是unit8或uint16则是【0,255】或者是【0，65535】，若是double 类，则像素取值是浮点数 二值图像只有：0和1的逻辑数组！ 、 简单操作： 读图并显示详细情况 >> f=imread('E:\image\book.pgm');whos Name Size Bytes Class Attributes f 289x338 97682 uint8 将图像垂直翻转： >> f=imread('E:\image\book.pgm');fp=f(end:-1:1, : );imshow(fp) 将图像上下左右翻转： f=imread('E:\image\book.pgm');fc=f(end:-1:1,end:-1:1);imshow(fc) 将图像二次采样并显示详情： >> fs=f(1:2:end,1:2:end);imshow(fs) >> whos fs Name Size Bytes Class Attributes fs 145x169 24505 uint8 将图像取出一部分： >> fg=f(200:250,200:300);imshow(fg) 显示图像中的一条水平扫描线： >> plot(f(200, : ) 将两幅图像进行相乘： f=imread('c:\image\liangdian.jpg');g=imread('c:\image\shuiguo.jpg'); g=g(300:715,500:1149);f=f(1:416,1:650);f d=double(f);gd=double(g); p=fd.*gd;数组乘！ pmax=max(p(:));pmin=min(p(:));取最大最小值！ pn=mat2gray(p);figure,imshow(pn) 亮度变化： 函数imadjust是对灰度图像进行亮度变化的基本ipt工具： g=imadjust（f，[low-in high-in]，[low-in high-in]，gamma） Gamma为1线性映射，大于1，则映射被加权至更低（更暗的）输出值，小于一，加权至更高的输出值 明暗反转图像（负片）参数不同： >> f=imread('E:\image\book.pgm');g=imadjust(f, [0 1],[1 0 ]);imshow(g) >> f=imread('E:\image\book.pgm');g=imadjust(f, [0 1],[1 0 ],2);imshow(g) >> f=imread('E:\image\book.pgm');g=imadjust(f, [0 1],[1 0 ],0.5);imshow(g) 另外也可以这样：进行明暗反转： g=imcomplement(f);imshow(g) 将0.5到0.75之间的灰度级拓展到0-1，可用于突出我们感兴趣的亮度带

Matlab2012教程--经典教程

第1章基础准备及入门 1.1 最简单的计算器使用法 为易于学习，本节以算例方式叙述，并通过算例归纳一些MATLAB最基本的规则和语 法结构。建议读者，在深入学习之前，先读一读本节。 2 【例1.3-1】求[122(74)]3的算术运算结果。本例演示：最初步的指令输入形式 和必需的操作步骤。 （1）用键盘在MA TLAB指令窗中输入以下内容 >> (12+2*(7-4))/3^2 （2）在上述表达式输入完成后，按[Enter] 键，该指令被执行，并显示如下结果。 ans = 2 〖说明〗 本例在指令窗中实际运行的情况参见图 1.3-1。 指令行“头首”的“>>”是“指令输入提示符”，它是自动生成的。本书在此后的输入指令前将不再带提示符“>>”。理由是：（A）为使本书简洁；（B）本书用MATLAB 的M-book写成，而在M-book中运行的指令前是没有提示符的。 5

MATLAB的运算符（如+、- 等）都是各种计算程序中常见的习惯符号。 一条指令输入结束后，必须按[Enter] 键，那指令才被执行。 由于本例输入指令是“不含赋值号的表达式”，所以计算结果被赋给MATLAB的一个默认变量“ans”。它是英文“answer”的缩写。 【例1.3-2】“续行输入”法。本例演示：或由于指令太长，或出于某种需要，输入指令行必 须多行书写时，该如何处理。 S=1-1/2+1/3-1/4+ ... 1/5-1/6+1/7-1/8 S = 0.6345 〖说明〗 MA TLAB用3个或3个以上的连续黑点表示“续行”，即表示下一行是上一行的继续。 本例指令中包含“赋值号”，因此表达式的计算结果被赋给了变量S。 指令执行后，变量S被保存在MA TLAB 的工作空间（Workspace）中，以备后用。如果用户不用clear 指令清除它，或对它重新赋值，那么该变量会一直保存在工作空间中， 直到本MATLAB 指令窗被关闭为止。 1.3.3数值、变量和表达式 前节算例只是表演了“计算器”功能，那仅是MA TLAB全部功能中小小一角。为深入 学习MA TLAB，有必要系统介绍一些基本规定。本节先介绍关于变量的若干规定。 一数值的记述 MATLAB的数值采用习惯的十进制表示，可以带小数点或负号。以下记述都合法。 3 -99 0.001 9.456 1.3e-3 4.5e33 在采用IEEE浮点算法的计算机上，数值通常采用“占用64位内存的双精度”表示。 其相对精度是eps （MATLAB的一个预定义变量），大约保持有效数字16位。数值范围大308308 致从10到10。 二变量命名规则 变量名、函数名是对字母大小写敏感的。如变量myvar和MyVar表示两个不同的变量。 sin是MATLAB定义的正弦函数名，但SIN,Sin等都不是。 变量名的第一个字符必须是英文字母，最多可包含63个字符（英文、数字和下连符）。 如myvar201是合法的变量名。 变量名中不得包含空格、标点、运算符，但可以包含下连符。如变量名my_var_201是合法的，且读起来更方便。而my,var201由于逗号的分隔，表示的就不是一个变量名。 6

Matlab经典案例

1、三维曲线 >> t=0:pi/50:10*pi; >> plot3(sin(2*t),cos(2*t),t) >> axis square >> grid on 2、一窗口多图形 >> t=-2*pi:0.01:2*pi; >> subplot(3,2,1) >> plot(t,sin(t)) >> subplot(3,2,2) >> plot(t,cos(t)) >> subplot(3,2,3) >> plot(t,tan(t)) >> axis([-pi pi -100 100]) >> subplot(3,2,4) >> plot(t,cot(t)) >> axis([-pi pi -100 100]) >> subplot(3,2,5) >> plot(t,atan(t)) >> subplot(3,2,6) >> plot(t,acot(t)) 3、图形样式、标注、题字(也可以利用菜单直接 Insert) >> x=0:pi/20:2*pi; >> plot(x,sin(x),'b-.') >> hold on >> plot(x,cos(x),'r--') >> hold on >> plot(x,sin(x)-1,'g:')

>> hold on >> plot(x,cos(x)-1) >> xlabel('x'); >> xlabel('x轴'); >> ylabel('y轴'); >> title('图形样式、标注等'); >> text(pi,sin(pi),'x=\pi'); >> legend('sin(x)','cos(x)','sin(x)-1','cos(x)-1'); >> [x1,y1]=ginput(1) %利用鼠标定位查找线上某点的值x1 = 2.0893 y1 = -0.5000 >> gtext('x=2.5') %鼠标定位放置所需的值在线上 4、 >> fplot('[sin(x),cos(x),sqrt(x)-1]',[0 2*pi]) M文件：myfun.m 内容如下： function y=myfun(x) y(:,1)=sin(x); y(:,2)=cos(x); y(:,3)=x^(1/2)-1; 再运行：>> fplot('myfun',[0 2*pi]) 同样可以得到右图 5、 >> [x,y]=fplot('sin',[0 2*pi]); >> [x1,y1]=fplot('cos',[0 2*pi]); >> plot(x,y,'-r',x1,y1,'-.k') >> legend('y=sinx','y=cosx') 6、

matlab经典代码大全

哈哈哈 MATLAB 显示正炫余炫图：plot(x,y1,'* r',x,y2,'o b') 定义【0,2π】；t=0:pi/10:2*pi; 定义函数文件：function [返回变量列表]=函数名（输入变量列表） 顺序结构：选择结构 1）if-else-end语句 其格式为： if 逻辑表达式 程序模块1； else 程序模块2； End 图片读取：%选择图片路径 [filename, pathname] = ... uigetfile({'*.jpg';'*.bmp';'*.gif'},'选择图片'); %合成路径+文件名 str=[pathname,filename]; %为什么pathname和filename要前面出现的位置相反才能运行呢？？？%读取图片 im=imread(str); %使用图片 axes(handles.axes1); %显示图片 imshow(im); 边缘检测： global im str=get(hObject,'string'); axes (handles.axes1); switch str case ' 原图' imshow(im); case 'sobel' BW = edge(rgb2gray(im),'sobel'); imshow(BW); case 'prewitt' BW = edge(rgb2gray(im),'prewitt');

imshow(BW); case 'canny' BW = edge(rgb2gray(im),'canny'); imshow(BW);Canny算子边缘定位精确性和抗噪声能力效果较好，是一个折中方案 end; 开闭运算： se=[1,1,1;1,1,1;1,1,1;1,1,1]; %Structuring Element I=rgb2gray(im); imshow(I,[]);title('Original Image'); I=double(I); [im_height,im_width]=size(I); [se_height,se_width]=size(se); halfheight=floor(se_height/2); halfwidth=floor(se_width/2); [se_origin]=floor((size(se)+1)/2); image_dilation=padarray(I,se_origin,0,'both'); %Image to be used for dilation image_erosion=padarray(I,se_origin,256,'both'); %Image to be used for erosion %%%%%%%%%%%%%%%%%% %%% Dilation %%% %%%%%%%%%%%%%%%%%% for k=se_origin(1)+1:im_height+se_origin(1) for kk=se_origin(2)+1:im_width+se_origin(2) dilated_image(k-se_origin(1),kk-se_origin(2))=max(max(se+image_dilation(k-se_origin(1):k+halfh eight-1,kk-se_origin(2):kk+halfwidth-1))); end end figure;imshow(dilated_image,[]);title('Image after Dilation'); %%%%%%%%%%%%%%%%% %%% Erosion %%% %%%%%%%%%%%%%%%%% se=se'; for k=se_origin(2)+1:im_height+se_origin(2) for kk=se_origin(1)+1:im_width+se_origin(1) eroded_image(k-se_origin(2),kk-se_origin(1))=min(min(image_erosion(k-se_origin(2):k+halfwidth -1,kk-se_origin(1):kk+halfheight-1)-se)); end end figure;imshow(eroded_image,[]);title('Image after Erosion'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% Opening(Erosion first, then Dilation) %%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

MATLAB在化工数据处理中的应用

化工数据处理与实验设计 期中论文 Matlab在化工数据处理中的应用 姓名：陈奕含 学号：2012115134 班级：化学工程与工艺

Matlab在化工数据处理中的应用 一、MATLAB简介 Matlab软件包最初是1967年由Clere Maler用FORTRAN语言设计和编写的。1984年Mathworks公司用C语言完成了Matlab的商业化版本并推向市场。经过20余年的改进，Matlab已发展成为一个具有极高通用性的、带有众多实用工具的运算平台，成为国际上广泛认可的优秀科学计算软件。 Matlab是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国Mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 Matlab和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。Matlab可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 Matlab的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用Matlab来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且Matlab也吸收了像Maple等软件的优点，使Matlab成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JA V A的支持。其特点如下： 1.支持多平台操作系统（Windows、Unix等）。 2.是一种简单易学的编程语言。 3.Matlab程序很容易维护。 4.编程效率很高。由于用户程序可直接调用大量的Matlab函数，因此编程速度快。 5.用途广泛。可用于数值计算和符号计算、数据分析、工程与科学绘图、图形用户界面设计、建模和仿真、控制系统设计、数字图像信号处理以及财务工程等。 6.功能超强。包含600多个常用算法内建函数，有众多面向具体应用的工具箱（如偏微分方程、最优化方法、数理统计、样条函数、神经网络工具箱等）和simulink仿真模块。此外，其他产品延伸了Matlab的能力，包括数据采集和依靠Matlab语言编程产生独立的C/C++代码等等。其算法函数大多由国际知名专家完成，算法稳定可靠、效率高。 7.具有开放式结构，扩展功能强。Matlab的开放式结构使Matlab产品族很

基于MATLAB的EXCEL数据计算与分析

基于MATLAB的EXCEL数据计算与分析 潜刘方 摘要：再怎么样希望先看摘要，阅读本文需要一定的MA TLAB基础知识，不需要excel相关知识。结合本人近期工作上的需要测量计算，想偷懒就选择了利用MATLAB偷懒，于是便有了本文。本文首先利用MA TLAB读取数据，计算，将数据写入excel，然后花了很大的精力来根据实际需要画图，最后将图保存在excel所在的文件夹下。这个m文件可谓花了我不少的时间和精力。最后根据m文件的不足（不能将图形输入到excel文档当中），进一步弥补这不足，就有了exlink（也叫excel link），在网上搜索了相关的知识，发现很多关于exlink 的培训，觉得实在可笑，所以就将exlink的使用写的比较详细，以供读者自行分析体会。关键字：MATLAB excel exlink 接口 一、前沿 MATLAB是一款应用在各个领域的数学软件，最初叫做矩阵实验室，专用于矩阵的运算，后来的版本再各个领域都得到了很好的应用，比如：通信、电力电子、电机控制、运动控制、计算机控制、自动控制，DSP数字信号处理。但是MATLAB对于数据的处理与可视化是很多软件所不能及的。 EXCEL作为办公必备软件，能对简单数据分析计算与作图分析，但是处理复杂数据显得力不从心，比如三维作图就无法利用EXCEL作出；EXCEL本身的函数远远没有MATLAB 多，MATLAB作为数据有其独特的优势，集成了很多数学函数，包括数据拟合差值等。MATLAB 可以从EXCEL中读取数据，经过相关运算之后又可以将数据写入EXCEL，假如需要重复性的对excel可以利用MATLAB编写函数，每次只要运行MATLAB程序就可以完成，大大节省时间和精力。 另外，MATLAB还有与EXCEL的接口，叫做EXLINK，运用这个接口可以在excel中完成MATLAB函数的调用，还能传送数据给MATLAB，从MATLAB当中读取数据，从MATLAB 当中读取图形，使用方便，操作简单。 二、基于MATLAB的数据分析 数据分析操作流程主要分为三步：第一步，从excel中读取数据；第二部：利用MATLAB 大量函数对数据分析处理；第三步：将分析结果写入excel中。在整个过程中，不需要打开excel软件，操作十分方便，每次操作唯一要做就是修改excel所在的目录及文件名。主要函数如下（具体使用方法可在MATLAB命令窗口输入help +函数名查看）：Xlsread 从excel中读数据 Xlswrite 向excel中邪数据 num2str 将数字转换为字符串 strncmp 字符串比较 polyfit 数据拟合 polyval 具体数值代入求值 plot 作图

基于MATLAB的数据实时采集与处理的实现_梁湘

０．引言 ＭＡＴＬＡＢ／ＳＩＭＵＬＩＮＫ是现在流行的仿真软件。ＭＡＴＬＡＢ集数学计算结果可视化和编程于一体，能够方便地进行科学计算和大量工程运算的数学软件；ＳＩＭＵＬＩＮＫ是ＭＡＴＬＡＢ的常用组件，它是基于ＭＡＴＬＡＢ的语言环境下实现动态装置建模，仿真和分析的一个集成环境，支持连续、离散及两者混合的线性和非线性装置，也支持具有多种采样速率的多速率装置，被广泛的用于控制系统设计和系统仿真等诸多领域。但是ＭＡＴＬＡＢ不能直接对硬件进行读写操作，从而影响了在控制系统仿真上的应用范围。ＭＡＴＬＡＢ提供了众多外部函数接口，本文从中选择ＭＥＸ文件接口作为ＭＡＴＬＡＢ对外界进行读写的通道，首先对ＵＳＢＣＡＮ接口卡进行设计，接下来以这个接口卡为通道，实现了ＭＡＴＬＡＢ对于硬件的访问。 １．ＵＳＢ－ＣＡＮ转接卡的设计 首先对于ＵＳＢＣＡＮ转接卡进行设计，以便于设计出的转接卡能作为通信通道，让ＭＡＴＬＡＢ能够通过ＵＳＢＣＡＮ转换模块采集ＣＡＮ总线上的数据，然后进行仿真。 在基于ＣＡＮ总线的控制系统中，作为下层网络的ＣＡＮ总线与计算机之间的通信以往是通过基于ＲＳ２３２接口、ＰＣ机上的ＩＳＡ总线和ＰＣＩ总线的通信适配卡来完成的。这些种类的转接卡有传输速率较慢，设计复杂和不便于扩展等缺点。 本课题设计的ＣＡＮ－ＵＳＢ转接卡的功能是：采集ＣＡＮ总线上的数据，通过ＵＳＢ总线和上位控制站之间进行数据传输。在通讯方面，一方面要满足ＣＡＮ的协议标准，实现开放性，互操作性，在高速通信的情况下做到拥有较好的抗干扰性能；另一方面，转接卡兼容ＵＳＢ１．１总线，转接卡可以通过ＵＳＢ连接到ＰＣ机。其主要设计参数如下：数据传输速率：波特率在５Ｋｂｉｔ／ｓ～１Ｍｂｉｔ／ｓ范围内： ＵＳＢ总线标准：满足ＵＳＢ１．１协议的标准ＵＳＢ设备Ａ／Ｂ插座 ＣＡＮ总线接口：ＤＢ９针型插座，符合ＤｅｖｉｃｅＮＥＴ和ＣＡＮｏｐｅｎ标准ＣＡＮ协议：支持ＣＡＮ２．０Ｂ协议（兼容ＣＡＮ２．０Ａ协议） 供电方式：ＵＳＢ总线供电（＋５Ｖ）或者使用外接电源（＋９Ｖ～＋２５Ｖ，４００ｍＡ） 运行环境：ＷｉｎｄｏｗｓＸＰ 转接卡的方案设计的思想是在已经存在的ＣＡＮ总线网络中增加一个数据采集节点，用来采集各个节点发来的的数据，再通过ＵＳＢ总线上传到ＰＣ机，进行分析过和存储等操作。设计的具体实现是采用单片机通过控制ＵＳＢ控制芯片，与ＰＣ机进行通信，从ＰＣ机得到采集数据的指令和向ＰＣ机上传数据。同时，因为ＵＳＢ总线速度较快，并且是要从ＣＡＮ网络读取数据分析，所以单片机通过不同的片选信号控制两个ＣＡＮ控制芯片，分别能够从两路ＣＡＮ网络上采集数据，和将从ＰＣ机得到的数据通过两个ＣＡＮ控制芯片分别发送出去。所设计的系统结构如下图１所表示： 两路ＣＡＮ ＣＡＮ总线 ＵＳＢ协议具有１．１和２．０两个版本。从ＭＡＴＬＡＢ对数据进行读写出发，并且处于ＣＡＮ总线所能提供的速率限制，ＵＳＢ协议１．１版本已经能满足我们所需要的性能要求。ＵＳＢ１．１提供两种数据传输速率，低速传输为１．５Ｍｂｐｓ，全速传输为１２Ｍｂｐｓ，并支持所有ＵＳＢ的特性，如热插拔、具有统一的设备标准以及可以连接多个设备等。ＣＡＮ总线目前有ＣＡＮ协议２．０Ａ和ＣＡＮ协议２．０Ｂ两个版本。本次设计所选用的单片机和ＣＡＮ网络控制芯片均可支持ＣＡＮ协议２．０Ａ和ＣＡＮ协议２．０Ｂ两个版本。 图１ＵＳＢＣＡＮ转接卡的总体设计 本设计选用８９ＣＳ５２作为单片机来对ＣＡＮ控制芯片和ＵＳＢ控制芯片操作。８９ＣＳ５２是ＡＴＭＥＬ公司生产的低电压，高性能的单片机，兼容８０Ｃ５１构架。选择ＳＪＡ１０００做ＣＡＮ总线控制芯片。这款芯片是一个由飞利浦公司生产的独立的ＣＡＮ控制器，它在汽车和普通的工业应用上都具有较为先进的特征。它能适合多种应用，特别是在系统优化、诊断和维护方面，并且能和支持８０Ｃ５１构架的单片机兼容。选择ＰＣＡ８２Ｃ２５０作为ＣＡＮ总线接口芯片，这款芯片专为ＣＡＮ总线进行差分通信而设计。 根据设计要求，ＵＳＢＣＡＮ转接卡应该直接能与ＰＣ机相连接，再与ＣＡＮ总线进行通信。对比ＵＳＢ协议三层设备：主机（ＨＯＳＴ），集线器（ＨＵＢ）和设备（ＮＯＤＥ），ＰＣ机作为主机，ＵＳＢＣＡＮ转接卡是作为最下面一个级别：设备（ＮＯＤＥ），只需要与其相连而不需要再作为ＨＵＢ接入其它ＵＳＢ设备。所以需要挑选的是能够在充当设备（ＮＯＤＥ）的芯片，这款芯片应该可以与５１系列单片机相连通，并且支持ＵＳＢ协议１．１版本。市场上此类控制芯片主要有两种，一种是将单片机（ＭＣＵ）集成到ＵＳＢ芯片上，代表如ＥＺ－ＵＳＢＦＸ２；另一种是纯粹的ＵＳＢ接口芯片，如ＰＤＩＵＳＢＤ１２。考虑到前文提出的设计指标等因素，本文选择了飞利浦公司的ＰＤＩＵＳＢＤ１２芯片。ＰＤＩＵＳＢＤ１２是一款性价比很高的ＵＳＢ器件，广泛的采用与许多ＰＣ机的外设。 硬件设计完成后，对于固化在单片机ＲＯＭ中的程序进行设计。单片机ＲＯＭ中的程序主要采用中断驱动。它将ＵＳＢ总线和ＣＡＮ总线上的数据请求都当作中断处理，在相应的中断处理程序中再将所采集到的数据转发到另外一个总线。 最后对于转接卡的驱动程序进行设计。驱动程序采用Ｗｉｎｄｏｗｓ体系下ＷＤＭ驱动程序模型。在这种模型下，对于ＵＳＢ设备来说，驱动程序可以分为ＵＳＢ底层驱动程序和ＵＳＢ功能驱动程序。ＵＳＢ底层驱动程序在实际运行中对实际硬件进行操作，实现了复杂的底层通信；ＵＳＢ功能驱动程序则一般由设备开发者编写，逻辑位置位于ＵＳＢ底 基于ＭＡＴＬＡＢ的数据实时采集与处理的实现 梁湘 （同济大学机械工程学院中国上海２０００９２） 【摘要】本文的利用ＵＳＢ与ＣＡＮ两种技术的优点设计了通信转接卡，提高整个网络的通信质量。接下来，采用该转接卡作为ＭＡＴＬＡＢ与ＳＩＭＵＬＩＮＫ环境下实时仿真的通信转接卡，进行ＭＡＴＬＡＢ环境下对ＣＡＮ总线数据的实时采集、处理与仿真功能，完成ＭＡＴＬＡＢ与ＳＩＭＵＬＩＮＫ下对控制系统的仿真和实时控制，为类似的通信转接卡的研究提供一定的参考。 【关键词】ＵＳＢ－ＣＡＮ转接卡；半实物平台；ＭＡＴＬＡＢ；ＳＩＭＵＬＩＮＫ ＴｈｅＲｅａｌｉｚａｔｉｏｎｏｆＲｅａｌ－ｔｉｍｅＲｅａｄｉｎｇａｎｄＷｒｉｔｉｎｇｉｎＭＡＴＬＡＢ ＬｉａｎｇＸｉａｎｇ （ＳｃｈｏｏｌｏｆＭｅｃｈａｎｉｃａｌＥｎｇｉｎｅｅｒｉｎｇ，ＴｏｎｇｊｉＵｎｉｖｅｒｓｉｔｙ，Ｓｈａｎｇｈａｉ２０００９２） 【Ａｂｓｔｒａｃｔ】ＴｈｅｐａｐｅｒｔａｋｅｓａｄｖａｎｔａｇｅｏｆＣＡＮ－ｂｕｓａｎｄＵＳＢ－ｂｕｓｔｏｄｅｓｉｇｎＵＳＢＣＡＮａｄａｐｔｅｒ．Ｔｈｅａｄａｐｔｅｒｉｍｐｒｏｖｅｓｔｈｅｑｕａｌｉｔｙｏｆｎｅｔｗｏｒｋ’ｓｃｏｍｍｕｎｉｃａｔｉｏｎ．Ｆｏｌｌｏｗｉｎｇ，ｔｈｅｐａｐｅｒｔａｋｅｓＵＳＢＣＡＮａｄａｐｔｅｒａｓａｃｈａｎｎｅｌｗｉｔｈｗｈｉｃｈＭＡＴＬＡＢｃａｎｃｏｌｌｅｃｔｄａｔａｆｒｏｍＣＡＮ－ｂｕｓ．Ｉｎｔｈｉｓｗａｙ，ＭＡＴＬＡＢｃａｎｃｏｌｌｅｃｔ，ｈａｎｄｌｅａｎｄｓｉｍｕｌａｔｅｄａｔａｉｎｒｅａｌ ｔｉｍｅ．