clear
[AB]=imread（‘p.jpg‘）;
X=imnoise（A‘gaussian‘0.001）;
subplot（131）imshow（AB）title（‘原圖‘）??
subplot（132）imshow（XB）title（‘加高斯白噪聲‘）
X=double（X）;
[iijjzz]=size（X）;
if?zz==3
for?i=1:ii
for?j=1:jj
?X（ij）=0.3*X（ij1）+0.59*X（ij2）+0.11*X（ij