MATLAB编程基础符号积分变换符号方程求解可视化重点.pptVIP

符号积分变换、符号方程求解、可视化分析 ;4.5 积分变换 常见的积分变换有傅立叶变换、拉普拉斯变换。 4.5.1．傅立叶(Fourier)变换 在MATLAB中，进行傅立叶变换的函数是： fourier(f,t,w)：求函数f(t)的傅立叶像函数F(w)。 ifourier(F,w,t)：求傅立叶像函数F(w)的原函数f(t)。;例4-22 Fourier变换和反变换的简单示例 syms t a b c; f = exp(-t^2-a^2); % 创建符号函数 F = fourier(f,a,b) % 求时域函数 的傅立叶变换 g = ifourier(F,c) % 求频域函数 的傅立叶反变换，频域自变量由函数findsym自动确认为;傅立叶分析;fs=100;%设定采样频率N=128;n=0:N-1;t=n/fs;f0=10;%设定正弦信号频率%生成正弦信号x=sin(2*pi*f0*t);figure(1);subplot(131);plot(t,x);%作正弦信号的时域波形xlabel(t);ylabel(y);title(正弦信号y=2*pi*10t时域波形);grid;%进行FFT变换并做频谱图y=fft(x,N);%进行fft变换mag=abs(y);%求幅值m=length(y);f=(0:m/2-1)*fs/m;%进行对应的频率转换figure(1);subplot(132);plot(f,mag(1:m/2));%做频谱图axis([0,100,0,80]);xlabel(频率(Hz));ylabel(幅值);title(正弦信号y=2*pi*10t幅频谱图N=128);grid;%求均方根谱sq=abs(y);figure(1);subplot(133);plot(f,sq(1:m/2));xlabel(频率(Hz));ylabel(均方根谱);title(正弦信号y=2*pi*10t均方根谱);grid;;;4.5.2 Laplace变换及反变换 在MATLAB中，进行拉普拉斯变换的函数是： laplace(fx,x,t)：求函数f(t)的拉普拉斯像函数L(s)。 ilaplace(Fw,t,x)：求拉普拉斯像函数L(s)的原函数f(t)。;例4-23 Laplace变换和反变换例子 syms t a y; f = exp(a*t); % 创建符号函数 L = laplace(f) % 求函数f的拉普拉斯变换，参数缺省 g = ilaplace(L,y) % 求反拉普拉斯变换，f函数自变量为y 训练任务： 计算y=x3的拉普拉斯变换及其逆变换。;4.6 符号方程求解 4.6.1 符号代数方程求解 一、线性方程 利用sym转为符号后，左除即可。 二、一般方程 在MATLAB中，求解用符号表达式表示的代数方程可由函数solve实现，其调用格式为： solve(s)：求解符号表达式s的代数方程，求解变量为默认变量。 solve(s,v)：求解符号表达式s的代数方程，求解变量为v。 solve(s1,s2,…,sn,v1,v2,…,vn)：求解符号表达式s1,s2,…,sn组成的代数方程组，求解变量分别v1,v2,…,vn。 ;;例4-26 函数solve()求解方程 syms a b c x y; f = a*x^2+b*x+c; % 创建符号方程，方程以符号表达式的形式给出 sf1 = solve(f) % 求解方程f，指定变量缺省 sf2 = solve(f, b) % 求方程f对变量b的解 f2 = x+y=2*a; f3 = 3*x+y-1=7*a; % 创建符号方程组 sf3 = solve(f2,f3) % 求解方程组，指定变量缺省 % 变量自动确认为x和y，显示x和y的值 sf3.x sf3.y sf4 = solve(f2,f3,a,x,y) % 求方程组对指定变量的解 ;;4.6.2 符号常微分方程求解 在MATLAB中，用大写字母D表示导数。例如，Dy表示y，D2y表示y，Dy(0)=5表示y(0)=5。D3y+D2y+Dy-x+5=0表示微分方程y+y+y-x+5=0。符号常微分方程求解可以通过函数dsolve来实现，其调用格式为： dsolve(e,c,v) 该函数求解常微分方程e在初值条件c下的特解。参数v描述方程中的自变量，省略时按缺省原则处理，若没有给出初值条件c，则求方程的通解。 dsolve在求常微分方程组时的调用格式为：dsolve(e1,e2,…,en,c1,…,cn,v1,…,vn) 该函数求解常微分方程组e1,…,en在初值条件c