草庐IT

MATLAB——观察吉布斯效应

吉布斯效应:当原信号是脉冲信号时,其高频分量主要影响脉冲的跳变沿,低频分量主要影响脉冲的顶部,因此,输出信号波形总是要发生失真,该现象称为吉布斯现象。原理:任意周期信号表示为傅立叶级数时,需要无限多项才能逼近原信号,但在实际应用中经常采用有限项级数来代替无限项级数。所选项数越多越接近原信号。所以用一个以原点为中心奇对称的周期性方波,可以用奇次正弦波的叠加来逼近可以观察出吉布斯效应。代码解释:主要是下列代码不好理解N=round((nf+1)/2);y=zeros(N,max(size(t)));x=zeros(size(t));fork=1:2:nfx=x+sin(w1*k*t)/k;y((k

什么是吉布斯采样(Gibbs Sampling)

目录1蒙特卡洛方法1.1蒙特卡洛方法的作用1.2非均匀分布采样1.3分布p(x)不好采样怎么办?2什么是吉布斯采样2.1马尔可夫链2.1.1什么是马尔可夫链呢?2.1.2为什么我们要引入马尔可夫链?2.1.3对给定的分布π\piπ,怎么找到对应的P,使得其为平稳马尔可夫过程2.2MCMC采样2.3M-H采样2.4吉布斯采样(Gibbs)2.4.1吉布斯采样原理2.4.1.1二维情况2.4.1.2高维情况2.4.2吉布斯采样过程参考资料1蒙特卡洛方法介绍吉布斯采样前,我们先看一下蒙特卡洛方法。1.1蒙特卡洛方法的作用有很多函数我们无法直接得到他的积分值,但我们可以利用蒙特卡洛方法来进行估计。比如