用MATLAB的function文件想出三个结果，但是为啥只有一个？

jiulquq

我下面的MATLAB程序，我想输出三个结果，为啥结果只会出现一个，就是b,不会出现其他的结果呢？不知道是什么原因，请大神帮帮我，谢谢了。
, D/ z) Z$ b0 m# u) nfunction[b,e,h]=lvboshibie(a)
7 L" G) q" t3 a4 M7 y0 s* i+ `a=imread('lena.bmp');
%figure;imshow(a);
6 I( `9 X: V& Z' _" j3 g2 W%title('原图像');
a=im2double(a);%将数变为double型数
7 b7 T' O" |" r3 B5 Y: _b=imnoise(a,'gaussian',0.015);
C=im2double(b);
F=im2col(C,[3,3],'distinct');%将加躁图像分割成3*3的小窗口
[i,j]=size(F);
m=zeros(i,j);
- s" |; P2 o1 Qn=size(i,j);
4 d, ~( K7 M' d6 Jfor i=1:9
for j=1:5180
) J* z; S4 e" w0 ]m(i,j)=F(i,j)-F(5,j);%用3*3的窗口中的中心值对窗口中的数据求出估计的受躁程度
m(5,j)=0.001;
- D3 w  X3 w  B- Bn(i,j)=m(i,j)./F(5,j);
end
: Z1 M" t/ n& ~  t  f: Gend
%为用accumarray函数把数据变回原来的排列方式做准备
%制造c矩阵
K=[1 1 0;1 2 0;1 3 0;2 1 0;2 2 0;2 3 0;3 1 0;3 2 0;3 3 0];
A=K;
for ii=1:5179;
A=[A;K];
" q. @0 s( D9 y, G% c- C4 Tend
5 E$ z" j$ p6 ^" t! h# g% W& Gg=zeros(1,5180);
F=[1:5180];
" w0 ?( H5 G# _1 p, `F=[g,g,F];
$ P) b$ p' c+ L  y5 j( S1 SF=[F;F;F;F;F;F;F;F;F];
F=reshape(F,[46620,3]);
4 |( F* }$ @. o5 l0 L: ]/ Dc=A+F;
val=reshape(n,[46620,1]);
$ `( X, v# ?' aA=accumarray(c,val);%使用此函数把原来用im2col变成9*5180的矩阵变为一个个3*3的小矩阵
D=reshape(A,[3,15540]);
# R/ l9 @2 v8 C* W8 w4 I2 H  u%将矩阵变为原来的222*210的形式
* M; J$ a+ }- f3 S: }  ]M = [];
! Z$ G+ O6 b, T. T1 }; K. E3 ~B=[];
2 B1 G3 L4 c; U$ Pfor i = 1:3:15538
" o% ^6 q+ \2 E+ Ftemp=D(:,i:i+2);

M = [M; temp]; % 使矩阵变为15540*3
1 ]# b( l# _) Mend
! N0 n5 y6 T  U  y( g0 e' Yfor j=1:222:15539
" t- n$ q: Y+ s
  f* F8 s+ e3 b8 k9 v# XTEMP=M(j:j+221,;
4 c1 N3 i( s* Z; T: WB = [B, TEMP]; % 使矩阵变为220*210
1 x1 M" k' v  v, p# Q$ mend
%使矩阵为222*208
5 y4 w$ ^5 a& p, {& T6 vB=B(1:222,1:208);
: B' L+ d+ K2 I8 x# g. Z& t% g5 GE=B+100;
5 E2 [1 Z, }8 U) m8 v6 Gd=100.*a;%求出估计的灰度时用数据
r=d./E;
e=im2uint8(r);
[m,n]=size(e);
for i=1:m
for j=1:n
if (e(i,j)>=0)&&(e(i,j)<=50)
u=23;
  }9 Y$ A' H" g/ _' o4 fq=23/3;
v(i,j)=e(i,j)-u;
V(1)=2*q^2;
" W" \6 y9 n' q3 B" D6 \* p/ Kelseif (e(i,j)>=51)&&(e(i,j)<=100)
8 O8 y. A8 p! ~* y; D9 Eu=83;
( w+ p/ x8 s3 _5 tq=47/3;
) v6 r+ D: H' o, ?v(i,j)=e(i,j)-u;
9 A& D! G( v( }/ Z! KV(2)=2*q^2;
else
1 e7 ^% D& |: L: ~: I+ ?& _u=117;
q=138/3;
v(i,j)=e(i,j)-u;
1 b# R' ?' J9 rV(3)=2*q^2;
. ]! n9 |6 g- o; Z" p# eend
4 @& o  N, y7 v/ r9 v& q, Z( b1 yend
end
9 ]: y- Y6 F) F$ S4 Ev=im2double(v);
6 X2 s7 H. |+ e2 y" W5 `Q=v.^2;
for i=1:m
for j=1:n
if (e(i,j)>=0)&&(e(i,j)<=50)
8 ~! B# \* E$ d/ @H(i,j)=im2double(Q(i,j)/V(1));
5 ?" J8 E4 ^4 S, uN(i,j)=exp(-H(i,j));
elseif (e(i,j)>=51)&&(e(i,j)<=100)
H(i,j)=im2double(Q(i,j)/V(2));
! K* K7 w  p7 m6 fN(i,j)=exp(-H(i,j));
else
4 p! V4 E0 H: c% BH(i,j)=im2double(Q(i,j)/V(3));
N(i,j)=exp(-H(i,j));
6 ]0 f7 q6 A! @/ Yend
end
end
s=ones(224,210);
l=s*26;
l(2:223,2:209)=e;
( B. d8 e3 N& P# [o=zeros(224,210);
o(2:223,2:209)=N;
[m,n]=find(l<=25);
q=size(m);
* L& Z* ~" Z( a& Wfor i=1:q
" \& x6 a1 H" B  y1 p, r. tx(m(i),n(i))=l(m(i)-1,n(i)-1)*o(m(i)-1,n(i)-1)+l(m(i)-1,n(i))*o(m(i)-1,n(i))+l(m(i)-1,n(i)+1)*o(m(i)-1,n(i)+1)+l(m(i),n(i)-1)*o(m(i),n(i)-1)+l(m(i),n(i))*o(m(i),n(i))+l(m(i),n(i)+1)*o(m(i),n(i)+1)+l(m(i)+1,n(i)-1)*o(m(i),n(i)-1)+l(m(i)+1,n(i))*o(m(i)+1,n(i))+l(m(i)+1,n(i)+1)*o(m(i)+1,n(i)+1);
8 _- y. n2 B" y+ _1 x, cy(m(i),n(i))=o(m(i)+1,n(i))+o(m(i),n(i))+o(m(i)-1,n(i))+o(m(i)+1,n(i)-1)+o(m(i)-1,n(i)-1)+o(m(i),n(i)-1)+o(m(i)-1,n(i)+1)+o(m(i),n(i)+1)+o(m(i)+1,n(i)+1);
l(m(i),n(i))=x(m(i),n(i))/y(m(i),n(i));
end
/ U; L3 f) P; f1 Bh=l(2:223,2:209);
& ?9 t; l( f1 b: c6 [' Ch=round(h);
%figure,imshow(h,[]);
* F. m# J& [. A. V  f! n, h  h- V%title('第二次去噪效果');
' R' t  K+ y9 x( A, T' ^% l! ]b
e
# N7 B! s+ U( u+ j( d# Fh
end

. P; i: |# u$ T8 ?) @$ }6 v" v

sharkN

转转