随笔分类 - Matlab练习程序
摘要：clear all; close all; clc; n=40; point=load('point.mat'); current_point=point.point; plot(current_point(:,1),current_point(:,2),'r'); %系统方程：x(k+1)=fi*
阅读全文
摘要：zseven.mfunction out=zseven(nbhd) s=sum(nbhd(:))-nbhd(5); temp1=(2<=s)&(s<=6); p=[nbhd(1) nbhd(4) nbhd(7) nbhd(8) nbhd(9) nbhd(6) nbhd(3) nbhd(2)]; pp=[p(2:8) p(1)]; xp=sum((1-p).*pp); temp2=(xp==1); prod1=nbhd(4)*nbhd(8)*nbhd(2); prod2=nbhd(4)*nbhd(6)*nbhd(2); temp3=(prod...
阅读全文
摘要：clear all;close all;clc;randn('seed',0);%%一维高斯函数mu=0;sigma=1;x=-6:0.1:6;y=normpdf(x,mu,sigma);plot(x,y);figure;%%二维或多维高斯函数mu=[0 0];sigma=[0.3 0;0 0.35];[x y]=meshgrid(linspace(-8,8,80)',linspace(-8,8,80)');X=[x(:) y(:)];z=mvnpdf(X,mu,sigma);surf(x,y,reshape(z,80,80));hold on;%再生成一个mu
阅读全文
摘要：clear all;close all;clc;randn('seed',0);%{一维情况mu=0;N=100000;S=5;data=mvnrnd(mu,S,N);me=mean(data);S2=1/N*sum((data-me).^2);%}%%二维或多维情况mu=[0 0];S=[0.3 0;0 0.35];N=10000;data=mvnrnd(mu,S,N);me=mean(data)S1=1/N*sum((data(:,1)-me(1)).^2)S2=1/N*sum((data(:,2)-me(2)).^2)
阅读全文
摘要：clear all; close all; clc; randn('seed',0); mu1=[0 0]; S1=[0.3 0;0 0.35]; cls1_data=mvnrnd(mu1,S1,1000); plot(cls1_data(:,1),cls1_data(:,2),'+'); hold
阅读全文
摘要：clear all;close all;clc;num=4; %元素数量k=180; %迭代次数step=0.1; %迭代步长w(1,:)=[-0.5 1 1]; %权值x=[1 0 0; %输入的值，每行为一组 1 1 0; 1 0 1; 1 1 1];d=[1 0 1 0]; %实际输出的值y=zeros(1,num); %迭代输出的值 y=w0+w1x1+w2x2e=zeros(1,k); %均方误差for j=1:k tmp_2=0; for i...
阅读全文
摘要：clear all;close all;clc;num=7;X=[ 1 0.4 0.05; 1 -0.2 0.75; 1 0.5 0.1; 1 -0.1 0.7; 1 0.3 0.2; 1 -0.15 0.9; 1 0.5 0.2];d=[1 -1 1 -1 1 -1 1]; %同一类的的数据正负号相同for i=1:num plot(X(i,2),X(i,3),'ro'); hold on;endaxis([-0.5 0.5 0 1]);Y=X'*X;inv_Y=inv(Y);w=inv_Y*X'*...
阅读全文
摘要：clear all;close all;clc;step=0.1; %迭代步长num=15; %数据个数k=100; %迭代次数set=[ 0.4 0.05 1; - 0.2 0.75 1; 0.5 0.1 1; - 0.1 0.7 1; 0.3 0.2 1; - 0.15 0.9 1; 0 0.8 1; 0.4 0.3 1; -0.2 0.87 1; 0.4 ...
阅读全文
摘要：%LMS算法演示(matlab)%设置参数，N为采样个数，u为步长clear,clc;N=16;u=0.1;%设置迭代次数kk=1000;%pha为随机噪声的平均功率rk=randn(1,k)/2;%%正态分布的随机矩阵pha=mean(rk);%%求元素平均值%设置起始权值wk(1,:)=[0 0];%用LMS算法迭代求最佳权值for i=1:k xk(i,:)=[sin(2*pi*i/N) sin(2*pi*(i-1)/N)]+rk(i);%输入信号 yk(i)=xk(i,:)*wk(i,:)';%输出信号 dk(i)=2*cos(2*pi*i/N);%期望信号 e...
阅读全文
摘要：clear all;close all;clc;img1=imread('Corner.png');img2=imread('Corner1.png');img3=imread('Corner2.png');img1=imresize(img1,[35 90]); %Matlab的svd不支持太大的数据，我把数据变小了.img2=imresize(img2,[35 90]);img3=imresize(img3,[35 90]);[height width]=size(img2);a=reshape(double(img1),1,35*90);b
阅读全文
摘要：clear all;close all;clc;im=imread('2.png');im2=im;im1=rgb2ycbcr(im);%将图片的RGB值转换成YCbCr值%Lu=im1(:,:,1);Cb=im1(:,:,2);Cr=im1(:,:,3);[x y z]=size(im);tst=zeros(x,y);%计算Cb、Cr的均值Mb、Mr%Mb=mean(mean(Cb));Mr=mean(mean(Cr));%计算Cb、Cr的均方差%Db=sum(sum(Cb-Mb))/(x*y);Dr=sum(sum(Cr-Mr))/(x*y);%根据阀值的要求提取出near
阅读全文
摘要：clear all;close all;clc;img=imread('Lena_cai.jpg');imshow(img);mean1=mean(mean(img(:,:,1)));mean2=mean(mean(img(:,:,2)));mean3=mean(mean(img(:,:,3)));me=mean(mean(mean(img)));img(:,:,1)=me/mean1*img(:,:,1);img(:,:,2)=me/mean2*img(:,:,2);img(:,:,3)=me/mean3*img(:,:,3);figure,imshow(img);参考：ht
阅读全文
摘要：clear all;close all;clc;img=imread('lena.jpg');imshow(img);phi=pi/2;theta=2;sigma=0.65*theta;filterSize=6;G=zeros(filterSize);for i=(0:filterSize-1)/filterSize for j=(0:filterSize-1)/filterSize xprime=j*cos(phi); yprime=i*sin(phi); K=exp(2*pi*theta*sqrt(-1)*(xprime+yprime)); ...
阅读全文
摘要：clear all;close all;clc;img=imread('saber9.jpg');img2=imread('saber2.jpg');imshow(img)figure,imshow(img2);tmp=rgb2gray(img);tmp2=rgb2gray(img2);img_re=imresize(tmp,[8 8]);img_re2=imresize(tmp2,[8 8]);img_re=uint8(double(img_re)/4);img_re2=uint8(double(img_re2)/4);me=mean(mean(img_re)
阅读全文
摘要：clear all; close all; clc; img1=imread('Corner.png'); img2=imread('Corner1.png'); img3=imread('Corner2.png'); [height width]=size(img2); a=reshape(img
阅读全文
摘要：close all; clc; H=1; %索引pix中第一个元素，即高度 W=2; %索引pix中第二个元素，即宽度 left_right=0.3; %抬起左边或右边时值为0-1之间，不抬起时为0 up_down=0; %抬起上边或下边时值为0-1之间，不抬起时为0 img=imread('len
阅读全文
摘要：close all; clear all; clc; img=imread('rice.png'); imshow(img); [m n]=size(img); tmp=zeros(m+2,n+2); tmp(2:m+1,2:n+1)=img; Ix=zeros(m+2,n+2); Iy=zeros
阅读全文
摘要：close all;clear all;clc;img=imread('pic.bmp');imshow(img);[m n]=size(img);re=zeros(m,n);tmp=zeros(m+2,n+2);tmp(2:m+1,2:n+1)=img;for i=2:m+1 for j=2:n+
阅读全文
摘要：close all;clear all;clc;img=imread('1.jpg');imshow(img);[m n]=size(img);img_re=zeros(m,n);tmp=zeros(8,8);for i=1:8:m for j=1:8:n tmp=mc(img(i:i+7,j:j+7)); img_re(i:i+7,j:j+7)=tmp; endendfigure,imshow(uint8(img_re));figure,imshow(histeq(img));%{a=[238 238 236 237 236 228 227 234;...
阅读全文
摘要：cl;H=1; %索引pix中第一个元素，即高度W=2; %索引pix中第二个元素，即宽度jiaodu=0; %要旋转的角度，旋转方向为顺时针img=imread('Corner.png'); %这里v为原图像的高度，u为原图像的宽度imshow(img); %这里y为变换后图像的高度，x为变换后图像的宽度[v u]=size(img);h=v; ...
阅读全文