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

jiulquq

我下面的MATLAB程序，我想输出三个结果，为啥结果只会出现一个，就是b,不会出现其他的结果呢？不知道是什么原因，请大神帮帮我，谢谢了。
3 V1 l  h& o2 {- c9 X" Ifunction[b,e,h]=lvboshibie(a)
a=imread('lena.bmp');
  Y2 u8 Z( T: ^) p- k' m' [%figure;imshow(a);
( n, v; ?' b8 b! @) S- r  k%title('原图像');
( d- O8 }! j1 M9 ^a=im2double(a);%将数变为double型数
* b2 \' _, g* y6 g- Rb=imnoise(a,'gaussian',0.015);
C=im2double(b);
% W0 n/ u6 s2 q8 b7 `5 S6 B+ @; ~8 sF=im2col(C,[3,3],'distinct');%将加躁图像分割成3*3的小窗口
[i,j]=size(F);
m=zeros(i,j);
n=size(i,j);
for i=1:9
  u3 A9 k& t5 x$ n6 W- i# v. F! ?" vfor j=1:5180
m(i,j)=F(i,j)-F(5,j);%用3*3的窗口中的中心值对窗口中的数据求出估计的受躁程度
m(5,j)=0.001;
n(i,j)=m(i,j)./F(5,j);
4 [( W% ~7 g" K( E1 D$ Vend
end
%为用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];
3 y* k$ y% J' C5 w9 {; N2 L; U$ c- g$ P7 IA=K;
for ii=1:5179;
& g; Y4 p# E+ BA=[A;K];
2 x, m& g; \: T0 B& h5 I3 yend
g=zeros(1,5180);
5 i8 G7 S' g. q# k8 f1 D. lF=[1:5180];
F=[g,g,F];
1 B8 ]1 x/ c! [* S# N% Z, bF=[F;F;F;F;F;F;F;F;F];
2 K# K  [. t7 @/ AF=reshape(F,[46620,3]);
% [& g. h3 X2 |& qc=A+F;
% s% g0 o& G, `( x% d0 _% nval=reshape(n,[46620,1]);
- r8 }% v7 G1 G! r3 P0 e  @A=accumarray(c,val);%使用此函数把原来用im2col变成9*5180的矩阵变为一个个3*3的小矩阵
D=reshape(A,[3,15540]);
%将矩阵变为原来的222*210的形式
3 R- v0 p: ~) q4 e9 n$ g3 CM = [];
# ^0 H( w/ u0 y; IB=[];
for i = 1:3:15538
temp=D(:,i:i+2);

M = [M; temp]; % 使矩阵变为15540*3
; N+ j2 T8 W# H; K- xend
) Y2 `1 ^( R; O% d- R7 n! @for j=1:222:15539
- G; t% O$ }3 i8 o: L
/ x2 [) a$ i) p' ^  [: G- A) }TEMP=M(j:j+221,;
B = [B, TEMP]; % 使矩阵变为220*210
end
%使矩阵为222*208
2 M* \4 W4 u; F7 D+ ]5 X9 C( r7 kB=B(1:222,1:208);
E=B+100;
1 {6 b) Q5 x$ t* H7 |0 |) Ud=100.*a;%求出估计的灰度时用数据
: V$ p/ ]4 c( }" I( n/ Tr=d./E;
e=im2uint8(r);
& R% G+ K! P1 d; B+ D[m,n]=size(e);
* z, J: s: X4 N- E; ^( ~for i=1:m
" W! T9 ^3 K3 y+ ?5 v3 G6 k$ ]for j=1:n
/ p" i, |& i' W1 W+ [if (e(i,j)>=0)&&(e(i,j)<=50)
; X/ Q" m' H) zu=23;
2 w3 o9 b% ?) n9 v+ ~! Q3 S& Nq=23/3;
0 k2 R9 I% y. b: [4 uv(i,j)=e(i,j)-u;
2 m; k, {" ]. r) A( \$ p! LV(1)=2*q^2;
elseif (e(i,j)>=51)&&(e(i,j)<=100)
u=83;
7 w6 t! |8 q1 Z, v# Pq=47/3;
2 Y$ W8 I+ S3 e% O: Q" Sv(i,j)=e(i,j)-u;
V(2)=2*q^2;
else
9 n8 n3 _5 ?9 t" Z! v5 ?2 k. ^u=117;
q=138/3;
v(i,j)=e(i,j)-u;
V(3)=2*q^2;
+ s5 T- ?( K5 u4 V4 ~( Tend
end
2 F8 f) X  y1 _3 h! ]* Xend
& `2 t2 e; H" r  f: wv=im2double(v);
Q=v.^2;
* k( B/ f* c. S8 H; `for i=1:m
for j=1:n
2 [4 l% ]4 I( d/ H1 Zif (e(i,j)>=0)&&(e(i,j)<=50)
' v3 l1 v: [" Z: |6 \. R+ K, RH(i,j)=im2double(Q(i,j)/V(1));
N(i,j)=exp(-H(i,j));
elseif (e(i,j)>=51)&&(e(i,j)<=100)
5 s- d  a. m" c. b+ T. AH(i,j)=im2double(Q(i,j)/V(2));
' I5 O- M$ i- H9 t2 c& C- hN(i,j)=exp(-H(i,j));
& P$ d: ^4 f8 Y( n; k& F. Felse
H(i,j)=im2double(Q(i,j)/V(3));
; K1 P0 o. |. v  K" n+ iN(i,j)=exp(-H(i,j));
end
1 x/ c0 D2 {$ a9 X2 B/ v2 Bend
end
s=ones(224,210);
l=s*26;
l(2:223,2:209)=e;
o=zeros(224,210);
o(2:223,2:209)=N;
7 _. ~" S5 e$ w; N[m,n]=find(l<=25);
q=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);
8 C( S9 ^, b9 {; R4 d1 fy(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
9 }! ?7 }3 @1 c$ j7 uh=l(2:223,2:209);
) f* x( }) F( c: @' S0 _h=round(h);
: ?5 m7 i" ~* {%figure,imshow(h,[]);
  i) B) |" D9 @* j- c0 M' b4 E%title('第二次去噪效果');
2 Z$ u' A8 R  eb
+ B4 e. W8 M- a8 z+ U& w) q) ~' X) ne
5 S5 a4 L$ {/ _& `4 ch
end

sharkN

转转