matlab拉普拉斯滤波器程序

function projectdemo(n,m,p)

if mod(n,floor(n))~=0 | n<=2 | mod(n,2)==0

error('n必须是大于等于3的奇数整数')

end

if mod(m,floor(m))~=0 | m<=2 | mod(m,2)==0

error('m必须是大于等于3的奇数整数')

end

if mod(p,floor(p))~=0 | p<=2 | mod(p,2)==0

error('p必须是大于等于3的奇数整数')

end

center1=-((n^2)-1)

lapmask1=ones(n)

lapmask1((n+1)/2,(n+1)/2)=center1

center2=-((m^2)-1)

lapmask2=ones(m)

lapmask2((m+1)/2,(m+1)/2)=center2

center3=-((p^2)-1)

lapmask3=ones(p)

lapmask3((p+1)/2,(p+1)/2)=center3

fin=im2double(imread('moon.tif'))

gout1=fin-imfilter(fin,lapmask1,'replicate')

gout2=fin-imfilter(fin,lapmask2,'replicate')

gout3=fin-imfilter(fin,lapmask3,'replicate')

S0=char('*')S1=num2str(n)S2=char('滤波图像')

S=strcat(S1,S0,S1,S2)

subplot(2,2,1)imshow(fin)axis on

title('原始图像','FontSize',14,'FontName','隶书','color','b')

subplot(2,2,2)imshow(gout1)axis on

title('3*3滤波','FontSize',14,'FontName','隶书','color','b')

subplot(2,2,3)imshow(gout2)axis on

title('5*5滤波','FontSize',14,'FontName','隶书','color','b')

subplot(2,2,4)imshow(gout3)axis on

title('7*7滤波','FontSize',14,'FontName','隶书','color','b')

t=(n-(n/1048576)+1048576)/1048576

令u1和u2＝0+（1-0）*t

若u1小等于0.5则输出 -（β+In （1.-u2））

若u1大于0则输出 （β*In （u2））

NT-D系统 即所谓的 Destroy Mode，当Destroy Mode开启时，机体会发生大幅度的变化，由先前的独角兽变身为毁灭者，机体内部的精神感应框架扩张，体型完全改变，装甲结合部露出散发着诡异光芒的精神力框体，全高由19.7M变为21.7M，头部的独角（即NT-D感应装置）一分为二成为“V”字形，封闭式面罩完全打开，转变为传统意义上的“GUNDAM”造型。虽然Destroy Mode带来巨大的性能提高，不过也使得驾驶员的身体与大脑的双重负担超乎常人想象。因此Destroy Mode的运作极限也只有5分钟左右。根据现有官方公开的资料，独角兽在毁灭模式下与独角模式的注重远程攻击相反，此模式极大的强化了近战攻击力和机动力，一般将弃用光束步q等传统远程武器

---