日历
改变窗函数的长度，观察Gibbs现象
2009-5-30 20:06:11 阅读(1) 评论(0)
clear;
clc;
wc=0.6*pi;a=[1];Nf=256;
Hd1=ones(1,Nf-128);
Hd2=zeros(1,128);
Hd=[Hd1,Hd2]';
% 改变窗函数的长度，观察Gibbs现象
%
M=10; % for M=10
N=2*M+1;
for n=1:N
if (n-M)==0;
h(n)=wc/pi;
else
Gravity reduction formulae
2009-5-24 10:10:42 阅读(2) 评论(0)
Definitions
| go | = | observed gravity |
| gt | = | theoretical gravity at the surface of the reference ellipsoid |
| θ | = | latitude (degrees) |
| δg/ δz | = | + 0.3086 mGal m-1 (average vertical gravity gradient per metre of elevation above sea level) |
| G | = | 6.672 x 10-11m3kg-1s-2 or 6.672 |
ind2sub 和sun2ind
2009-5-22 13:12:25 阅读(0) 评论(0)
[i,j]=ind2sub([3 3],5)
% 3×3矩阵的第 5个元素的全下标
ind=sub2ind([3 3],3,3)
%矩阵第三行、第三列元素的序号
discrete fourier transform
2009-4-22 17:39:23 阅读(1) 评论(0)
C
C Subroutine to perform the discrete Fourier transform with
C FR and FI as the real and imaginary parts of the input and
C GR and GI as the corresponding output.
C
IFLAG=1 :FORWARD DISCRETE FOURIER TRANSFORM
IFLAG=-1:INVERSION DISCRETE FOURIER TRANSFORM
LP instruments
2009-4-12 13:17:46 阅读(1) 评论(0)
The science goals of LP were to map the Moon's surface composition and its magnetic and gravity fields, to determine the frequency and location of gas release events, and to search for polar ice deposits. To meet these objectives, LP had five science instruments located on three booms: a gamma ray spectrometer, a neutron spectrometer,
bessel transform (3)
2009-3-23 11:50:22 阅读(6) 评论(0)
C
NOFUN1=0
C-----ERROR CHECKS
IF(NB.LT.1.OR.NREL.LT.1.OR.BMAX.LE.0.0D0) GO TO 9999
Y=BMAX*ER**(NB-1)
IF(Y.LE.0.0D0) GO TO 9999
IERR=0
C-----INITIALIZE RELATED CONVOLUTION WITHIN LAGGED CONVOLUTION LOOPS
DO 10 I=1,801
bessel transform (2)
2009-3-23 11:49:10 阅读(10) 评论(0)
*WT0(761)/ 1.119169236611821D-10/,WT0(762)/-1.086836597292220D-10/,
*WT0(763)/ 1.055438043302322D-10/,WT0(764)/-1.024946589050432D-10/,
*WT0(765)/ 9.953360285547502D-11/,WT0(766)/-9.665809129205519D-11/,
*WT0(767)/ 9.386565284680683D-11/,WT0(768)/-9.115388754922578D-11/,
FAST REAL HANKEL TRANSFORM
2009-3-23 11:46:46 阅读(38) 评论(1)
欢迎继续探讨
C ALGORITHM 588, COLLECTED ALGORITHMS FROM ACM.
C ALGORITHM APPEARED IN ACM-TRANS. MATH. SOFTWARE, VOL.8, NO. 4,
C DEC., 1982, P.369-370, P. 344-468.
SUBROUTINE DHANKL(BMAX,NB,NREL,TOL,NTOL,NORD,FUN1,IJREL,DWORK, DHA 10
* DANS,ARG,NOFUN1,IERR)
第一类bessel 函数 fortran程序
2009-3-22 9:50:38 阅读(9) 评论(0)
C =====================================================
C Purpose: 计算第一类bessel函数及其导数 ( n= 0,1,。。。。)
C Input: x --- Argument of Jn(x) ( x 〉0 )
C n --- Order of Jn(x)
C Output: BJ(n+1) --- Jn(x)
C
matlab 对多组变量同时设置属性
2009-3-10 17:12:46 阅读(4) 评论(0)
t=0:0.05:5;
y=sin(t)';
k=1:6;
set(gcf,'DefaultAxesLineStyleOrder','-|:|-.');
set(gcf,'DefaultAxesColorOrder',[0 0 0;0.6 0.6 0.6]);
plot(t,y*k)
该种设置方法可以很好的克服matlab 中线型不足的问题。