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

jiulquq

我下面的MATLAB程序，我想输出三个结果，为啥结果只会出现一个，就是b,不会出现其他的结果呢？不知道是什么原因，请大神帮帮我，谢谢了。
function[b,e,h]=lvboshibie(a)
4 `" v6 K& n. N7 }9 E0 ^a=imread('lena.bmp');
%figure;imshow(a);
%title('原图像');
, h/ g: l' Q' i5 ?a=im2double(a);%将数变为double型数
b=imnoise(a,'gaussian',0.015);
C=im2double(b);
- ]* g/ v/ u5 {5 v, }' F/ |F=im2col(C,[3,3],'distinct');%将加躁图像分割成3*3的小窗口
- o% h; d- ~& o: c[i,j]=size(F);
m=zeros(i,j);
n=size(i,j);
for i=1:9
for j=1:5180
7 _0 {+ O! l2 u% x7 Om(i,j)=F(i,j)-F(5,j);%用3*3的窗口中的中心值对窗口中的数据求出估计的受躁程度
m(5,j)=0.001;
' b& }& U1 V+ ^4 n) w: Rn(i,j)=m(i,j)./F(5,j);
* `" R& V7 X4 v' ^end
end
4 N' H& {! J. K# G3 s%为用accumarray函数把数据变回原来的排列方式做准备
%制造c矩阵
2 Z1 s8 o' U/ p* F* DK=[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;
1 O0 A8 W) {2 s) N& ?for ii=1:5179;
A=[A;K];
end
+ w+ |5 ^3 Z# q7 ^/ r! f, bg=zeros(1,5180);
: S0 _5 \" {+ `F=[1:5180];
F=[g,g,F];
3 ^1 H% S3 A; y4 n. kF=[F;F;F;F;F;F;F;F;F];
F=reshape(F,[46620,3]);
c=A+F;
7 V/ n; V# K6 \& X; g) b* }- Uval=reshape(n,[46620,1]);
' P2 u: @! j* [/ ~: X) @0 aA=accumarray(c,val);%使用此函数把原来用im2col变成9*5180的矩阵变为一个个3*3的小矩阵
/ k1 ?2 c5 w- T) b8 g4 xD=reshape(A,[3,15540]);
; y0 t4 E7 B" V% w6 h3 z( ?3 l%将矩阵变为原来的222*210的形式
$ S- r- A9 `  ?* fM = [];
B=[];
- h# z* S" h) i$ Vfor i = 1:3:15538
5 [8 Q+ R, ?: {& ]temp=D(:,i:i+2);

M = [M; temp]; % 使矩阵变为15540*3
end
, V1 V; U! E7 P" tfor j=1:222:15539

4 k3 x7 C+ ^. a0 X8 A% w3 c9 ZTEMP=M(j:j+221,;
& w" y2 B4 a/ W5 Q* gB = [B, TEMP]; % 使矩阵变为220*210
. r0 h0 O' O. D. n$ _1 g& ~end
%使矩阵为222*208
9 ^& ^& c! J$ wB=B(1:222,1:208);
E=B+100;
5 b* C  A' @, \( X( k0 t* q+ Id=100.*a;%求出估计的灰度时用数据
r=d./E;
e=im2uint8(r);
[m,n]=size(e);
for i=1:m
& U% A. j$ ?/ h2 ]: hfor j=1:n
; Z! o5 D( M6 bif (e(i,j)>=0)&&(e(i,j)<=50)
& I; e/ O0 ^1 ]; U7 R+ P/ Zu=23;
* W' t1 l& q9 f. Y4 \6 f0 m, ^; p9 ^q=23/3;
v(i,j)=e(i,j)-u;
V(1)=2*q^2;
elseif (e(i,j)>=51)&&(e(i,j)<=100)
7 `* M: M! P! E* E3 Iu=83;
q=47/3;
v(i,j)=e(i,j)-u;
+ C. f$ S4 O7 L! BV(2)=2*q^2;
else
. T  u% x2 J0 w5 |8 j: R7 S4 cu=117;
9 b4 W! R4 y# x( F, ^q=138/3;
; _# n% t" f) H7 q; h* @v(i,j)=e(i,j)-u;
V(3)=2*q^2;
end
% T& c" F0 v2 \9 W6 M8 X/ i; send
end
v=im2double(v);
* ^0 \7 _+ r+ P3 K3 MQ=v.^2;
for i=1:m
+ V# ]* Y+ o% ]/ Z  x- Ffor j=1:n
if (e(i,j)>=0)&&(e(i,j)<=50)
H(i,j)=im2double(Q(i,j)/V(1));
% A2 h) ^3 L$ ?N(i,j)=exp(-H(i,j));
elseif (e(i,j)>=51)&&(e(i,j)<=100)
) t) R0 ^# |0 BH(i,j)=im2double(Q(i,j)/V(2));
N(i,j)=exp(-H(i,j));
5 K1 `2 G* \2 @" r8 N% qelse
3 r+ S" h$ X; y+ d  |H(i,j)=im2double(Q(i,j)/V(3));
N(i,j)=exp(-H(i,j));
end
end
8 a- O3 H. G4 ?* v0 yend
s=ones(224,210);
8 f6 x7 O9 e5 s$ q; X; F. |l=s*26;
$ Z2 U* y4 c3 b# Z6 vl(2:223,2:209)=e;
o=zeros(224,210);
o(2:223,2:209)=N;
$ O3 X6 c# u9 L6 D[m,n]=find(l<=25);
& }, I0 P5 g- q* j+ ?q=size(m);
for i=1:q
/ ~- f9 o: v8 z  }' S- zx(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);
( N+ u. s' E8 k$ p2 vy(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);
" U; J0 Y8 P/ E' Z, a/ O- Hl(m(i),n(i))=x(m(i),n(i))/y(m(i),n(i));
end
h=l(2:223,2:209);
h=round(h);
%figure,imshow(h,[]);
4 E& E7 h9 e! o: M- N& y5 b: j%title('第二次去噪效果');
2 o( H- }8 x7 x4 z2 C  @) Vb
e
& S; J2 b$ U  ^7 z: q: f* Ch
end

. E- |+ P! Z# a

sharkN

转转