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

jiulquq

我下面的MATLAB程序，我想输出三个结果，为啥结果只会出现一个，就是b,不会出现其他的结果呢？不知道是什么原因，请大神帮帮我，谢谢了。
9 n/ {4 J0 ^+ D8 X1 o2 Jfunction[b,e,h]=lvboshibie(a)
a=imread('lena.bmp');
& T! D% C2 @, m- I2 p2 f1 O8 h%figure;imshow(a);
3 U+ C; G0 p( J8 G. a& V%title('原图像');
( m0 B5 A9 @/ v& p" B9 v! R" Na=im2double(a);%将数变为double型数
b=imnoise(a,'gaussian',0.015);
$ O$ f* r* X3 o* k5 Y7 kC=im2double(b);
( d5 [+ i: y5 R6 P$ JF=im2col(C,[3,3],'distinct');%将加躁图像分割成3*3的小窗口
0 C3 {; S8 S$ J$ G[i,j]=size(F);
" W# v% n1 L" X8 l" N$ V5 @m=zeros(i,j);
/ J" I) O* L8 n. P. a% qn=size(i,j);
$ D. J4 U3 }; m% `$ ifor i=1:9
- @: K0 f* T; @8 b" Z; m1 Kfor j=1:5180
9 A: \/ A3 a9 A4 A6 Qm(i,j)=F(i,j)-F(5,j);%用3*3的窗口中的中心值对窗口中的数据求出估计的受躁程度
m(5,j)=0.001;
n(i,j)=m(i,j)./F(5,j);
2 e5 d) f6 f, c' g1 a6 A) \/ wend
end
%为用accumarray函数把数据变回原来的排列方式做准备
%制造c矩阵
! H. }$ U. A% z; G' qK=[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];
9 t: P$ F$ {9 ~. f8 vA=K;
for ii=1:5179;
) l) v: v2 b- ]* L4 e9 A+ @. EA=[A;K];
) H$ }! u4 ]+ U4 F# ]$ w) [+ ?5 D9 Y, y2 |end
g=zeros(1,5180);
6 A+ f& B6 i; k2 k! c' E3 kF=[1:5180];
# L5 U0 L. W! G# n) AF=[g,g,F];
3 `2 O: w' ?- n4 j- x4 k" ~F=[F;F;F;F;F;F;F;F;F];
. B0 M  l+ R, h* WF=reshape(F,[46620,3]);
; P5 z( O# V/ l) U8 l( Hc=A+F;
3 \. F! i. {, [# m* o, r/ Qval=reshape(n,[46620,1]);
A=accumarray(c,val);%使用此函数把原来用im2col变成9*5180的矩阵变为一个个3*3的小矩阵
D=reshape(A,[3,15540]);
%将矩阵变为原来的222*210的形式
  W1 C* o+ @9 M/ l- fM = [];
6 ^9 u& ~$ \" }9 HB=[];
1 y0 `% v- F: s  `for i = 1:3:15538
temp=D(:,i:i+2);

M = [M; temp]; % 使矩阵变为15540*3
+ ^# d4 \9 ~8 g$ Uend
! F2 u% K: U9 K0 _. lfor j=1:222:15539
7 G7 x; l8 I) A$ r( _- u
TEMP=M(j:j+221,;
B = [B, TEMP]; % 使矩阵变为220*210
end
0 }5 |& I8 M- J. A% _3 n%使矩阵为222*208
" T' E- M8 s! L& Y$ z% D! G  [/ KB=B(1:222,1:208);
: K5 t( ~1 E* z% \. `% d% G+ iE=B+100;
d=100.*a;%求出估计的灰度时用数据
r=d./E;
e=im2uint8(r);
[m,n]=size(e);
for i=1:m
for j=1:n
1 h, A& p2 T* T. |3 X8 eif (e(i,j)>=0)&&(e(i,j)<=50)
u=23;
q=23/3;
v(i,j)=e(i,j)-u;
2 A  A, v% {: D! k* [( ?0 LV(1)=2*q^2;
$ j/ _! @# l* M# [' m6 h% ]; welseif (e(i,j)>=51)&&(e(i,j)<=100)
( ^2 f# W+ {1 }u=83;
9 x3 m( F# A4 r# tq=47/3;
v(i,j)=e(i,j)-u;
. ?9 S  o2 A- UV(2)=2*q^2;
" C$ y+ a4 u! C- Qelse
u=117;
q=138/3;
v(i,j)=e(i,j)-u;
V(3)=2*q^2;
end
/ n* F+ G/ P' x4 |) K7 i; Nend
9 D; U* j( o4 \6 Y6 Eend
v=im2double(v);
Q=v.^2;
for i=1:m
for j=1:n
$ v0 S8 `1 B$ Q; o& C; zif (e(i,j)>=0)&&(e(i,j)<=50)
H(i,j)=im2double(Q(i,j)/V(1));
. s) {5 m8 D7 A8 Y3 IN(i,j)=exp(-H(i,j));
. u( j; y* I! f6 ^elseif (e(i,j)>=51)&&(e(i,j)<=100)
H(i,j)=im2double(Q(i,j)/V(2));
* _/ f$ Y$ z" n: LN(i,j)=exp(-H(i,j));
, m- m# d7 J- B+ C( z6 c- v% oelse
H(i,j)=im2double(Q(i,j)/V(3));
4 k/ w  v: x- [  ^1 oN(i,j)=exp(-H(i,j));
# M6 ~# E! y/ ?4 }# nend
end
end
& g: q9 C5 u* b. G2 H1 Is=ones(224,210);
( A( J' p) l- i6 G2 G" ?l=s*26;
l(2:223,2:209)=e;
* W. _3 x, k- D7 t, R6 io=zeros(224,210);
o(2:223,2:209)=N;
2 h2 n# w  s3 p4 T. A[m,n]=find(l<=25);
q=size(m);
+ N) B) R2 \% i: m' C$ i) E7 w9 o+ Sfor i=1:q
- e0 ~# O0 S7 U( ?/ G  ?- 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);
$ H/ p  t4 c0 s6 ^& l( Py(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);
7 c7 y0 K' E/ z; Z" q8 W0 a/ r/ wl(m(i),n(i))=x(m(i),n(i))/y(m(i),n(i));
end
; Y& J6 @& d/ q8 R& Lh=l(2:223,2:209);
h=round(h);
6 E% G# S* j% C- F%figure,imshow(h,[]);
- J2 {7 v8 Q7 k! \: x2 f%title('第二次去噪效果');
b
e
" F2 ^+ |+ x( U6 z9 Uh
/ D. L$ \. w) C8 rend

. S* b  a( V0 j4 {9 j) C: D

sharkN

转转