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

jiulquq

我下面的MATLAB程序，我想输出三个结果，为啥结果只会出现一个，就是b,不会出现其他的结果呢？不知道是什么原因，请大神帮帮我，谢谢了。
9 G2 u9 E6 P) b2 W6 B, U- rfunction[b,e,h]=lvboshibie(a)
a=imread('lena.bmp');
% |- ]5 p, k  t( A/ T%figure;imshow(a);
%title('原图像');
a=im2double(a);%将数变为double型数
b=imnoise(a,'gaussian',0.015);
C=im2double(b);
F=im2col(C,[3,3],'distinct');%将加躁图像分割成3*3的小窗口
[i,j]=size(F);
4 o& b( g% n' y& x  f! tm=zeros(i,j);
' u. S; o  T; E0 kn=size(i,j);
for i=1:9
for j=1:5180
) E* v6 \. B, X) h- ?m(i,j)=F(i,j)-F(5,j);%用3*3的窗口中的中心值对窗口中的数据求出估计的受躁程度
- f- |# f) T7 f, R' Im(5,j)=0.001;
n(i,j)=m(i,j)./F(5,j);
end
end
1 G5 H6 D" A/ ]  H) b* a%为用accumarray函数把数据变回原来的排列方式做准备
%制造c矩阵
K=[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;
+ O9 p6 m" z( {4 n7 x( h- Xfor ii=1:5179;
! w2 Z) h9 A# l. |' x) ^- ^A=[A;K];
/ h3 J- t: S5 x7 d# M- P+ Jend
g=zeros(1,5180);
F=[1:5180];
F=[g,g,F];
" ~; Y: V$ p9 p* M& z# Y3 pF=[F;F;F;F;F;F;F;F;F];
F=reshape(F,[46620,3]);
+ x9 x& o% H4 v1 x& j: o: Pc=A+F;
val=reshape(n,[46620,1]);
A=accumarray(c,val);%使用此函数把原来用im2col变成9*5180的矩阵变为一个个3*3的小矩阵
: A4 D* t; A$ u- tD=reshape(A,[3,15540]);
" B* R( w- T" D2 g1 i$ I%将矩阵变为原来的222*210的形式
M = [];
B=[];
/ a+ u$ u# D( w6 [% _3 P6 j) Qfor i = 1:3:15538
7 g7 e# V" S5 K8 x# j1 stemp=D(:,i:i+2);

M = [M; temp]; % 使矩阵变为15540*3
5 p3 a$ A) L3 o1 B" H& q! Q9 iend
2 I% y4 O5 [% _# b# o! h) c# N1 ?/ hfor j=1:222:15539
' U$ R+ K2 g6 L# C- ~0 U- Q
! m) }9 |3 e6 L) ^3 e% D' z* ZTEMP=M(j:j+221,;
& X8 P) N3 s' U  v1 ^B = [B, TEMP]; % 使矩阵变为220*210
2 y: ?  Q- H7 F, ~; a" ?end
%使矩阵为222*208
$ D% u/ A' X. o, [B=B(1:222,1:208);
E=B+100;
6 _, Q. O- x* B3 E. `( a! N) Id=100.*a;%求出估计的灰度时用数据
r=d./E;
" K, {& o8 ?/ \  v8 }' C9 ue=im2uint8(r);
[m,n]=size(e);
for i=1:m
for j=1:n
if (e(i,j)>=0)&&(e(i,j)<=50)
u=23;
0 s9 y  }7 x( X; h+ b$ i" g$ iq=23/3;
v(i,j)=e(i,j)-u;
' F5 n: l" ^  gV(1)=2*q^2;
* _$ \. A( {6 _0 felseif (e(i,j)>=51)&&(e(i,j)<=100)
' s8 c$ M' F( y  x; U% F' O! c6 t( ~8 ku=83;
q=47/3;
v(i,j)=e(i,j)-u;
V(2)=2*q^2;
5 k4 M8 ]7 w6 V- P3 g7 S/ S- t2 B3 xelse
8 Q2 j5 ]/ H8 W- ku=117;
& S0 F+ N7 P6 Z6 P6 L& x' Gq=138/3;
v(i,j)=e(i,j)-u;
V(3)=2*q^2;
1 ?' k5 z! g- ?- D9 s1 p3 `end
( P+ k7 v& @! A' d) D2 p: u& kend
0 M; E, A7 E5 D# y7 s) n! zend
v=im2double(v);
/ n- O& @' U2 ^! L5 V( e$ S* ^6 P& d$ Q$ _Q=v.^2;
3 s; w$ d# }; k7 p3 F( G: Nfor i=1:m
for j=1:n
% J$ K+ }  z6 {- kif (e(i,j)>=0)&&(e(i,j)<=50)
H(i,j)=im2double(Q(i,j)/V(1));
% y3 f: N, s! aN(i,j)=exp(-H(i,j));
elseif (e(i,j)>=51)&&(e(i,j)<=100)
2 s" t: @) U/ [/ j6 E, Y1 iH(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));
! Y+ N1 s$ q1 iend
/ \' r: W; P5 k! P5 x% w) ~0 u, nend
end
6 ?% o! O) E5 i: \5 W4 G2 C% L- Ts=ones(224,210);
, M# p5 s2 f1 Z0 Q" p- J* Sl=s*26;
l(2:223,2:209)=e;
3 y; ^" f! F+ z6 b3 eo=zeros(224,210);
o(2:223,2:209)=N;
: a9 V& U% H5 D9 @* S: k[m,n]=find(l<=25);
q=size(m);
3 w: v/ I+ Y% i4 ?8 ?4 S7 Nfor i=1:q
6 |# z# J& H0 D! X$ |. T; Ux(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));
* y& m* Q# X& {/ r& qend
3 F; h1 j/ q$ i7 o# N1 H# E& d. I+ Th=l(2:223,2:209);
5 r; G- K. |: f( Z4 C4 Ph=round(h);
%figure,imshow(h,[]);
; e4 C4 m' p4 u; V6 u%title('第二次去噪效果');
3 N! l( U' D, `/ N' \/ w7 cb
& p; [+ c9 Z- He
h
0 f( h+ \  L$ V& S1 d6 ~3 o* {end

* ]3 I5 `, u/ V9 W; ~

sharkN

转转