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

jiulquq

我下面的MATLAB程序，我想输出三个结果，为啥结果只会出现一个，就是b,不会出现其他的结果呢？不知道是什么原因，请大神帮帮我，谢谢了。
function[b,e,h]=lvboshibie(a)
a=imread('lena.bmp');
%figure;imshow(a);
" I/ [6 }; ?, C%title('原图像');
a=im2double(a);%将数变为double型数
  B6 N" S( u- c; ]2 ?b=imnoise(a,'gaussian',0.015);
C=im2double(b);
! g4 I* R& D. r! J4 bF=im2col(C,[3,3],'distinct');%将加躁图像分割成3*3的小窗口
7 e# K( o% n( A; v[i,j]=size(F);
$ P) B( ?- i* ?, ^m=zeros(i,j);
n=size(i,j);
" _% S- `' w! ]2 n6 V' y( \for i=1:9
for j=1:5180
m(i,j)=F(i,j)-F(5,j);%用3*3的窗口中的中心值对窗口中的数据求出估计的受躁程度
m(5,j)=0.001;
2 o0 g; e/ {# y' a# f9 e& F7 jn(i,j)=m(i,j)./F(5,j);
end
! S. I. o8 i' H9 ~end
* E5 C  H6 T2 N  ]0 v0 w%为用accumarray函数把数据变回原来的排列方式做准备
%制造c矩阵
$ J9 i4 [2 G4 V! }3 _# I. iK=[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;
: q) L: Z* t6 \8 A1 E2 Tfor ii=1:5179;
$ a$ S& i& |5 q# pA=[A;K];
end
2 o8 i% j$ a. h! E3 {) m5 y9 Eg=zeros(1,5180);
. L0 W% n7 M, |: ]$ [* A! z4 \F=[1:5180];
F=[g,g,F];
F=[F;F;F;F;F;F;F;F;F];
. X8 t  F" Y  {/ ^5 M9 c  c5 @+ YF=reshape(F,[46620,3]);
c=A+F;
7 s! e! R, `  I% ^3 u# z6 Q; hval=reshape(n,[46620,1]);
8 i& S4 ~8 C+ DA=accumarray(c,val);%使用此函数把原来用im2col变成9*5180的矩阵变为一个个3*3的小矩阵
, E" ?5 t' |4 b$ M/ u7 hD=reshape(A,[3,15540]);
' k5 Q3 W+ Y: W1 _%将矩阵变为原来的222*210的形式
+ _6 P& T4 T7 \; M4 T& dM = [];
B=[];
for i = 1:3:15538
+ f( k+ ?. c+ ttemp=D(:,i:i+2);

M = [M; temp]; % 使矩阵变为15540*3
end
for j=1:222:15539
# `) A  r( c; V% l2 A" Z4 O7 r- `
  k. r8 _; {- d$ z' I+ d5 QTEMP=M(j:j+221,;
B = [B, TEMP]; % 使矩阵变为220*210
end
%使矩阵为222*208
B=B(1:222,1:208);
, K" w8 [$ }8 i( r2 UE=B+100;
d=100.*a;%求出估计的灰度时用数据
r=d./E;
e=im2uint8(r);
[m,n]=size(e);
for i=1:m
& b1 A7 T+ M- f5 w' S- _6 _for j=1:n
if (e(i,j)>=0)&&(e(i,j)<=50)
- u/ x- S9 h; I. I- u% au=23;
q=23/3;
v(i,j)=e(i,j)-u;
V(1)=2*q^2;
elseif (e(i,j)>=51)&&(e(i,j)<=100)
5 e& q7 Z' S) e* q8 O: g/ p* V$ bu=83;
2 P1 s, n. P9 l( j, u& Gq=47/3;
v(i,j)=e(i,j)-u;
V(2)=2*q^2;
9 T" n- I8 h, r+ ?% }else
, ?9 Z6 L& A) s% M) h" j7 E: Eu=117;
: Z3 r) P5 a2 ~/ X% jq=138/3;
6 `" J. g, L- ?4 Z( b+ p3 h; x8 Cv(i,j)=e(i,j)-u;
V(3)=2*q^2;
4 h$ y3 G3 `7 `6 yend
end
end
v=im2double(v);
Q=v.^2;
7 b) ?7 c" H8 R  Pfor i=1:m
for j=1:n
- g' M- K8 N: W7 Qif (e(i,j)>=0)&&(e(i,j)<=50)
H(i,j)=im2double(Q(i,j)/V(1));
N(i,j)=exp(-H(i,j));
elseif (e(i,j)>=51)&&(e(i,j)<=100)
H(i,j)=im2double(Q(i,j)/V(2));
3 E9 s; g1 M- W0 b8 V! ]N(i,j)=exp(-H(i,j));
else
H(i,j)=im2double(Q(i,j)/V(3));
( a: r7 j  f! `3 JN(i,j)=exp(-H(i,j));
end
end
4 s" T! i! `! N7 {4 S; v& nend
8 i5 ^) H7 E5 D" Ns=ones(224,210);
. x8 N4 J0 R+ ]7 ll=s*26;
. o* h+ A+ j* n4 v  rl(2:223,2:209)=e;
o=zeros(224,210);
2 f' U4 s6 {  Vo(2:223,2:209)=N;
7 T, E9 |/ s' g1 V" O. j0 s[m,n]=find(l<=25);
6 G& U, S" B# Rq=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));
end
" b* i: O4 j' ^h=l(2:223,2:209);
h=round(h);
%figure,imshow(h,[]);
%title('第二次去噪效果');
3 P7 Y' u9 _) K5 pb
' l& U9 |! A" T/ T7 s  s$ ?! Ze
2 |6 c$ B# Z; M- z3 A% oh
4 t" R+ }" ^! K. u# O" Xend

- g& a1 ]0 H& D$ s3 O- G% s+ Y

sharkN

转转