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

jiulquq

我下面的MATLAB程序，我想输出三个结果，为啥结果只会出现一个，就是b,不会出现其他的结果呢？不知道是什么原因，请大神帮帮我，谢谢了。
+ I. B6 E6 N0 @+ u+ {/ @2 w, ?5 Bfunction[b,e,h]=lvboshibie(a)
a=imread('lena.bmp');
%figure;imshow(a);
1 x* D3 k! ]$ }  J%title('原图像');
a=im2double(a);%将数变为double型数
b=imnoise(a,'gaussian',0.015);
) Y0 y' t; e& U/ k2 ^C=im2double(b);
6 e  W% O7 v4 [1 X' hF=im2col(C,[3,3],'distinct');%将加躁图像分割成3*3的小窗口
. y) Z, V1 R0 @) Z. r: L/ \[i,j]=size(F);
m=zeros(i,j);
n=size(i,j);
for i=1:9
% s( ?  s4 X3 R2 f+ B- Kfor 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);
end
end
$ \6 @* x* S$ Z8 y4 y. X%为用accumarray函数把数据变回原来的排列方式做准备
" W* }- ^0 ]# P" \3 h: w8 S9 ?0 s+ J%制造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];
( |2 d3 J4 {' z1 \/ B2 L* q6 _A=K;
& [( @1 r/ Z4 T5 J' j. l/ ]for ii=1:5179;
8 L, R: |5 \# N% g5 s) W- \) I! PA=[A;K];
5 ~; R. m' Y1 P9 G2 ]" Aend
/ y- k, a& S% _3 z+ ^g=zeros(1,5180);
F=[1:5180];
. d; i/ o$ x& k2 u% c$ c; j1 YF=[g,g,F];
: d4 _9 h& I8 @9 x3 OF=[F;F;F;F;F;F;F;F;F];
F=reshape(F,[46620,3]);
' [$ _. l3 g" p4 @: s* _c=A+F;
3 p+ r$ k2 R9 o5 U" Q4 _0 }% Pval=reshape(n,[46620,1]);
2 V3 z5 n. B! X) {9 E5 m, DA=accumarray(c,val);%使用此函数把原来用im2col变成9*5180的矩阵变为一个个3*3的小矩阵
D=reshape(A,[3,15540]);
%将矩阵变为原来的222*210的形式
M = [];
& D8 b  R5 a- k9 @6 lB=[];
for i = 1:3:15538
. J' J" u5 E0 U! O- ftemp=D(:,i:i+2);

M = [M; temp]; % 使矩阵变为15540*3
1 p7 `; B* p5 B4 X9 e% Zend
for j=1:222:15539
6 y% ~9 r- g- @6 O/ u1 x
TEMP=M(j:j+221,;
B = [B, TEMP]; % 使矩阵变为220*210
8 [# J7 W9 U# T, w1 Send
! f" O. r, b5 D4 m5 |) F5 D5 M%使矩阵为222*208
B=B(1:222,1:208);
+ w9 |4 \5 W1 E! q" `/ }E=B+100;
d=100.*a;%求出估计的灰度时用数据
r=d./E;
e=im2uint8(r);
[m,n]=size(e);
( i! j( V. w9 b8 K9 m9 a1 Rfor i=1:m
# S/ t; v/ t& `$ [+ [7 O2 Nfor j=1:n
) h' h7 C) w5 t: S# Vif (e(i,j)>=0)&&(e(i,j)<=50)
u=23;
$ D* Z( B6 w1 q* G8 y8 ~  xq=23/3;
4 n1 \2 Q$ G, Rv(i,j)=e(i,j)-u;
V(1)=2*q^2;
elseif (e(i,j)>=51)&&(e(i,j)<=100)
4 B( i; o# T* Qu=83;
q=47/3;
v(i,j)=e(i,j)-u;
V(2)=2*q^2;
else
: o" \" ]$ Y( p' W/ [* g5 Nu=117;
3 }( J  Y7 U2 m7 B" Tq=138/3;
% ^5 ~( y' |: ev(i,j)=e(i,j)-u;
( E% U, C- H" ~1 ]V(3)=2*q^2;
end
end
6 w; Y( ~* d  n) C6 Z! Mend
v=im2double(v);
Q=v.^2;
- F) \3 E- A' [* ~; K# pfor i=1:m
for j=1:n
if (e(i,j)>=0)&&(e(i,j)<=50)
# j3 z8 x' r/ ]+ d9 D+ FH(i,j)=im2double(Q(i,j)/V(1));
8 `2 Q7 {( _1 P- N# c4 \/ U/ pN(i,j)=exp(-H(i,j));
1 S( `$ L$ A% ]1 Eelseif (e(i,j)>=51)&&(e(i,j)<=100)
+ q9 k. k1 |8 X# V. x, O6 HH(i,j)=im2double(Q(i,j)/V(2));
N(i,j)=exp(-H(i,j));
else
' S$ s5 S8 S4 j8 hH(i,j)=im2double(Q(i,j)/V(3));
N(i,j)=exp(-H(i,j));
" K  q* T5 L2 C3 r! u" x5 {9 Gend
7 s2 y4 l* X* W& F- |8 v9 Oend
end
2 V: J- U4 q2 m! e: A- hs=ones(224,210);
l=s*26;
. a; d# `0 w  il(2:223,2:209)=e;
" ~$ y3 y& j; u6 to=zeros(224,210);
o(2:223,2:209)=N;
[m,n]=find(l<=25);
q=size(m);
for i=1:q
0 y% r. S/ i3 jx(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
h=l(2:223,2:209);
h=round(h);
$ \' e8 ^, A3 F1 k# m. S%figure,imshow(h,[]);
%title('第二次去噪效果');
b
e
h
end

6 H" [. g1 }2 ~$ n8 u" `$ w  Q

sharkN

转转