X0=zeros(1,19);
Y0=zeros(1,19);
Z0=zeros(1,19);
i=1;
for t=0:pi/3:6*pi
   X0(i)=cos(t);
   Y0(i)=sin(t);
   Z0(i)=t;
   i=i+1;
end
plot3(X0,Y0,Z0,'k-')
hold on

fftimage=fftshift(fft2(double(I))); % 对条纹图像 I 进行二维傅里叶变换
figure,plot(mat2gray(abs(fftimage(257,:)))); % 显示频谱图
arrinput = ginput; % 鼠标选取高低频之间的波谷位置，点一下，按回车
arrtemp = fftimage(257,:);
[FringeMax,FringeMaxIndex] = max(arrtemp(arrinput(1,1):512));
% 得到从选取点开始向右的最大值(波峰)，并记录位置
FringeMaxIndex = round(arrinput(1,1) + FringeMaxIndex - 1);
% 得到 1 级谱(基频)的波峰位置
maxvalue=max(fftimage(257,FringeMaxIndex:512));
% 取出 1 级谱的位置。257 表示第 257 行(中间行，谱的位置)
[m,n]=find(fftimage==maxvalue); % 取出 1 级谱最大值位置
fringe_num=abs(n-m); % 确定条纹个数
......................

%-------------------------------------------------------------------------
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};
......................

在MATLAB 命令中输入主函数
Syms t
f=t^4-t^2-2*t+5;
[x1,x2] = minJT(f,0,0.1)
进退法确定搜索区间函数文件minJT如下：
function [minx,maxx]= minJT(f,x0,h0,eps)
%目标函数: f ;
%初始点 : x0;
%初始步长: h0;
%精度 :esp;
%区间左端点 : minx;
%区间右端点: maxx;
......................

%***************  生成输入信号  *****%
    f = 2000                  %设定基频率为2000HZ
    Fs=44100;                 %设定采样频率Fs=44.1kHz
    N=256;                    %取的样本点数N
    n=0:N-1;                  %建立矩阵
    t=n/Fs;                   %采样时间t
    w = 2*pi*f;
    T = 2*pi/w;                   %求出正弦函数的周期
    x1 = sin(w*t);
    x2 =awgn(x1,5,'measured');      %加入SNR为5dB的高斯白噪声得到输入信号x2
    xn =x2-x1;                    %间接获得白噪声xn

......................

%***************  生成输入信号  *****%
    f = 2000                  %设定基频率为2000HZ
    Fs=44100;                 %设定采样频率Fs=44.1kHz
    N=256;                    %取的样本点数N
    n=0:N-1;                  %建立矩阵
    t=n/Fs;                   %采样时间t
    w = 2*pi*f;
    T = 2*pi/w;                   %求出正弦函数的周期
    x1 = sin(w*t);
    x2 =awgn(x1,5,'measured');      %加入SNR为5dB的高斯白噪声得到输入信号x2
    xn =x2-x1;                    %间接获得白噪声xn

......................