MatLab在函数的求解方法
matlab函数求极限

matlab函数求极限使用Matlab函数求极限极限是数学中的一个重要概念,它描述了函数在某个点无限接近某个值的情况。
求解极限可以帮助我们理解函数的性质和行为,对于数学建模和问题求解也具有重要意义。
在Matlab中,我们可以利用一些函数来求解极限,本文将介绍一些常用的方法和技巧。
一、符号计算工具箱Matlab中的符号计算工具箱(Symbolic Math Toolbox)提供了一系列函数,可以用于进行符号计算和代数运算。
在求解极限时,可以使用其中的limit函数。
limit函数的调用格式为:limit(f, x, a)其中,f是要求极限的函数,x是自变量,a是自变量趋近的点。
通过调用limit函数,我们可以得到函数f在x趋近于a时的极限值。
例如,我们要求函数f(x) = sin(x)/x 在x趋近于0时的极限。
可以使用以下代码:syms xf = sin(x)/x;limit(f, x, 0)运行以上代码,可以得到极限值为1。
这是由于当x趋近于0时,sin(x)/x的值趋近于1。
二、数值计算工具箱除了符号计算工具箱,Matlab还提供了数值计算工具箱(Numerical Computing Toolbox),可以进行数值计算和数值优化。
在求解极限时,我们可以利用数值计算工具箱中的一些函数来进行近似计算。
例如,我们要求函数f(x) = (1+x)^(1/x) 在x趋近于0时的极限。
由于这个函数在x趋近于0时的极限不存在,我们可以通过数值计算来逼近这个极限值。
可以使用以下代码:x = linspace(-1, 1, 1000);f = (1+x).^(1./x);limit_value = f(end)运行以上代码,可以得到极限值为2.7183。
这是通过将x取一个足够小的范围,并计算函数在该范围内的取值来进行近似计算得到的。
三、图形方法除了数值计算,我们还可以利用图形方法来求解极限。
在Matlab中,可以使用plot函数绘制函数的图像,并观察函数在某个点附近的行为。
MATLAB计算方法3解线性方程组计算解法

MATLAB计算方法3解线性方程组计算解法线性方程组是数学中的一个重要问题,解线性方程组是计算数学中的一个基本计算,有着广泛的应用。
MATLAB是一种功能强大的数学软件,提供了多种解线性方程组的计算方法。
本文将介绍MATLAB中的三种解线性方程组的计算方法。
第一种方法是用MATLAB函数“linsolve”解线性方程组。
该函数使用高斯消元法和LU分解法求解线性方程组,可以处理单个方程组以及多个方程组的情况。
使用该函数的语法如下:X = linsolve(A, B)其中A是系数矩阵,B是常数向量,X是解向量。
该函数会根据A的形式自动选择求解方法,返回解向量X。
下面是一个使用“linsolve”函数解线性方程组的例子:A=[12;34];B=[5;6];X = linsolve(A, B);上述代码中,A是一个2×2的系数矩阵,B是一个2×1的常数向量,X是一个2×1的解向量。
运行代码后,X的值为[-4.0000;4.5000]。
第二种方法是用MATLAB函数“inv”求解逆矩阵来解线性方程组。
当系数矩阵A非奇异(可逆)时,可以使用逆矩阵求解线性方程组。
使用“inv”函数的语法如下:X = inv(A) * B其中A是系数矩阵,B是常数向量,X是解向量。
该方法先计算A的逆矩阵,然后将逆矩阵与B相乘得到解向量X。
下面是一个使用“inv”函数解线性方程组的例子:A=[12;34];B=[5;6];X = inv(A) * B;上述代码中,A是一个2×2的系数矩阵,B是一个2×1的常数向量,X是一个2×1的解向量。
运行代码后,X的值为[-4.0000;4.5000]。
第三种方法是用MATLAB函数“mldivide”(或“\”)求解线性方程组。
该函数使用最小二乘法求解非方阵的线性方程组。
使用“mldivide”函数的语法如下:X=A\B其中A是系数矩阵,B是常数向量,X是解向量。
matlab函数积分

matlab函数积分在MATLAB中,可以使用多种方法进行函数积分。
下面将详细介绍几种常用的方法。
1.基于符号计算的积分MATLAB的Symbolic Math Toolbox提供了一个功能强大的符号计算引擎,可以用于解析函数并求解积分。
首先,需要定义一个符号变量,然后使用int函数对其进行积分。
```matlabsyms x;f=x^2+3*x+2;integral_f = int(f, x);```这将返回一个符号表达式,表示函数f的积分。
如果要计算具体的数值积分,可以使用double函数对符号表达式进行求值。
```matlabnumerical_integral_f = double(integral_f);```这将返回函数f在积分区间上的数值积分结果。
2.数值积分对于无法通过符号方法求解的复杂函数,可以使用数值积分方法。
MATLAB提供了多种数值积分函数,其中最常用的是quad和quadl函数。
这些函数可以用于计算定积分和自适应积分。
```matlabintegral_f = quad(f, a, b);```这将返回函数f在积分区间[a, b]上的定积分结果。
quadl函数与quad函数类似,但可以处理更广泛的函数类型。
3.数值积分的误差控制在使用数值积分方法时,可以通过指定误差容限来控制积分的准确性。
例如,可以使用quad函数的相对误差容限选项来指定积分结果的相对误差范围。
```matlabintegral_opts = quadOptions('RelTol', 1e-6);integral_f = quad(f, a, b, integral_opts);```这将返回函数f在积分区间[a,b]上的定积分结果,并确保相对误差小于1e-64.多重积分MATLAB的Symbolic Math Toolbox还支持多重积分。
可以通过嵌套多个符号积分来进行多重积分的计算。
matlab求连续信号的频谱函数和离散信号频谱函数的方法

matlab求连续信号的频谱函数和离散信号频谱函数的方法Matlab提供了多种方法来求解连续信号和离散信号的频谱函数。
在本文中,我们将分步骤介绍这些方法。
一、连续信号频谱函数的方法连续信号的频谱函数是通过对连续信号进行傅里叶变换得到的。
而在Matlab中,傅里叶变换可以通过fft函数实现。
下面是求解连续信号频谱函数的步骤:1. 定义连续信号首先,我们需要定义一个连续信号,用一个函数来表示。
例如,我们定义一个简单的三角波信号:matlabt = linspace(0, 1, 1000); 定义时间范围x = sawtooth(2*pi*5*t); 定义三角波信号2. 进行傅里叶变换接下来,我们使用fft函数对连续信号进行傅里叶变换。
傅里叶变换将信号从时域转换到频域。
matlabX = fft(x);3. 计算频谱函数通过进行傅里叶变换,我们得到了频谱函数X。
然而,频谱函数X是一个复数数组,其中包含了信号的幅度和相位信息。
为了获得真正的频谱,我们需要计算幅度谱。
matlabP2 = abs(X/length(x));P1 = P2(1:length(x)/2+1);P1(2:end-1) = 2*P1(2:end-1);在上述代码中,我们将频谱函数除以信号长度,然后计算幅度,并使用对称性将频谱函数变换为正频率部分。
最后,我们将频谱函数的第一个和最后一个值乘以2。
4. 绘制频谱图最后,我们可以使用plot函数将频谱函数可视化。
matlabfs = 1000; 采样频率f = fs*(0:(length(x)/2))/length(x);plot(f,P1)xlabel('Frequency (Hz)')ylabel('Amplitude')以上步骤可以用于求解任何连续信号的频谱函数。
二、离散信号频谱函数的方法离散信号的频谱函数可以通过对信号进行离散傅里叶变换(Discrete Fourier Transform,DFT)来获得。
matlab中求函数的导数

matlab中求函数的导数MATLAB提供了几种不同的方法来计算函数的导数。
本文将介绍三种常用的方法:符号求导、数值求导和有限差分法。
1.符号求导符号求导是一种利用符号计算来找到函数导数的方法。
MATLAB的Symbolic Math Toolbox提供了符号计算的功能。
使用符号计算,可以求出任意复杂函数的导数。
以下是一个示例,展示了如何使用符号求导计算函数f(x)=x^2的导数:```matlabsyms xf=x^2;diff(f,x)```输出结果为:`2*x`符号求导的优点是可以得到一个精确的导数表达式,适用于数学函数和解析函数。
然而,计算符号导数可能需要大量的计算资源和时间,尤其是对于复杂的函数和高阶导数。
2.数值求导数值求导是一种使用数值方法计算函数导数的方法。
它基于函数在一些点的变化率来近似导数。
在MATLAB中,可以使用函数`diff`或`gradient`来进行数值求导。
以下是一个使用`diff`函数计算函数f(x) = x^2在x=1处的导数的示例:```matlabx=1;f=x^2;h=1e-6;%步长df = (f(x+h)-f(x))/h;```在数值求导中,步长h的选择对结果精度起着重要作用。
通常,较小的步长会导致较高的精度,但也会增加运算时间。
因此,需要在精度和效率之间找到一个平衡。
3.有限差分法有限差分法是一种数值计算方法,用于近似函数的导数。
它通过计算函数在邻近点上的差异来估计导数。
MATLAB中也有一些内置的函数用于计算导数,如`diff`, `gradient`和`diffusehess`等。
以下是一个使用`diff`函数计算函数f(x) = x^2在x=1处的导数的示例:```matlabx=1;f=x^2;h=1e-6;df = diff(f)/h;```有限差分法适用于函数没有解析表达式或难以求解的情况,它的运算速度相对符号求导和数值求导较快。
但是,有限差分法的精度受到步长h的约束,需要进行适当的调整以获得更精确的结果。
使用Matlab进行公式推导、求解

使用Matlab进行公式推导、求解MATLAB是一种强大的数值计算和科学编程软件,在MATLAB 中,可以使用其丰富的数学函数和符号计算工具进行公式推导和求解。
本文将以案例的形式介绍如何使用MATLAB进行公式推导和求解,包括符号计算、方程求解、微分和积分等方面的应用。
案例1:对以下公式进行“去括号展开”和“幂级数形式整理”(1)采用符号工具箱syms进行方程写入注意:代码建议写在实时编辑器中,这里语句都没有加分号,运行之后会直接显示结果如上图(2)“去括号展开”采用expand,“幂级数形式整理”采用series函数说明:"expand" 是 MATLAB 中用于展开代数表达式的一个函数。
它是 Symbolic Math Toolbox 的一部分。
当应用于符号表达式时,"expand" 函数会通过展开括号和简化项来展开表达式。
这在处理复杂的数学表达式或尝试简化和操作方程时非常有用。
"series" 是 MATLAB 中用于展开符号表达式的级数形式的一个函数。
它也属于 Symbolic Math Toolbox,提供了处理符号表达式和符号数学的功能,可以将符号表达式展开为指定的级数形式(‘order’后加展开到的最高指数),通常是泰勒级数。
它可以帮助在数学建模、分析和求解问题时处理复杂的表达式,并在需要时进行近似计算。
(3)案例完整代码clearsymsa0a1a2a3b0b1b2b3tt_ea=[a2;a1;a0];b=[b2;b1;b0] ;x(t)=[t^2t1]*a+a3*t*(t-t_e)^2y(t)=[t^2t1]*b+b3*t*(t-t_e)^2x(t)=expand(x)y(t)=expand(y)x(t)=series(x,t,'Ord er',4)y(t)=series(y,t,'Order',4)案例2:建立等式方程组,转换为线性矩阵表达式,然后求解(1)利用实时编辑器建立方程组我们定义了一些符号变量:a、b、x1、x2、y1、y2。
用Matlab求解函数的导数

用Matlab求解函数的导数标题:使用MATLAB求解函数的导数摘要:MATLAB是一种强大的数学软件,可用于解决各种数学问题。
本文将探讨如何使用MATLAB求解函数的导数。
我们将从简单的数值方法开始,逐步介绍MATLAB中提供的不同工具和技术,以获得更精确和高效的导数计算结果。
此外,我们还将分享对导数概念及其在数学和科学领域中的实际应用的理解。
导论:导数是微积分中一个重要的概念,用于描述函数在某一点的变化率。
求解函数的导数在数学、物理、工程等领域中具有广泛的应用。
使用MATLAB可以更方便地进行导数计算,并得到高质量的结果。
I. 数值导数方法在MATLAB中,最简单的求解导数的方法是使用数值差商近似。
通过计算函数在两个非常接近的点上的斜率来估计导数。
我们将演示如何使用数值差商近似来计算函数的导数,并讨论其精度和收敛性。
II. 符号导数计算MATLAB还提供了符号计算工具箱,可以通过符号表达式来求解函数的导数。
我们将介绍如何使用符号计算工具箱来获取函数的符号导数,并讨论符号计算与数值方法的比较。
III. 数值优化方法对于复杂的函数或需要高精度的导数计算,数值优化方法可以提供更准确的结果。
我们将介绍MATLAB中的几种高级数值优化方法,如梯度法和拟牛顿法,并演示如何在MATLAB中应用它们来求解函数的导数。
IV. 应用实例在本节中,我们将通过一些实际的应用示例来展示导数的重要性。
我们将通过MATLAB来解决一些典型的问题,如最小二乘拟合、优化问题和微分方程求解,以展示导数在不同领域中的实际应用。
总结与展望:通过本文,我们了解了如何使用MATLAB求解函数的导数。
我们从数值方法开始,逐步介绍了符号计算和数值优化方法,并演示了导数在实际问题中的应用。
MATLAB提供了丰富的工具和函数,能够满足不同需求的导数计算,并提供高质量的结果。
在今后的研究中,我们可以进一步探索MATLAB在数学建模、优化和控制等领域中的导数求解能力。
matlab如何求解三角函数方程

matlab如何求解三角函数方程(最新版)目录一、引言二、MATLAB 求解三角函数方程的方法1.使用符号计算函数2.使用 solve 函数3.使用 vpasolve 函数4.使用 fsolve 函数5.使用数值分析算法三、具体示例1.求解带有三角函数的方程组2.求解复杂三角方程四、结论正文一、引言在数学问题中,三角函数方程是一种常见的问题类型。
求解这类问题,我们可以使用 MATLAB 这一强大的数学软件。
MATLAB 提供了多种方法来求解三角函数方程,本文将对这些方法进行详细介绍。
二、MATLAB 求解三角函数方程的方法1.使用符号计算函数在 MATLAB 中,可以使用 symbolic 计算函数来求解三角函数方程。
例如,对于方程 sin(x) = y,可以使用以下命令求解:```matlabsyms x y;eq = sin(x) - y;solve(eq, x)```2.使用 solve 函数solve 函数是 MATLAB 中求解方程的常用函数。
对于三角函数方程,也可以使用 solve 函数求解。
例如,对于方程 sin(x) = y,可以使用以下命令求解:```matlabf = @(x) sin(x) - y;x0 = pi/4;[x, fval] = solve(f, x0)```3.使用 vpasolve 函数vpasolve 函数是 MATLAB 中用于求解带有参数的方程的函数。
对于三角函数方程,可以使用 vpasolve 函数求解。
例如,对于方程 a*sin(x) + b*cos(x) = c,可以使用以下命令求解:```matlabvpasolve("a*sin(x) + b*cos(x) = c", [a, b, c])```4.使用 fsolve 函数fsolve 函数是 MATLAB 中用于求解非线性方程的函数。
对于三角函数方程,可以使用 fsolve 函数求解。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
MatLab & 数学建模第四讲数值计算一、方程求解求解单个代数方程MATLAB具有求解符号表达式的工具,如果表达式不是一个方程式(不含等号),则在求解之前函数solve将表达式置成等于0。
>> solve( ' a*x^2+b*x+c ' ) % solve for the roots of the eqution ans=[1/2/a*(-b+(b^2-4*a*c)^1/2)][1/2/a*(-b-(b^2-4*a*c)^1/2)]结果是符号向量,其元素是方程的2个解。
如果想对非缺省x变量求解,solve 必须指定变量。
>> solve( ' a*x^2+b*x+c ' , ' b ' ) % solve for bans=-(a*x^2+c)/x带有等号的符号方程也可以求解。
>> f=solve( ' cos(x)=sin(x) ' ) % solve for xf=1/4*pi>> t=solve( ' tan(2*x)=sin(x) ' )t=[ 0][acos(1/2+1/2*3^(1/2))][acos(1/2=1/2*3^(1/2))]并得到数值解。
>> numeric(f)ans=0.7854>> numeric(t)ans=0 + 0.8314i1.9455注意在求解周期函数方程时,有无穷多的解。
在这种情况下,solve对解的搜索范围限制在接近于零的有限范围,并返回非唯一的解的子集。
如果不能求得符号解,就计算可变精度解。
>> x=solve( ' exp(x)=tan(x) ' )x=1.306326940423079代数方程组求解可以同时求解若干代数方程,语句solve(s1,s2,.....,sn)对缺省变量求解n个方程,语句solve(s1,s2,...,sn,' v1,v2,...,vn ')对n个' v1,v2,...vn '的未知数求解n个方程。
solve(f) 解符号方程式f。
solve(f1,…,fn) 解由f1,…,fn组成的联立方程式。
我们先定义以下的方程式:>>eq1 = 'x-3=4'; % 注意也可写成'eq1=x-7'>>eq2 = 'x*2-x-6=0'; % 注意也可写成'eq2=x*2-x-6'>>eq3 = 'x2+2*x+4=0';>>eq4 = '3*x+2*y-z=10';>>eq5 = '-x+3*y+2*z=5';>>eq6 = 'x-y-z=-1';>>solve(eq1)ans=7>>solve(eq2)ans=[[3],[-2]]' % 原方程式有二个根3, -2>>solve(eq3)ans=[[-1+i*3^(1/2)],[-1-i*3^(1/2)]]' % 注意实根和虚根的表示式>>solve(eq4,eq5,eq6) % 解三个联立方程式ans=x = -2, y = 5, z = -6如何处理中小学典型的代数问题?黛安娜(Diane)想去看电影,她从小猪存钱罐倒出硬币并清点,她发现:∙10美分的硬币数加上5美分的硬币总数的一半等于25美分的硬币数。
∙1美分的硬币数比5美分、10美分以及25美分的硬币总数多10。
∙25美分和10美分的硬币总数等于1美分的硬币数加上1/4的5美分的硬币数∙25美分的硬币数和1美分的硬币数比5美分的硬币数加上8倍的10美分的硬币数多1。
如果电影票价为3.00美元,爆米花为1.00美元,糖棒为50美分,她有足够的钱去买这三样东西?首先,根据以上给出的信息列出一组线性方程,假如p,n,d和q分别表示1美分,5美分,10美分,和25美分的硬币数dn pq p n d q q d pnq p n d++==++-+=++=+-210481然后,建立MATLAB符号方程并对变量求解。
>> eq1= ' d+(n+p)/2=q ' ;>> eq2= ' p=n+d+q-10 ' ;>> eq3= ' q+d=p+n/4 ' ;>> eq4= ' q+p=n+8*d-1 ' ;>>[pennies ,nickles ,dimes ,quarters]=solve(equ1,equ2,equ3,equ4,' p ,n ,d ,q ' ) pennies= 16 nickles= 8 dimes= 3quarters= 15所以,黛安娜有16枚1美分的硬币,8枚5美分的硬币,3枚10美分的硬币,15枚25美分的硬币,这就意味着>> money=.01*16+.05*8+.10*3+.25*15 money= 4.6100她就有足够的钱去买电影票,爆米花和糖棒并剩余11美分。
【例】求解二元函数方程组⎩⎨⎧=+==-=0)cos(),(0)sin(),(21y x y x f y x y x f 的零点。
(0)从三维坐标初步观察两函数图形相交情况x=-2:0.05:2;y=x;[X,Y]=meshgrid(x,y); %产生x-y 平面上网点坐标 F1=sin(X-Y);F2=cos(X+Y); F0=zeros(size(X)); surf(X,Y,F1),xlabel('x'),ylabel('y'), view([-31,62]),hold on, surf(X,Y,F2),surf(X,Y,F0), shading interp, %(间隔补齐) hold off图 5.6.3-0 两函数的三维相交图(1)在某区域观察两函数0等位线的交点情况clear;x=-2:0.5:2;y=x;[X,Y]=meshgrid(x,y); %产生x-y平面上网点坐标F1=sin(X-Y);F2=cos(X+Y);v=[-0.2, 0, 0.2]; %指定三个等位值,是为了更可靠地判断0等位线的存在。
contour(X,Y,F1,v) %画F1的三条等位线。
hold on,contour(X,Y,F2,v),hold off %画F2的三条等位线。
图 5.6.3-1 两个二元函数0等位线的交点图(2)从图形获取零点的初始近似值在图5.6.3-1中,用ginput获取两个函数0等位线(即三线组中间那条线)交点的坐标。
[x0,y0]=ginput(2); %在图上取两个点的坐标disp([x0,y0])-0.7926 -0.78430.7926 0.7843(3)利用fsolve求精确解。
以求(0.7926,7843)附近的解为例。
本例直接用字符串表达被解函数。
注意:在此,自变量必须写成x(1), x(2)。
假如写成xy(1), xy(2),指令运行将出错。
fun='[sin(x(1)-x(2)),cos(x(1)+x(2))]'; %<12> xy=fsolve(fun,[x0(2),y0(2)])%<13>xy =0.7854 0.7854(4)检验fxy1=sin(xy(1)-xy(2));fxy2=cos(xy(1)+xy(2));disp([fxy1,fxy2])1.0e-006 *-0.0994 0.2019〖说明〗●指令<12><13>可用以下任何一组指令取代。
(A)内联函数形式指令fun=inline('[sin(x(1)-x(2)), cos(x(1)+x(2))]', 'x'); %项'x'必须有。
xy=fsolve(fun,[x0(2), y0(2)]);(B)M函数文件形式及指令先用如下fun.m表示被解函数(并在搜索路径上)[fun.m]function ff=fun(x)ff(1)=sin(x(1)-x(2));ff(2)=cos(x(1)+x(2));然后运行指令xy=fsolve('fun',[x0(2),y0(2)]) 。
●第四步检验中的结果表明:所找零点处的函数值小于610 ,是一个十分接近零的小数。
该精度由options.TolFun控制。
options.TolFun的缺省值是1.0000e-006。
它可以用下列指令看到options=optimset('fsolve');options.TolFunans =1.0000e-006线性方程求解a= [ 7 2 1 -29 15 3 -2-2 -2 11 51 32 13]b=[4 7 -1 0]'x=a\bx =0.49790.14450.0629-0.0813单个微分方程常微分方程有时很难求解,MATLAB提供了功能强大的工具,可以帮助求解微分方程。
函数dsovle计算常微分方程的符号解。
因为我们要求解微分方程,就需要用一种方法将微分包含在表达式中。
所以,dsovle句法与大多数其它函数有一些不同,用字母D来表示求微分,D2,D3等等表示重复求微分,并以此来设定方程。
任何D后所跟的字母为因变量。
MATLAB解常微分方程式的语法是dsolve('equation','condition'),其中equation代表常微分方程式即y'=g(x,y),且须以Dy代表一阶微分项y'D2y代表二阶微分项y'',condition则为初始条件。
方程d y dx/=0用符号表达式D2y=0来表示。
独立变量可以指定或由symvar 22规则选定为缺省。
例如,一阶方程dy/dx=1+y2的通解为:>> dsolve( ' Dy=1+y^2 ' ) % find the general solutionans=-tan(-x+C1)其中,C1是积分常数。
求解初值y(0)=1的同一个方程就可产生:>> dsolve(' Dy=1+y^2 ',' y(0)=1 ') % add an initialconditiony=tan(x+1/4*pi)独立变量可用如下形式指定:>> dsolve(' Dy=1+y^2 ',' y(0)=1 ',' v ') % findsolution to dy/dvans=tan(v+1/4*pi)让我们举一个二阶微分方程的例子,该方程有两个初始条件:d ydx22=cos(2x)-y dy dx (0)=0 y(0)=1>> y=dsolve(' D2y=cos(2*x)-y ',' Dy(0)=0 ',' y(0)=1 ') y=-2/3*cos(x)^2+1/3+4/3*cos(x)>> y=simple(y) % y looks like it can be simplified y=-1/3*cos(2*x)+4/3*cos(x)通常,要求解的微分方程含有一阶以上的项,并以下述的形式表示:d y dx22-2dydx -3y=0 通解为:>> y=solve( 'D2y-2Dy-3*y=0 ')y=C1*exp(-x)+C2*exp(3*x)加上初始条件:y(0)=0和y(1)=1可得到:>> y=solve( ' D2y-2Dy-3*y=0 ' , ' y(0)=0,y(1)=1 ' ) y=1/(exp(-1)-exp(3))*exp(-x)-1/(exp(-1)-exp(3))*exp(3*x)>> y=simple(y) % this looks like a candidate for simplification y=-(exp(-x)-exp(3*x))/(exp(3)-exp(-1))>> pretty(y) % pretty it up exp(-x)-exp(3 x) - --------------------- exp(3) -exp(-1)现在来绘制感兴趣的区域内的结果。