如何用matlab求解二阶微分方程，以及程序实例

1、首先看一下matlab求解方程的方法，首先指明所解方程的变量，然后指明方程，未知数和限制条件，最后求解方程。

2、在matlab命令行窗口中输入symsx[x,params,conds]=solve(sin(x)==1,'ReturnConditions', true)，按回车键可以得到方程解。

3、转换一下，可以看到sin(x)=1方程的解是如下图所示 。

4、也可以求解下面的一个方程。

5、输入symsabcyx；[x,y]=solve([a*x^2+b*y+c==0,a*x+2*y==4],[x,y])。

6、按回车键可以得到方程解。

常用的窗函数：

bartlett，巴特利特窗口

调用格式：w = bartlett(L)，%L在列向量中返回一个点的Bartlett窗口w，其中L 必须是一个正整数。

blackman，布莱克曼窗口

调用格式：w = blackman(N)，%返回N列向量中的点对称Blackman窗口w，其中N是一个正整数。

bohmanwin，Bohman窗口

调用格式：w = bohmanwin(L)，%L在列向量中返回一个点Bohman窗口w。Bohman窗口是两个半持续时间余弦波瓣的卷积。在时域中，它是一个三角形窗口和一个余弦单周期的乘积，加上一个术语可以将边界处的一阶导数设置为零。Bohman窗脱落为1 / 瓦特^4。

chebwin，切比雪夫窗口

调用格式：w= chebwin(L,r)，%返回w包含长度LChebyshev窗口的列向量，  其傅里叶变换旁瓣幅度r 低于主瓣幅度dB。默认值为r100.0 dB。

gausswin，高斯窗口

调用格式：w = gausswin(N,Alpha)，%返回一个N与Alpha标准差的倒数成比例的高斯点高斯窗口。窗口的宽度与α的值成反比。α值越大，窗口越窄。的值α的默认值至2.5。

hamming，海明窗口

调用格式：w = hamming(L) ，%返回一个L点对称的Hamming窗口。

hann，汉恩（汉宁）窗口

调用格式：w = hann(L)，%返回一个L对称的Hann窗口。

kaiser，凯撒窗口

调用格式：w = kaiser(L,beta)，%返回L列向量中的-point Kaiser窗口w。beta是Kaiser窗口参数，它影响窗口傅里叶变换的旁瓣衰减。默认值为beta0.5。

rectwin，矩形窗口

调用格式：w = rectwin(L)，%返回L列向量中  的矩形窗口长度w。该功能是为了完整性而提供的 矩形窗口相当于没有窗口。

taylorwin，泰勒窗口

调用格式：w = taylorwin(n)，%n在列向量中返回一个点泰勒窗口，w。这个向量中的值是窗口权重或系数。

triang，三角窗口

调用格式：w = triang(L)，%返回L列向量中的一个点三角窗口。

具体应用场合可以查看有关数字信号处理方面的书籍。

我刚刚编了一个小程序，希望能帮助你。

clc

all=50

rent=1000

add=50

fix=100

In_max=0

for i=0:50

In=(rent+add*i)*(all-i)-fix*(all-i)

if In>In_max

In_max=Innumber=i

end

end

realrent=rent+add*number

In_max

结果输出：

realrent =

1800

In_max =

57800

---