信息工程学院电子科学与技术系徐利民.ppt

课题3：股票预测 低频部分有两个峰值，局部放大后确定峰值位置为k=9和k=25 上图是全部频谱，下图是低频部分放大 课题3：股票预测 由k=9和k=25可以确定周期信号的频率： 用两个正弦信号近似股价曲线： 其中R1、I1、R2、I2分别是对应k=9和k=25时的傅立叶系数的实部和虚部。 课题3：股票预测 k=1:N; R1=real(Xc(9));I1=imag(Xc(9)); R2=real(Xc(25));I2=imag(Xc(25)); cc=(2/N)*(R1*cos(2*9*pi*k/2048)-I1*sin(2*9*pi*k/2048)+(R2*cos(2*25*pi*k/2048)-I2*sin(2*25*pi*k/2048))); plot(c,:k);hold;plot(cc,k); 拟和正弦曲线 2007/02/26-2008/04/18股价变化曲线(虚线)与股价预测曲线(实线) 真实股价曲线(2007/03/15-2008/07/22) 课题3：股票预测 课题4：乐音分析与合成 (1)利用MATLAB的声音处理的相关函数合成一段乐音。 1=C 4/4拍的乐曲《上海滩》的简谱： 课题4：乐音分析与合成 分析： 乐音的基本特征可以用基波频率，谐波频率和包络波形三个方面来表述。 乐音基波构成有一定规律。用C，D，E，F，G，A，B大写英文字母表示每个音的“音名”（或称为“音调”），音名对应固定的基波信号频率。 钢琴琴键分布图 每个音符频率值按“十二平均律”计算导出，计算规则略。 倍频程 (小字一组) 2倍关系 22/12倍 基波——音调；谐波——泛音列，音色 ω 0 Ω An 基波 基音 定音调 2Ω nΩ… 谐波 泛音(列) 定音色 设定乐音采样频率为 8kHz。 对这一段乐谱，根据前述的音乐知识，可确定C调中央频率是261.63Hz，依次推出乐谱每个音符的频率及持续时间并编程奏出乐曲 。 谐波决定音色，称为泛音列。 在合成的时候添加谐波可以使乐音音色听起来饱满优美。 用矩形、锯齿波代替单频正弦波做练习。 更真实的乐音合成还需要考虑乐音波形的包络以及音符的叠接。 添加包络： 实现叠接。 课题4：乐音分析与合成 思考与研究： 1）频谱图不能体现音符随时间变化的特点，怎么办？ 实际分析中常采用短时傅立叶变换(STFT)，绘制“语谱图”进行信号分析。 specgram(ys,[ ],Fs);%采用缺省数据绘制上述合成的乐曲《上海滩》语谱图 课题4：乐音分析与合成 合成乐曲《上海滩》的语谱图，清晰地刻画了每个音符的持续时间和频率结构 课题4：乐音分析与合成 2）什么叫“滑动窗”？什么叫“窄带”和“宽带”语谱图？ 3）利用miditoolbox进行midi音乐处理(miditoolbox工具箱及相关pdf文档已发给学员)。 4）研究音频编码 wav文件和mp3文件有何不同？什么是“知觉编码” ？ 课题4：乐音分析与合成 1 2 3 4 二、信号与系统的数值化处理 基础2：微分的数值计算 利用以下公式将微分化为差分： 后向Euler法 二、信号与系统的数值化处理 例：数值法求解微分方程， 初始条件为：y(0-)=y‘(0-)=1， 步长取TS=0.1，求解范围为区间[0,4]。 解：把微分方程转化为差分方程为 这些系数的转化如何不用人工计算？ 二、信号与系统的数值化处理 根据此方程求出差分方程的初始条件： 求解差分方程需要初始条件:y(-1),y(-2),怎么办？ 上例解得y(-1)=1；y(-2)=0.9。编程如下： 二、信号与系统的数值化处理 Ts=0.1; %时间间隔 n=1:4/Ts+2; %坐标1,2代表初始条件坐标-2,-1 F(n)=1;F(1)=0;F(2)=0; %激励信号初始条件 Y(1)=0.9;Y(2)=1; %输入初始条件y(-2)=0.9;y(-1)=1 for k=3:length(n) %迭代求解 Y(k)=(13*F(k)-10*F(k-1)-100*Y(k-2)+240*Y(k-1))/143; end n1=(n-3)*Ts; subplot(2,1,1);plot(n1,Y); xlabel(t);ylabel(y(t));title(全响应(数值计算));grid on; %下图用实际响应进行对比 subplot(2,1,2);ezplot(-exp(-3*t)+exp(-t)+1,[0,4]); xlabel(t);ylabel(y(t));title(全响应);grid on; 二、信号与