有关chirp信号的frft程序

有关chirp信号的FRFT，包括FRFT卷积，FRFT乘法，以及chirp信号的参数估计等

%FRFT
%50次统计平均

clc
clear all
close all
%常量设置
f0=2000;
B=500:200:2500;
fs=4*f0;

ts=0.5;
t=-1*ts/2:1/fs:ts/2;
alfa=pi/6;
alfa11=30;


t=-1*ts/2:1/fs:ts/2;
S=sqrt（ts/fs）;
dx=sqrt（ts*fs）;
tt=-1*dx/2:1/dx:dx/2;%归一化之后的数据
point=ts*fs+1;


snr=5;

%滤波器
% b=fir1（256[（f0-B/2）/（fs/2）（f0+B/2）/（fs/2）]）;
for k=1:length（B）%k为snr;
    k
    for k1=1:50%k1为统计平均次数
%造信号
        signal=exp（i*2*pi*（f0）*tt*S+i*pi*（B（k）/ts）*tt.*tt*S*S）;
        noise=normrnd（011point）;
        nx=normrnd（011point）;
        ny=normrnd（011point）;

        %加信噪比
        %信噪比db/Hz
        b1=10^（snr/10）;  
        a1=std（signal）;
        a2=std（noise）;
        n=sqrt（b1*a2*B（k）/a1/fs）;

        signal1=n*signal;
        p=signal1+noise;

        %通过滤波器
%         p1=conv（b

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2668  2010-10-26 11:01  BFRFT.m
     文件        2373  2010-06-19 09:37  chirpallfft.m
     文件        2237  2010-06-02 19:13  chirpcexiang.m
     文件        1182  2010-06-04 11:26  Chirpexample.m
     文件        2135  2010-07-03 12:33  chirpFrft.m
     文件        2527  2010-05-11 11:38  chirpFrfthuadongchuang.m
     文件         973  2010-12-28 15:09  corror.m
     文件        1959  2010-05-20 20:59  CWcexiang.m
     文件        1560  2010-06-06 15:24  cWchirp.m
     文件        2541  2013-06-06 22:05  Disfrft.m
     文件         141  2010-04-19 16:15  fconv.m
     文件        1620  2010-04-19 16:40  frft.m
     文件         192  2010-04-19 16:21  interp.m
     文件        2727  2010-04-25 10:42  jiaoduguji.m
     文件         310  2010-05-19 09:42  matlabalfa2.mat
     文件         316  2010-05-19 09:41  matlaberror2.mat
     文件        2643  2010-07-03 22:09  PJSQQfft.m
     文件        1276  2011-01-05 16:57  pujisnr.m
     文件         715  2010-04-23 09:22  RectExample.m
     文件         977  2010-12-30 16:05  twodom.m