(完整word版)matlab经典代码大全.doc

哈哈哈
MATLAB

显示正炫余炫图：plot(x,y1,'*r',x,y2,'ob')
定义【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);
边缘检测：
globalim

str=get(hObject,'string');

axes(handles.axes1);

switchstr
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];%StructuringElement
I=rgb2gray(im);
imshow(I,[]);title('OriginalImage');
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');%Imagetobeusedfordilation
image_erosion=padarray(I,se_origin,256,'both');%Imagetobeusedforerosion
%%%%%%%%%%%%%%%%%%
%%%Dilation%%%
%%%%%%%%%%%%%%%%%%
fork=se_origin(1)+1:im_height+se_origin(1)
forkk=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+halfheight-1,kk-se_origin(2):kk+halfwidth-1)));
end
end
figure;imshow(dilated_image,[]);title('ImageafterDilation');
%%%%%%%%%%%%%%%%%
%%%Erosion%%%
%%%%%%%%%%%%%%%%%
se=se';
fork=se_origin(2)+1:im_height+se_origin(2)
forkk=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('ImageafterErosion');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%Opening(Erosionfirst,thenDilation)%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
se=se';
image_dilation2=eroded_image;%Imagetobeu