ahp

前言 主要是通过matlab求解微分方程的实例，先来看一个小小问题： \[ x^2+y+(x-2y)y^{'}=0 \] 可以通过 matlab 的 dsolve 函数求得通解,两行代码就能够解决问题，关于 dsolve 的具体用法，参考官方文档的介绍 clc,clear syms x %定义符号变量 y=dsolve('x^2+y+(x-2*y)*Dy==0','x') %运行结果 % y = % x/2 + ((4*x^3)/3 + x^2 + C1)^(1/2)/2 % x/2 - ((4*x^3)/3 + x^2 + C1)^(1/2)/2 % C1是常数 问题一 求解常微分方程组的通解 代码： clc,clear syms f(x) g(x) [f1,g1]=dsolve('D2f+3*g==sin(x)','Dg+Df==cos(x)'); f1=simplify(f1) %化简 g1=simplify(g1) % 运行结果 % f1 = % C2 - sin(x)/3 + t*cos(x) + (3^(1/2)*C3*exp(3^(1/2)*t))/3 - (3^(1/2)*C4*exp(-3^(1/2)*t))/3 % g1 = % % sin(x)/3 - (3^(1/2)*C3*exp(3^(1/2)*t))/3 + (3^(1/2)*C4*exp(-3^(1

前言 主要是通过matlab求解微分方程的实例，先来看一个小小问题： \[ x^2+y+(x-2y)y^{'}=0 \] 可以通过 matlab 的 dsolve 函数求得通解,两行代码就能够解决问题，关于 dsolve 的具体用法，参考官方文档的介绍 clc,clear syms x %定义符号变量 y=dsolve('x^2+y+(x-2*y)*Dy==0','x') %运行结果 % y = % x/2 + ((4*x^3)/3 + x^2 + C1)^(1/2)/2 % x/2 - ((4*x^3)/3 + x^2 + C1)^(1/2)/2 % C1是常数 问题一 求解常微分方程组的通解 代码： clc,clear syms f(x) g(x) [f1,g1]=dsolve('D2f+3*g==sin(x)','Dg+Df==cos(x)'); f1=simplify(f1) %化简 g1=simplify(g1) % 运行结果 % f1 = % C2 - sin(x)/3 + t*cos(x) + (3^(1/2)*C3*exp(3^(1/2)*t))/3 - (3^(1/2)*C4*exp(-3^(1/2)*t))/3 % g1 = % % sin(x)/3 - (3^(1/2)*C3*exp(3^(1/2)*t))/3 + (3^(1/2)*C4*exp(-3^(1

代码： %% ------------------------------------------------------------------------ %% Output Info about this m-file fprintf('\n***********************************************************\n'); fprintf(' <DSP using MATLAB> Problem 8.35 \n\n'); banner(); %% ------------------------------------------------------------------------ % Digital Highpass Filter Specifications: wphp = 0.6*pi; % digital passband freq in rad wshp = 0.4*pi; % digital stopband freq in rad Rp = 0.5; % passband ripple in dB As = 60; % stopband attenuation in dB Ripple = 10 ^ (-Rp/20) % passband ripple in absolute Attn = 10 ^ (-As