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

jiulquq

我下面的MATLAB程序，我想输出三个结果，为啥结果只会出现一个，就是b,不会出现其他的结果呢？不知道是什么原因，请大神帮帮我，谢谢了。
+ _% O' P* u2 c$ M% _, k) }function[b,e,h]=lvboshibie(a)
a=imread('lena.bmp');
3 B8 }" C, f: W% U0 B, s; f/ ?%figure;imshow(a);
( [1 l% r% T" L8 S" \%title('原图像');
9 U" A% w8 k$ c; fa=im2double(a);%将数变为double型数
b=imnoise(a,'gaussian',0.015);
  z5 k; s; _# h5 {' Z' ^* sC=im2double(b);
F=im2col(C,[3,3],'distinct');%将加躁图像分割成3*3的小窗口
[i,j]=size(F);
m=zeros(i,j);
n=size(i,j);
' V# A6 j3 y2 M0 @/ ~! Tfor i=1:9
" Q( P0 D* A2 G: z/ ~: Sfor j=1:5180
% M5 x* K/ j- P# lm(i,j)=F(i,j)-F(5,j);%用3*3的窗口中的中心值对窗口中的数据求出估计的受躁程度
m(5,j)=0.001;
# k# w5 N; p) Z" r/ [: zn(i,j)=m(i,j)./F(5,j);
end
end
%为用accumarray函数把数据变回原来的排列方式做准备
%制造c矩阵
! l5 M0 ~6 |4 h5 v3 s  v# V/ hK=[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];
) m- C. {7 ~% m( m( iA=K;
for ii=1:5179;
' C! x$ m# `3 NA=[A;K];
end
g=zeros(1,5180);
F=[1:5180];
- B$ f; A% {5 H' q/ ^! gF=[g,g,F];
F=[F;F;F;F;F;F;F;F;F];
' v9 P% l  `5 oF=reshape(F,[46620,3]);
c=A+F;
  g* e& S/ A' I/ d# P) n( aval=reshape(n,[46620,1]);
: S+ C- ]' G# X5 U: n: x: D* d, y0 dA=accumarray(c,val);%使用此函数把原来用im2col变成9*5180的矩阵变为一个个3*3的小矩阵
: d  c, {9 c1 x; WD=reshape(A,[3,15540]);
%将矩阵变为原来的222*210的形式
$ t2 E6 p( m3 LM = [];
B=[];
for i = 1:3:15538
8 H( n/ o) X" h5 D# x4 G. M2 Wtemp=D(:,i:i+2);

M = [M; temp]; % 使矩阵变为15540*3
end
* _  k8 \+ h3 j2 v! tfor j=1:222:15539

TEMP=M(j:j+221,;
+ `  W1 u- t8 H# UB = [B, TEMP]; % 使矩阵变为220*210
0 K1 ~7 N" ^) a0 f/ G) hend
%使矩阵为222*208
& Z) V: {& R8 ]# T/ K9 g- A& rB=B(1:222,1:208);
) z1 B! \. y; H9 W7 AE=B+100;
d=100.*a;%求出估计的灰度时用数据
# Q8 S$ f  h1 ^* Z5 m3 Qr=d./E;
e=im2uint8(r);
  O6 S) y' Q2 B8 d" B! ?[m,n]=size(e);
8 E) \5 o" d0 Dfor i=1:m
$ c, ]3 h4 d9 \8 l* r! O" qfor j=1:n
5 Y$ f0 e+ t) A. mif (e(i,j)>=0)&&(e(i,j)<=50)
u=23;
- J5 X' ?' s6 l4 v4 lq=23/3;
* I- {( i+ u( k5 H! a+ q4 _v(i,j)=e(i,j)-u;
; v$ A* S/ i' Y0 I0 yV(1)=2*q^2;
elseif (e(i,j)>=51)&&(e(i,j)<=100)
u=83;
q=47/3;
: e) f3 H3 }6 w- |v(i,j)=e(i,j)-u;
V(2)=2*q^2;
else
4 r( W  f) e/ n; C* [5 z* mu=117;
3 Y% H$ ]( h& Uq=138/3;
v(i,j)=e(i,j)-u;
V(3)=2*q^2;
end
. e6 Y: s) _5 K- G$ N+ |# u9 jend
end
v=im2double(v);
Q=v.^2;
3 k0 ^/ p- q+ w) Vfor i=1:m
7 d/ O3 C$ p1 T+ p) |0 v, afor j=1:n
if (e(i,j)>=0)&&(e(i,j)<=50)
H(i,j)=im2double(Q(i,j)/V(1));
N(i,j)=exp(-H(i,j));
* v. r1 K2 u# K, ]# w9 Yelseif (e(i,j)>=51)&&(e(i,j)<=100)
H(i,j)=im2double(Q(i,j)/V(2));
N(i,j)=exp(-H(i,j));
else
H(i,j)=im2double(Q(i,j)/V(3));
N(i,j)=exp(-H(i,j));
5 o2 ^- \1 l# Fend
end
. m$ c* a5 I$ h; Q6 q' hend
s=ones(224,210);
% X; n" o! u8 j% X  @8 \7 al=s*26;
" z( }+ Z5 ~: W* ml(2:223,2:209)=e;
4 I1 o1 A0 [* j2 l7 E1 lo=zeros(224,210);
o(2:223,2:209)=N;
  Z3 ~- A- B# R& O8 _. I6 N( M[m,n]=find(l<=25);
q=size(m);
for i=1:q
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);
y(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));
$ Q3 j5 Z+ s+ J& e9 \end
h=l(2:223,2:209);
h=round(h);
0 x9 F( f! `. @: _% S! e%figure,imshow(h,[]);
8 O0 m' ]2 U2 [5 v: y+ e$ c3 @# J%title('第二次去噪效果');
) Q% v) u' A8 b) U* n9 ^) Bb
e
h
/ p$ i, O7 r9 _2 `! G5 tend

# B1 s7 L1 G4 S( s1 d" K5 e6 v; @

sharkN

转转