hmar94.rar

hmar94.rar

function Fresnel_demo（）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 此文件为菲涅尔积分演示程序（菲涅尔积分是求线性调频信号频谱的必要手段）
% 输入自变量x调用菲涅尔积分函数，分别获得实部和虚部
% 输出三个figure演示菲涅尔积分的空间特征
% 参考文献：国鹏.基于单脉冲技术的雷达三维成像[D].北京:北京航空航天大学.2008.pp35-36
% 樊文贵/北京航空航天大学电子信息工程学院
% 5/12/2014
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x=0:0.02:10;            %输入自变量数据
[CS]=Fresnel_init（x）;  %调用菲涅尔积分

%% 分别获得积分实部和虚部
figure（1）;
plot（xC‘b‘‘linewidth‘2）;hold on;
plot（xS‘r‘‘linewidth‘2）;hold off;
title（‘C（x）与S（X）震荡曲线‘）;
legend（‘C（x）‘‘S（x）‘）;
grid on;

%% 以x为自变量获得二维菲涅尔积分曲线
figure（2）;
plot（CS‘linewidth‘2）;
title（‘二维菲涅尔积分曲线‘）;
grid on;

%% 以x为自变量获得空间菲涅尔积分曲线
figure（3）;
n=zeros（length（x））;
plot3（xCS‘b‘‘linewidth‘2）;hold on;
%%绘三维图形比较费时间
% plot3（nCS‘k‘‘linewidth‘2）;hold on;
% plot3（xnS‘r‘‘linewidth‘2）;hold on;
% plot3（xCn‘y‘‘linewidth‘2）;hold off;
grid on;xlabel（‘x‘）;ylabel（‘C（x）‘）;zlabel（‘S（x）‘）;
title（‘空间菲涅尔积分曲线‘）;
end

function [CS] = Fresnel_init（t）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 菲涅尔积分函数
% 利用符号表达式求得菲涅尔积分表示
% 利用替换函数将符号替换为数值并利用vpa函数求得对应数值解
% 输入 x 积分上限
% 输出 C 积分值实部
%      S 积分值虚部
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
syms x;                                  %设置表达式变量
freint=int（exp（1i*pi*x.^2/2）x0x）;     %求得0~x的菲涅尔积分
fr=vpa（subs（freint‘x‘t）8）;            %利用VPA（）函数获得解析解
C=real（fr）;                              %获得实部
S=imag（fr）;                              %获得虚部
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1884  2014-05-12 10:53  菲涅尔积分\Fresnel_demo.m

     文件      47424  2014-05-12 11:16  菲涅尔积分\菲涅尔积分.docx

     目录          0  2014-05-12 11:17  菲涅尔积分

----------- ---------  ---------- -----  ----

                49308                    3