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

jiulquq

我下面的MATLAB程序，我想输出三个结果，为啥结果只会出现一个，就是b,不会出现其他的结果呢？不知道是什么原因，请大神帮帮我，谢谢了。
function[b,e,h]=lvboshibie(a)
a=imread('lena.bmp');
%figure;imshow(a);
) T+ w( w7 j: B# J%title('原图像');
. r- V, G  f' na=im2double(a);%将数变为double型数
b=imnoise(a,'gaussian',0.015);
, K5 d2 Q. g( l$ aC=im2double(b);
, u2 c) H, z6 G5 ?' \9 H6 EF=im2col(C,[3,3],'distinct');%将加躁图像分割成3*3的小窗口
9 L; Q/ N/ i5 L0 X6 o. n9 {. y- O[i,j]=size(F);
m=zeros(i,j);
n=size(i,j);
# l- z3 N- ]- S! D1 ]! Xfor i=1:9
for j=1:5180
; _' r, ?# d# a% vm(i,j)=F(i,j)-F(5,j);%用3*3的窗口中的中心值对窗口中的数据求出估计的受躁程度
m(5,j)=0.001;
2 F) A- y; u! N* C* R& |4 cn(i,j)=m(i,j)./F(5,j);
end
& H# y0 Q* a* W5 |; Q* e1 Tend
0 T8 A, o2 T) v/ x6 Q. ]% g%为用accumarray函数把数据变回原来的排列方式做准备
% L9 H7 u) T9 y$ |. \%制造c矩阵
8 Z: [8 x0 L) K8 D9 @* i8 X2 XK=[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];
% X$ }5 n$ P3 D, `" N# c8 }A=K;
for ii=1:5179;
* R" _: _5 K* L" |% R: t4 mA=[A;K];
end
$ I) L+ [( t' wg=zeros(1,5180);
F=[1:5180];
' I3 i' d1 k  t; TF=[g,g,F];
F=[F;F;F;F;F;F;F;F;F];
F=reshape(F,[46620,3]);
3 A: I+ H3 v7 h& F$ F3 {6 uc=A+F;
! f' L+ c- v- x9 ~! m) T' y6 Aval=reshape(n,[46620,1]);
A=accumarray(c,val);%使用此函数把原来用im2col变成9*5180的矩阵变为一个个3*3的小矩阵
* K& w0 y$ `3 E& v! S' hD=reshape(A,[3,15540]);
%将矩阵变为原来的222*210的形式
' L- _! R) ]. r0 J6 oM = [];
, \/ i3 \: V+ d& f! }B=[];
for i = 1:3:15538
temp=D(:,i:i+2);

M = [M; temp]; % 使矩阵变为15540*3
end
for j=1:222:15539

TEMP=M(j:j+221,;
' p; ^0 s* D3 x& F) OB = [B, TEMP]; % 使矩阵变为220*210
end
%使矩阵为222*208
B=B(1:222,1:208);
3 M6 x  a8 o2 _, CE=B+100;
/ ~& D  }; i# e9 C* e; Z  Ad=100.*a;%求出估计的灰度时用数据
  F8 @+ ]9 P/ F, T; D  C  ]0 Hr=d./E;
e=im2uint8(r);
[m,n]=size(e);
for i=1:m
2 n0 r' O" Z9 l# N( E3 N+ Dfor j=1:n
% m! E8 \2 U4 i2 Zif (e(i,j)>=0)&&(e(i,j)<=50)
u=23;
' h9 A' \: D0 l7 M/ Fq=23/3;
& k* i2 i5 [% d7 T3 R, r* Q, B# I& l- gv(i,j)=e(i,j)-u;
V(1)=2*q^2;
elseif (e(i,j)>=51)&&(e(i,j)<=100)
u=83;
q=47/3;
v(i,j)=e(i,j)-u;
0 [: G$ N* C0 Z% |; U% @$ E7 ~V(2)=2*q^2;
else
# `  u$ z( y& ]u=117;
q=138/3;
- g3 g9 [- _5 q/ v5 `v(i,j)=e(i,j)-u;
' g' A: n6 _7 Q* QV(3)=2*q^2;
3 |% t2 s) z! X/ H6 I  Wend
$ I: i# q" m) i$ A( _0 ^end
: r/ u' y1 l7 {9 X0 dend
( t% @3 i) r' zv=im2double(v);
+ M: E$ j: F5 RQ=v.^2;
! |' @$ p$ f( m8 H. wfor i=1:m
for j=1:n
/ M& u/ _; z) L" _if (e(i,j)>=0)&&(e(i,j)<=50)
' R. V5 Z# X1 {3 ~1 sH(i,j)=im2double(Q(i,j)/V(1));
N(i,j)=exp(-H(i,j));
" |  A$ b0 d" ~+ V% u2 L& selseif (e(i,j)>=51)&&(e(i,j)<=100)
9 f# g7 B. U- d( _H(i,j)=im2double(Q(i,j)/V(2));
N(i,j)=exp(-H(i,j));
/ h: r  q8 Z, Zelse
H(i,j)=im2double(Q(i,j)/V(3));
N(i,j)=exp(-H(i,j));
end
  z3 M, _1 I* eend
end
) w" u. `+ j. D! U* H# h8 ms=ones(224,210);
2 G0 p% r& z3 i' F* }l=s*26;
l(2:223,2:209)=e;
8 P8 W! W5 ^( x; j8 w: @! a6 yo=zeros(224,210);
/ a4 l8 t$ S1 `( i4 Jo(2:223,2:209)=N;
3 S- f' S9 g+ m3 g* M[m,n]=find(l<=25);
6 E% h* G2 t* d. G& oq=size(m);
; ~0 O) x. S9 Dfor i=1:q
  B1 @" I, |, A& ]! l4 |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);
! ~0 a% X, U7 p, Y6 e8 v& Gy(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));
4 b# [+ ~7 _) J+ |end
h=l(2:223,2:209);
h=round(h);
4 o2 M, V0 j/ W$ j%figure,imshow(h,[]);
8 b* R8 ^% \0 q" B%title('第二次去噪效果');
  x* \; H0 e5 lb
/ L5 k. n7 z/ w/ T# i6 |e
; l1 S5 I' J' X0 Kh
* V, N- E0 U* I% o; hend

sharkN

转转