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

我下面的MATLAB程序，我想输出三个结果，为啥结果只会出现一个，就是b,不会出现其他的结果呢？不知道是什么原因，请大神帮帮我，谢谢了。
5 J  n. Q0 {+ a7 ffunction[b,e,h]=lvboshibie(a)6 Y, c; M/ k, v* Y* I4 g9 R# W
a=imread('lena.bmp');
! V4 e' x/ `. ^2 m%figure;imshow(a);
5 m8 K# A5 |3 Z4 I%title('原图像');
3 k9 y* ~8 h3 O3 Z. P8 Ta=im2double(a);%将数变为double型数# g1 X( y2 r" g$ U4 f
b=imnoise(a,'gaussian',0.015);
6 ?2 I2 ^+ I! r( [C=im2double(b);$ {1 b$ e% N$ ]1 a
F=im2col(C,[3,3],'distinct');%将加躁图像分割成3*3的小窗口
' ]: w" n8 d0 [; t( q[i,j]=size(F);
; |4 b% r, D; U  nm=zeros(i,j);
* ?, z2 L& x. r4 e/ x8 {$ q' yn=size(i,j);; W2 m) h% `, m3 T$ E
for i=1:9 2 n+ Z+ V7 p5 G0 w  H- q" e; Q6 I/ {
for j=1:5180
  W+ e( x+ \8 i( Sm(i,j)=F(i,j)-F(5,j);%用3*3的窗口中的中心值对窗口中的数据求出估计的受躁程度9 k+ ^2 D! P/ a! `5 h
m(5,j)=0.001;
1 ^  x$ p) L) X$ t! K. H2 `  rn(i,j)=m(i,j)./F(5,j);
' u4 a/ l9 G# |; X+ N4 M/ qend) F# l% ^  F5 j8 B0 \6 n  ]
end, G" b& f6 V9 D- x: U& m
%为用accumarray函数把数据变回原来的排列方式做准备
5 R6 i) ~& i7 `: U" H# Q%制造c矩阵
' o1 d) B: b- S/ BK=[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];
- c$ I7 h+ H! U4 c- s4 L/ i7 W9 yA=K;
3 I; E( J1 v4 n- S3 s4 I- qfor ii=1:5179;7 A* I; z; C& n" K' D2 a- w
A=[A;K];5 d) l# I2 b$ a2 x% H% p
end
4 r5 I7 W* a( `0 Zg=zeros(1,5180);. x- h. `" n0 D% x0 P, ?; w
F=[1:5180];
$ F1 @& a' X( y9 t2 ?F=[g,g,F];  F* N7 X6 |& X) U* b, \/ F
F=[F;F;F;F;F;F;F;F;F];- g! ^4 Z" U9 A0 b
F=reshape(F,[46620,3]);+ i2 w( U; T3 C7 b( E" R
c=A+F;
, I4 y8 a3 _* a  d$ o7 x: ?9 Y# |val=reshape(n,[46620,1]);
( }7 E/ C6 P) |1 K) x* ?A=accumarray(c,val);%使用此函数把原来用im2col变成9*5180的矩阵变为一个个3*3的小矩阵4 e' {2 j. q2 N/ h& M6 v2 [5 K
D=reshape(A,[3,15540]);
8 \/ D5 l* @0 i" S: u%将矩阵变为原来的222*210的形式4 f! \$ Q8 a( f! \# A& M9 P: O8 [
M = [];
" }/ B* U9 M% P, d2 M$ @B=[];4 A. t9 M5 Q6 k. i. H2 [
for i = 1:3:15538
; ]* X1 j3 ?4 V4 w7 @  p# ttemp=D(:,i:i+2);

M = [M; temp]; % 使矩阵变为15540*3
; T# _+ \$ Y" |1 {2 `; pend
8 n+ y# R* u; N9 xfor j=1:222:15539* ?- ~# J* I, y! @" @  _
  K7 z" ~* O" ~! |
TEMP=M(j:j+221,;
8 Z8 a& k- s/ i; B# W% Y) IB = [B, TEMP]; % 使矩阵变为220*210
! X4 f+ @: x2 V0 m8 oend
) R4 u/ r2 \8 l6 H( p- L%使矩阵为222*208% l0 n. I$ m1 c9 f3 F$ d
B=B(1:222,1:208);# _( k1 h9 C# d
E=B+100;+ }8 f- R- ^  G/ X, U- ]
d=100.*a;%求出估计的灰度时用数据
7 Z4 Q8 I/ Y& r( `6 V9 H! z& x$ Sr=d./E;! ~& z0 A: s1 u6 [, l6 E; g
e=im2uint8(r);
. X+ ~# E( W# A+ {[m,n]=size(e);$ J9 r/ N  c! _( K/ Z
for i=1:m/ q& L# V! x& _2 \
for j=1:n
9 K; ~& e3 h9 c2 D- W, Jif (e(i,j)>=0)&&(e(i,j)<=50)
3 R6 e2 {. q% }u=23;/ b' B, g1 s# N) u
q=23/3;9 ~+ d( d. t2 N1 g* B
v(i,j)=e(i,j)-u;4 }$ s1 ^) Z# Q) n
V(1)=2*q^2;
! k4 q$ ~- {3 J" D7 [& Gelseif (e(i,j)>=51)&&(e(i,j)<=100)
& l# k; F, _3 K4 z9 Q7 {: {u=83;; @/ q3 j' L' J' Q; h. t
q=47/3;
' M# ^& D# V/ T) W1 T2 ^v(i,j)=e(i,j)-u;6 _/ |# h( P8 c
V(2)=2*q^2;7 Y) C% F! S+ Q9 F' R# {* f( }
else # ?! K2 ]( M6 `, c1 z% v7 o& P
u=117;
. O, n8 C+ D  q4 x, H; {9 m& nq=138/3;
/ V5 I. U' _8 N, r. V3 ?v(i,j)=e(i,j)-u;
! `+ R" y) s' ]- g% R! ~7 S, `V(3)=2*q^2;! F, L- L( k) y* J
end
2 [9 C, `% N5 cend4 t2 v; ~7 W: X) M" C2 C! ~3 s
end$ v1 u& s' |. x. M( L" E' Z
v=im2double(v);
2 a2 N; J  r9 r/ tQ=v.^2;
! J) z6 l/ e# r8 m# ~7 n9 r; Jfor i=1:m
2 a9 t/ z& o4 {! A5 |/ afor j=1:n- i0 G. ]6 u/ }
if (e(i,j)>=0)&&(e(i,j)<=50)
8 L! c6 a6 [' _: C- l! B0 cH(i,j)=im2double(Q(i,j)/V(1));3 H4 t% Q! Z3 u
N(i,j)=exp(-H(i,j));
1 V, G" W: ^( yelseif (e(i,j)>=51)&&(e(i,j)<=100)
& l9 J* M# v, Q- \H(i,j)=im2double(Q(i,j)/V(2));
" ?3 [- Q% H( }) z7 {, |N(i,j)=exp(-H(i,j));
& o* N: n* u# c. r- f  ?/ Z9 celse
, @( ]7 F0 C4 n  l0 m* PH(i,j)=im2double(Q(i,j)/V(3));
) z4 q! t- F; d( v: Q, K0 ]N(i,j)=exp(-H(i,j));: Z+ L2 U) @9 k# o7 q- F1 Y7 V+ j' B
end" _) M  \4 Q) d
end
2 m* f- Q/ K9 ]: Rend5 t+ s6 j# @0 T0 C4 j
s=ones(224,210);+ V. V  z2 k# L
l=s*26;+ K8 C) Y. a- ]# ?4 q2 C3 E
l(2:223,2:209)=e;6 [3 O7 F4 l0 h( h; u
o=zeros(224,210);- |7 F8 A% J& q; o' }
o(2:223,2:209)=N;' g$ z; U5 T( A* N3 |' X
[m,n]=find(l<=25);- Q2 F% A, a$ C
q=size(m);9 `7 v" j$ {, W6 ?3 I. ]5 S! N
for i=1:q7 |: `' c& {; ~- L
x(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);
; U' b* M5 r9 F0 L2 qy(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);
; I6 I) G3 C6 f: _% V8 `9 Ml(m(i),n(i))=x(m(i),n(i))/y(m(i),n(i));
& R- Y" \' T" i& O6 ^end
* R: [2 k8 q- j& w$ k& |- z) r; e$ _h=l(2:223,2:209);
, @. W9 K) k6 R4 Th=round(h);
) r5 F3 q' F9 h3 [9 j3 g, |%figure,imshow(h,[]);
/ |: ?8 r* @! _' o%title('第二次去噪效果');1 x$ ^$ j$ g+ W% G: C8 w+ X" N
b
* c) K1 v( j8 A# M6 y# `8 U: i3 u# we# e3 x5 Q  D5 t4 |* K! a. t+ _
h6 H5 d) g7 h# v% x6 l
end