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

jiulquq

我下面的MATLAB程序，我想输出三个结果，为啥结果只会出现一个，就是b,不会出现其他的结果呢？不知道是什么原因，请大神帮帮我，谢谢了。
1 k% p( s. P, ?$ Ufunction[b,e,h]=lvboshibie(a)
a=imread('lena.bmp');
0 R% \* f2 S3 D%figure;imshow(a);
5 X9 S& h; g# x2 y* ^! e%title('原图像');
a=im2double(a);%将数变为double型数
) M* K, |2 ^# Z5 tb=imnoise(a,'gaussian',0.015);
' H: m4 L) [; K6 v0 x) D: \C=im2double(b);
F=im2col(C,[3,3],'distinct');%将加躁图像分割成3*3的小窗口
" w, v4 y5 w+ x+ d" @0 P" q* B[i,j]=size(F);
$ Y: Y/ a* W8 o- C( \" {$ j$ em=zeros(i,j);
n=size(i,j);
for i=1:9
for j=1:5180
) A# V: @" _5 z6 i, Bm(i,j)=F(i,j)-F(5,j);%用3*3的窗口中的中心值对窗口中的数据求出估计的受躁程度
$ G9 j8 Z  x* s% C- Z' b2 k( r* lm(5,j)=0.001;
! p+ n6 ]8 d' Q) K3 C6 An(i,j)=m(i,j)./F(5,j);
6 X8 g9 T1 X4 ]7 y) ?4 nend
9 Q4 c7 I4 x: O9 r6 q3 d4 eend
2 q* B( K6 |9 o%为用accumarray函数把数据变回原来的排列方式做准备
( ^1 ^  ]1 y+ p; Z9 \% e4 t) s%制造c矩阵
1 D$ ^, b, M# U5 H  I8 F3 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];
A=K;
  ^' d+ Y& H: R# o1 hfor ii=1:5179;
A=[A;K];
end
9 j, U/ g: s4 A5 c4 qg=zeros(1,5180);
F=[1:5180];
% }" [5 N$ @- N% R* ?# e! a& ~% EF=[g,g,F];
& g( p, t% N* f5 [5 z$ ZF=[F;F;F;F;F;F;F;F;F];
F=reshape(F,[46620,3]);
c=A+F;
val=reshape(n,[46620,1]);
A=accumarray(c,val);%使用此函数把原来用im2col变成9*5180的矩阵变为一个个3*3的小矩阵
D=reshape(A,[3,15540]);
4 B. A7 ~- l1 ^: y  }%将矩阵变为原来的222*210的形式
; x4 s5 L& l5 DM = [];
B=[];
" ~7 @9 Q- ], K  T3 u5 D8 wfor i = 1:3:15538
temp=D(:,i:i+2);

M = [M; temp]; % 使矩阵变为15540*3
" \, K7 E! ^- h! G% W* f) F7 V% Oend
for j=1:222:15539

TEMP=M(j:j+221,;
1 G" }) A5 C0 e$ F- I( h! sB = [B, TEMP]; % 使矩阵变为220*210
end
, R; {4 V% A7 ]6 q+ J! ]%使矩阵为222*208
B=B(1:222,1:208);
7 u7 i$ \! G; L. x- k1 GE=B+100;
d=100.*a;%求出估计的灰度时用数据
r=d./E;
e=im2uint8(r);
4 V% u8 p, c  [6 h[m,n]=size(e);
for i=1:m
for j=1:n
if (e(i,j)>=0)&&(e(i,j)<=50)
u=23;
* [+ m7 b! J2 Tq=23/3;
v(i,j)=e(i,j)-u;
) B0 W. u$ S  c2 {# @  i3 lV(1)=2*q^2;
6 a# f& M8 |- ?; G9 g5 lelseif (e(i,j)>=51)&&(e(i,j)<=100)
u=83;
q=47/3;
: s: @1 [- S. r$ \) o: `3 mv(i,j)=e(i,j)-u;
V(2)=2*q^2;
4 b: u: ]2 i" d5 o5 U+ Jelse
' T) z( u( ]/ [5 w( M- _u=117;
$ K- c- {% K& c, e0 D" F# Qq=138/3;
. i/ N6 W# j& @) r6 Pv(i,j)=e(i,j)-u;
V(3)=2*q^2;
7 W: K2 r% s% Q9 V, L3 s' S+ |end
7 u: O, y) E  r& i. N9 z) g$ dend
end
. Q8 {( N: e; y( F' l. T* w; w. Z' }v=im2double(v);
# g& s7 Q( k* `$ u/ V4 S# XQ=v.^2;
for i=1:m
for j=1:n
if (e(i,j)>=0)&&(e(i,j)<=50)
H(i,j)=im2double(Q(i,j)/V(1));
* u) F5 ^" k  o, ^2 t" V( `1 h& B' M% xN(i,j)=exp(-H(i,j));
elseif (e(i,j)>=51)&&(e(i,j)<=100)
$ N3 t$ c0 |8 {: t: j, mH(i,j)=im2double(Q(i,j)/V(2));
N(i,j)=exp(-H(i,j));
; s, |: g# y4 c. ]% M3 ^else
H(i,j)=im2double(Q(i,j)/V(3));
N(i,j)=exp(-H(i,j));
0 q4 g& D! j( q: [( ]" qend
1 I9 m2 p  I# }# ?end
end
s=ones(224,210);
; M: F# K& I! N2 k! Q- G+ bl=s*26;
l(2:223,2:209)=e;
1 I( b" ], ?% p: }% d# o$ Go=zeros(224,210);
! J) T4 N% F) N! K& k1 @% D2 Yo(2:223,2:209)=N;
[m,n]=find(l<=25);
q=size(m);
4 D/ p, _+ B$ Q; yfor i=1:q
" ]0 _1 V0 u& \* U! M* ex(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);
, u8 \( O( m8 j  C! }+ P+ ]l(m(i),n(i))=x(m(i),n(i))/y(m(i),n(i));
4 ~' B, l7 {# ^! D: ]( w3 pend
h=l(2:223,2:209);
, a; A; ^$ L8 Z* @1 \9 N* oh=round(h);
%figure,imshow(h,[]);
6 U, y* V, |; F% _. c3 v%title('第二次去噪效果');
0 r2 x, Q. l0 {5 z) m; S' _, |. E5 I! Ib
, o' S7 X, a" R. B" K7 ce
h
( M, `. {9 ]/ t$ G" X& }. xend

sharkN

转转