基于MATLAB的金融工程方法及应用 资料_ 附录.pptVIP

15.5 插值 15.5.1 一维插值（interp1） 【例】 p=1:1:10;q=3*p q = 3 6 9 12 15 18 21 24 27 30 p1=3.5 p1 = 3.5000 q1=interp1(p,q,p1,spline) q1 = 10.5000 15.65插值 15.5.2 样条插值（spline） 样条插值对应的函数是 spline，使用方法如下： qq = spline (p,q,pp) 其中，自变量p是单调增加的。 【例】 q1=spline(p,q,p1) q1 = 10.5000 15.5 插值 15.5.3 Hermite插值 在实际应用中我们可能还需要节点处的导数是相等的，而Hermite插值方法可以保持节点处的导数相等，在编写程序中为实现这种插值方法所需调用的函数是 pchip。 【例】 q1=pchip(p,q,p1) q1 = 10.5000 上文中自变量p1也可为向量变量，如p1=[3.5,4.5,5.5]。 15.6 符号计算 15.6.1 符号对象 建立符号变量和符号常数 (1)sym函数 sym函数用来建立单个符号量，例如，x=sym(x)建立符号变量x，此后，用户可以在表达式中使用变量x进行各种运算。 【例】 a=sym(a);b=sym(b);c=sym(c);d=sym(d); %定义4个符号变量 A=[a,b;c,d]; %建立符号矩阵A det(A); %计算符号矩阵A的行列式 ans = a*d-b*c 15.6 符号计算 15.6.1 符号对象 建立符号变量和符号常数 (2)syms函数 运用syms函数可以批量定义符号变量var1,var2,…,varn等。用这种格式定义符号变量时不要在变量名上加字符分界符()，变量间用空格而不要用逗号分隔,调用格式如下： syms var1 var2 … varn 15.6 符号计算 15.6.1 符号对象 建立符号表达式 在MATLAB中建立符号表达式主要有以下3种方法： (1)利用单引号来生成符号表达式。 (2)用sym函数建立符号表达式。 (3)使用已经定义的符号变量组成符号表达式。 15.6 符号计算 15.6.2 符号表达式的四则运算 符号表达式的四则运算和其他表达式的运算并无不同，但要注意运算结果依然是一个符号表达式。符号表达式的加、减、乘、除运算可分别由函数symadd、symsub、symmul和symdiv来实现，幂运算可以由sympow来实现。 15.6 符号计算 15.6.3 因式分解与展开 factor(S) 对S分解因式，S是符号表达式或符号矩阵。 expand(S) 对S进行展开，S是符号表达式或符号矩阵。 collect(S) 对S合并同类项，S是符号表达式或符号矩阵。 collect(S,v) 对S按变量v合并同类项，S是符号表达式或符号矩阵。 【例】 factor(p^2+2*p*q+q^2) ans = (q+p)^2 collect(x^5+2*x*y+x*y^3) ans = x^5+(2*y+y^3)*x 15.6 符号计算 15.6.4 符号矩阵 transpose(S) 返回S矩阵的转置矩阵。 determ(S) 返回S矩阵的行列式值。 colspace(S) 返回S矩阵列空间的基。 ? 15.6 符号计算 15.6.5 符号函数求导及其应用 MATLAB中可以实现求导的函数为： diff(f,x,n) diff函数求函数f对变量x的n阶导数。参数x的用法同求极限函数limit，可以缺省，缺省值与limit相同，n的缺省值是1。 ? 15.6 符号计算 15.6.5 符号函数求导及其应用 【例】在曲线y=x3+3x-2上哪一点的切线与直线y=5x+1平行。 命令如下： x=sym(x); y=y=2*x^3+4*x^2-5; %定义曲线函数 d=diff(y); %对曲线求导数 r=f-5; solve? %求方程f-4=0的根，即求曲线何处的导数为4 ans = -2/3+1/6