matlab代码大全

1、超效率 DEA clearX= ;%用户输入多指标输入矩阵 XY= ;%用户输入多指标输入矩阵 Yn=size(X, 1);m=size(X,1);s=size(Y,1);epsilon=10-10;%定义非阿基米德无穷小量 =10-10f=zeros(1,n) -epsilon*ones(1,m+s) 1;A=zeros(1,n+m+s+1);b=0LB=zeros(n+m+s+1,1);UB=;LB(n+m+s+1)=-inf;for i=1:n;Aeq= X(:,1:i-1),zeros(m,1),X(:,i+1:n) eye(m)zeros(m,s)-X(:,i)Y(:,1:i-1),zeros(s,1),Y(:,i+1:n) zeros(s,m) -eye(s) zeros(s,1)beq=zeros(m,1)Y(:,i) ;w(:,i)= linprog(f,A,b,Aeq,beq,LB,UB);%解线性规划, 得 DMUi 的最佳权向量 wiendw % 输 出 最 佳。

2、HTML 代码教程教程一、基本标志1.2.3.4.Html 是英文 HyperText Markup Language 的缩写，中文意思是“超文本标志语言” ，用它编写的文件(文档) 的扩展名是.html 或.htm ，它们是可供浏览器解释浏览的文件格式。您可以使用记事本、写字板或 FrontPage Editor 等编辑工具来编写 Html 文件。Html 语言使用标志对的方法编写文件，既简单又方便，它通常使用来表示标志的开始和结束(例如标志对)，因此在 Html 文档中这样的标志对都必须是成对使用的。在此教程中，我先讲一下 Html 的基本标志：1.标志用于 Html 文档的最前边，用来标识 Html 文。

3、clear ;close all;%Step1 获取图像 装入待处理彩色图像并显示原始图像Scolor = imread(C:UserslenovoDesktop程序chepail2.jpg);%将彩色图像转换为黑白并显示Sgray = rgb2gray(Scolor);figure,imshow(Scolor),title(原始彩色图像);figure,imshow(Sgray),title(原始灰度图像);%bw3=medfilt2(Sgray,5 5);%figure,imshow(bw3);title(图像中值滤波);%bw4=imfilter(bw3,fspecial(laplacian);%figure,imshow(bw4);title(图像中值滤波);%Egray=imsubtract(bw3,bw4);%figure,imshow(Egray);%Step2 图像预处理 对 Sgray 原始黑白图像进行开操作得到。

4、1.色彩空间转换function r,g=rgb_RGB(Ori_Face)R=Ori_Face(:,:,1);G=Ori_Face(:,:,2);B=Ori_Face(:,:,3);R1=im2double(R); % 将 uint8 型转换成 double 型G1=im2double(G);B1=im2double(B);RGB=R1+G1+B1;row=size(Ori_Face,1); % 行像素column=size(Ori_Face,2); % 列像素for i=1:rowfor j=1:columnrr(i,j)=R1(i,j)/RGB(i,j);gg(i,j)=G1(i,j)/RGB(i,j);endendrrr=mean(rr);r=mean(rrr);ggg=mean(gg);g=mean(ggg);2.均值和协方差t1=imread(D:matlab皮肤库 1.jpg);r1,g1=rgb_RGB(t1);t2=imread(D:matlab皮肤库 2.jpg);r2,g2=rgb_RGB。

5、%function C=DES(M,K)% -% 数据加密标准 DES 算法的 Matlab 实现% 输入明文 M 和密钥 K 用 16 进制数表示，字符串格式%-% -第一步：数据初始化-%以字符串形式输入的明文转换成 01 数字串M=123456ABCD132536;MB=;for i=1:16Mi=M(i);MBi=0000,dec2bin(hex2dec(Mi);MBi=MBi(end-3:end);MBi=str2num(MBi(1),str2num(MBi(2),str2num(MBi(3),str2num(MBi(4);MB=MB,MBi;endM=。

6、对于正弦序列 进行 16 点的 dft 仿真)(8sin)(3RxNk=16; n1=0:1:15; xa1=sin(2*pi*n1/k); subplot(2,2,1) plot(n1,xa1) xlabel(t/T);ylabel(x(n); xk1=fft(xa1);xk1=abs(xk1); subplot(2,2,2) stem(n1,xk1) xlabel(k);ylabel(X(k); n2=0:1:7; xa2=sin(2*pi*n2/k); subplot(2,2,3) plot(n2,xa2) xlabel(t/T);ylabel(x(n); xk2=fft(xa2);xk2=abs(xk2); subplot(2,2,4) stem(n2,xk2) xlabel(k);ylabel(X(k);function y=myfft(xr,n) p=0:n-1 ; % 开始倒位序。

7、采用 Newton 法求方程 的根的代码13xffunction NewtonMethodx=1.0;it=0;table=;while 1f,df=fdfCal(x);table=table;it,x,f,dfx=x-f/df;f2,df2=fdfCal(x);it=it+1;if f21e-10break;end; end;return;function f,df=fdfCal(x)f=x3-x-1;df=3*x2-1;return;。

8、说明：本代码为2048游戏matlab代码，程序未进行游戏结束判定，节目如下。function g2048(action)global totalscore flag score_boardif nargin0score_board(num_0(ceil(l*rand)=2+2*(rand0.1);end endendfunction Y=addnum(X)num_0=find(X=0);l=length(num_0);X(num_0(ceil(l*rand)=2+2*(rand0.1);Y=X;endfunction B score=move(A)score=0;for k=1:2for i=1:3if A(i)=0for j=i:3A(j)=A(j+1);endA(4)=0;endendendif A(1)=A(2)if A(3)=A(4)A(1)=A(1)+A(2);A(2)=A(3)+A(4);A(3)=0;A(4)=0;score=A(1)+A(2);else A(1)=A(1)+A(2);A(2)=A(3);A(3)。

9、遗传算法function Zp,Y1p,Y2p,Y3p,Xp,LC1,LC2=JSPGA(M,N,Pm,T,P)%-% JSPGA.m% 车间作业调度问题遗传算法%-% 输入参数列表% M 遗传进化迭代次数% N 种群规模(取偶数)% Pm 变异概率% T mn 的矩阵，存储 m 个工件 n 个工序的加工时间% P 1n 的向量，n 个工序中，每一个工序所具有的机床数目% 输出参数列表% Zp 最优的 Makespan 值% Y1p 最优方案中，各工件各工序的开始时刻，可根据它绘出甘特图% Y2p 最优方案中，各工件各工序的结束时刻，可根据它绘出甘特图% Y3p 最优方案中，各工件各工序使用的机器编号% Xp 最优决策变量的。

10、%傅里叶变换clc;clear all;close all;ticFs=128;%采样频率，频谱图的最大频率T=1/Fs;%采样时间，原始信号的时间间隔L=256;%原始信号的长度，即原始离散信号的点数t=(0:L-1)*T;%原始信号的时间取值范围x=7*cos(2*pi*15*t-pi)+3*cos(2*pi*40*t-90*pi/180)+3*cos(2*pi*30*t-90*pi/180);z=7*cos(2*pi*15*t-pi)+3*cos(2*pi*40*t-90*pi/180);z1=6*cos(2*pi*30*t-90*pi/180);z1(1:L/2)=0;z=z+z1;y=x;%+randn(size(t);figure;plot(t,y)title(含噪信号)xlabel(时间（ s）)hold onplot(t,z,r-)N=2nextpow2(L);%N 为使 2N=L 的最小幂Y=fft(y,N)/N*2;Z。

11、注：显示图像后，单击鼠标左键连出完全背景区域，双击左键结束；图像改变后，再单击鼠标左键连出完全前景区域，双击左键结束。即可得到结果。A=imread(ppmm2.bmp);imshow(A);hold on;Threshold = 1;UniBack=00255;% choose background object, and B results into the Foreground Object and the% Boundary Region% choose the outer line of the boundary region. Single left click the mouse to specify vertice.Double left click or single right click to finish;OuterBW = roipoly(A);B(:,:,1) = immultiply(A(:,:,1),OuterBW);B(:。

12、 8 第 2章 图像 获取 2.3.2 二维连续傅里叶变换 例 2.2 figure(1); %建立图形窗口 1 u,v = meshgrid(-1:0.01:1); %生成二维频域网格 F1 = abs(sinc(u.*pi); F2 = abs(sinc(v.*pi); F=F1.*F2; %计算幅度频谱 F=|F(u,v)| surf(u,v,F); %显示幅度频谱，如图 2.3(b) shading interp; %平滑三维曲面上的小格 axis off; %关闭坐标系 figure(2); %建立图形窗口 2 F1=histeq(F); %扩展 F 的对比度以增强视觉效果 imshow(F1); %用图像来显示幅度频谱，如图 2.3(c) 第 3章 图像变换 3.4.4 二维 FFT 的 MATLAB 实现 例 3.2 简单图像及其傅里叶变换 MA。

13、1.图像反转MATLAB程序实现如下：I=imread(xian.bmp);J=double(I);J=-J+(256-1); %图像反转线性变换H=uint8(J);subplot(1,2,1),imshow(I);subplot(1,2,2),imshow(H);2.灰度线性变换MATLAB程序实现如下：I=imread(xian.bmp);subplot(2,2,1),imshow(I);title(原始图像);axis(50,250,50,200);axis on; %显示坐标系I1=rgb2gray(I);subplot(2,2,2),imshow(I1);title(灰度图像);axis(50,250,50,200);axis on; %显示坐标系J=imadjust(I1,0.1 0.5,); %局部拉伸，把0.1 0.5内的灰度拉伸为0 1subplot(2,2,3),imshow(J);title(线性变换图。

14、例 Error No text of specified style in document 1 周期信号 方波 的展开 fb jinshi m close all clear all N 100 取展开式的项数为2N 1项 T 1 fs 1 T N sample 128 为了画出波形 设置每个周期的采样点数 dt T N sample t 0 dt 10 T dt n N N Fn sinc n。

15、附 录Part1：对 cameraman 原始图像处理的仿真程序：clcclear allclose allA = imread(cameraman.bmp); % 读入图像subplot(2,4,1);imshow(A);title(原图);x_mask = 1 0;0 -1; % 建立 X 方向的模板y_mask = rot90(x_mask); % 建立 Y 方向的模板I = im2double(A); % 将图像数据转化为双精度dx = imfilter(I, x_mask); % 计算 X 方向的梯度分量dy = imfilter(I, y_mask); % 计算 Y 方向的梯度分量grad = sqrt(dx.*dx + dy.*dy); % 计算梯度grad = mat2gray(grad); % 将梯度矩阵转换为灰度图像level = graythresh(grad); % 计算灰度阈值BW = i。

16、abs 绝对值、模、字符的 ASCII 码值acos 反余弦 acosh 反双曲余弦 acot 反余切acoth 反双曲余切acsc 反余割 acsch 反双曲余割align 启动图形对象几何位置排列工具all 所有元素非零为真angle 相角ans 表达式计算结果的缺省变量名any 所有元素非全零为真area 面域图argnames 函数 M 文件宗量名asec 反正割asech 反双曲正割asin 反正弦 asinh 反双曲正弦assignin 向变量赋值atan 反正切 atan2 四象限反正切atanh 反双曲正切autumn 红黄调秋色图阵 axes 创建轴对象的低层指令axis 控制轴刻度和风格的高层指令B bbar 二维直方图bar3 三维直方图。

17、1程序代码：%普通的灰色预测GM1.mclearX=input(请输入原始数据:,s);%原始数据（可以多行，每一行为一类原始数据，即可多类原始数据）X=str2num(X);m1m2=size(X);%m1和m2分别表示X的行数和列数k0=input(请输入所要预测的阶数:);%GM(1,1)模型fori=1:m1n=i;x0=X(i,:);%将原始数据X中的第i行数据赋给x0，即取出一类原始数据disp(1.原始数据:);Y=;forz=1:m2Y=strcat(Y,(,num2str(x0(z),);enddisp(Y);%1.利用一次累加（1-AGO）生成新数列E=triu(ones(m2);%E表示元素为1右上三角阵x1=x0*E;%对原始数据进行一次累加（1-AGO）生成新数列x1disp(2.一。

18、哈哈哈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);边缘检测：globa。

19、MATLAB主要命令汇总 MATLAB函数参考 附录1 1 管理用命令 函数名 功能描述 函数名 功能描述 addpath 增加一条搜索路径 rmpath 删除一条搜索路径 demo 运行Matlab演示程序 type 列出 M文件 doc 装入超文本文档 version 显示Matlab的版本号 help 启动联机帮助 what 列出当前目录下的有关文件 lasterr 显示最后一条信息 wha。

20、8第2章 图像获取2.3.2 二维连续傅里叶变换例2.2 figure(1); %建立图形窗口 1u,v = meshgrid(-1:0.01:1); %生成二维频域网格F1 = abs(sinc(u.*pi); F2 = abs(sinc(v.*pi);F=F1.*F2; %计算幅度频谱 F=|F(u,v)|surf(u,v,F); %显示幅度频谱，如图 2.3(b)shading interp; %平滑三维曲面上的小格axis off; %关闭坐标系figure(2); %建立图形窗口 2F1=histeq(F); %扩展 F 的对比度以增强视觉效果imshow(F1); %用图像来显示幅度频谱，如图 2.3(c)第3章 图像变换3.4.4 二维 FFT的 MATLAB实现例3.2 简单图像及其傅里叶变换 MATLAB程序：%建立简单图像。