%-------------------------------------------------------------------------
clear all
x=imread('test_image_9.jpg');  %读取原始图像
figure(1);subplot(1,4,1);imshow(x);title('原图像1');
y=rgb2gray(x);
subplot(1,4,2);imshow(y);title('图像1的灰度图');           % 图1灰度图
u1=imnoise(y,'salt & pepper',0.07);
subplot(1,4,3);imshow(u1);title('图像1加噪声图');     %给图1加椒盐噪声
zz=medfilt2(u1,[3 3]);z=medfilt2(zz,[5 5]);%(2次中值滤波)
subplot(1,4,4);imshow(z);title('图像1中值滤波');  %图1中值滤波（3*3窗口）
% figure(6);subplot(1,2,1);imshow(zz);title('一次滤波');
% subplot(1,2,2);imshow(z);title('二次滤波');
......................

%给图像加扰动，扰动方式以zernike多项式的形式加入，Image为源图像，Co为要加的zernike系数，为行向量，
%输出分别为扰动后的图像，扰动PV值，扰动RMS值
function [OUT,dis_pv,dis_rms]=Image_disturb(Image,Co)
Image=imread('1.bmp');
M=size(Image,1);N=size(Image,2);
%--------所加扰动，用zernike多项式表示--------------------
Co=[1,1,1,0.02,0.02,0.02,0.02,0.02,0.02];
N_Zer=size(Co,2);%所用zernike阶数
disturb=zeros(M,N);
Zer=zernike_poly(M,N);%1-36阶
for q=1:N_Zer 
    disturb=disturb+Co(1,q)*Zer{q};
......................