chd_lkl的博客

【代码】matlab生成线性等间隔离散的梯形波。

clear; clcfm = 500; %主频dt = 1e-4; %时间域采样间隔number = 100; %采样点数t = -number/2+1:number/2;R = (1-2*(pi*fm*(t-number/4)*1e-3).^2) .* exp(-(pi*fm*(t-number/4)*1e-3).^2);plot(t,R,'linewidth', 1.5);title('Ricker-wavelet');xlabel('t(ms)');ylabel('A');s.

x = [-pi:0.01:pi]';y = sin(10*x);[~,loc1] = findpeaks(y); % 极大值[~,loc2] = findpeaks(max(y)-y); % 极小值

X = linspace(0,2*pi,25)';Y = sin(2*X);stem(X,Y,'LineStyle','-.',... 'LineWidth',1.25,... 'MarkerFaceColor','red',... 'MarkerEdgeColor','green')

计算卷积的物理意义与公式网上与书上很多，这里不再赘述。下面是计算两个向量u和v的fortran代码program test_conv !// 计算向量u和v的卷积 implicit none integer :: i, k integer, parameter :: m = 3, n = 8 !// 向量...

下面代码以2-30为数字区间，判断这个区间的哪些数字是素数logical :: log = .false.do n = 2, 30 do i = 2, n-1 if ( mod(n,i) == 0 ) then log = .true. print*...

使用ivf编译器Program main implicit none integer(kind=8) :: nn = 1000000000 !.. 设置抛硬币的总次数 integer(kind=8) :: m, n, i real(kind=8) :: a !.. a > 0.5: 设为硬币朝上；反之朝下 m = 0 !...