金融MATLAB实验报告三答案详解
MATLAB原理及应用实验报告第三章答案
![MATLAB原理及应用实验报告第三章答案](https://img.taocdn.com/s3/m/1099673f6c85ec3a87c2c5de.png)
《MATLAB 原理及应用》实验报告9.课后练习题答案1.提取符号表达式222m z u ar +++的自由变量(当符号表达式中含有多余一个符号变量时,只有一个变量是独立变量,其余的符号当作常量。
如果不指定那一个变量当作是自由变量,matlab 将基于一定原则选择一个自由变量。
) 提示:findsym (s ,n )程序如下:>> syms a r u z mh=sym('a*r^2+u^2+z+m^2');>> findsym (h)ans =a, m, r, u, z 2.在MATLAB 中计算多项式2352212-+++-x x x x 的父母和分子 (提示:使用[n ,d]=numden (A ))程序如下:>> s=sym('(x^2-1)/(x+2)+(2*x+5)/(3*x-2)');[n,d]=numden(s)n =3*x^3+6*x+12d =(x+2)*(3*x-2)2.1、建立符号函数5332232-+x x (1)提取该表达式的分子和分母,并分别付给两个变量(2)对这两个变量分别进行代数运算(加减乘除及乘方)程序如下:>> s=sym('3/2*x^2+2/3*x-3/5');>> [n,d]=numden(s)n =45*x^2+20*x-18d =30>> s1=n-d;>> s2=n+d;>> s3=n*d;>> s4=n/d;>> s5=n^d;>> s=[s1;s2;s3;s4;s5]s =45*x^2+20*x-4845*x^2+20*x+121350*x^2+600*x-5403/2*x^2+2/3*x-3/5(45*x^2+20*x-18)^303.在MA TLAB 中,按照不同的方式合并表达式)()(23x e y x y x e x y y --+*-+的参数类项。
MATLAB实验三参考答案
![MATLAB实验三参考答案](https://img.taocdn.com/s3/m/a1d30c07866fb84ae45c8db2.png)
how =collect(x)
4、求下列函数的极限(写出命令) (1) lim
cos x e x 0 x4
x2 2
syms x; limit('(cos(x)-exp(-1/2*x^2))/(x^4)',x,0) -1/12 (2) lim
2 x ln 2 x 1 x 0 1 cos x
syms n; S=symsum(1/((3*n-2)*(3*n+1)),n,1,inf) 8、试求出函数 f ( x )
sin x 的麦克劳林幂级数展开式的前 9 项,并求出关于 x=2 x 3x 2
2
的 Taylor 幂级数展开式的前 5 项。(命令 taylor 或者 taylortool)
河南财经政法大学数学与信息科学学院 1
实验报告
结果: EXPR =(x^2+x*exp(-t)+1)*(x+exp(-t)) expr1 =x^3+2*exp(-t)*x^2+(1+exp(-t)^2)*x+exp(-t) expr2 =x*exp(-t)^2+(2*x^2+1)*exp(-t)+(x^2+1)*x 3、factor(因式分解),simple(简化运算,对表达式尝试多种不同的算法进行简化,并以最 简化形式给出,How 中记录的为简化过程中使用的方法, )指令的使用 syms a x; f1=x^4-5*x^3+x^2+5*x-6; factor(f1) x^4-5*x^3+x^2+5*x-6 f2=x^2-a^2; factor(f2) (x-a)*(x+a) f3=2*sin(x)^2-cos(x)^2 [y , how]=simple(f3) y 为 f 的最优化简形式,How 中记录的为简化过程中使用的方法 y =-3*cos(x)^2+2 how =simplify [y , how]=simple(f1) y =x^4-5*x^3+x^2+5*x-6
matlab2011实验3参考答案
![matlab2011实验3参考答案](https://img.taocdn.com/s3/m/ae7ec19d6bec0975f465e240.png)
MATLAB 实验三MATLAB 计算可视化试验报告说明:1 做试验前请先预习,并独立完成试验和试验报告。
2 报告解答方式:将MATLAB 执行命令和最后运行结果从命令窗口拷贝到每题的题目下面,请将报告解答部分的底纹设置为灰色,以便于批阅。
3 在页眉上写清报告名称,学生姓名,学号,专业以及班级。
3 报告以Word 文档书写。
文档命名方式: 学号+姓名+_(下划线)+试验几.doc 如:110400220张三_试验1.doc 4 试验报告doc 文档以附件形式发送到maya_email@ 。
凡文档命名不符合规范,或者发送方式不正确,不予登记。
一 目的和要求(1)熟练掌握MATLAB 二维曲线的绘制(2)熟练掌握图形的修饰(3)掌握三维图形的绘制(4)熟练掌握各种特殊图形的绘制(5)掌握句柄图形的概念二 内容和步骤(说明:绘图题需在报告中粘贴图形结果。
)1 绘制二维曲线(必做)查阅stem 绘图指令,利用stem 绘制复指数序列n j e n x )3..01.0()(+-=(其中1010<<-n )的相关图形,要求:(1)把窗口分成2行2列。
(2)左上窗口绘制复指数序列的实部波形,窗口右上角绘制复指数序列的虚部波形,左下角绘制复指数序列的幅度波形,右下角绘制窗口的相位波形。
要求每个子图标注标题,X,Y 轴说明。
要求显示网格。
% 复指数序列程序n=-10:10;alpha=-0.1+0.3*j;x=exp(alpha*n);real_x=real(x); % 实部img_x=imag(x); % 虚部mag_x=abs(x); % 幅值phase_x=(180/pi)*angle(x); % 相角,转为角度subplot(2,2,1)stem(n,real_x);title('实数部分');xlabel('n');ylabel('real'); grid onsubplot(2,2,2)stem(n,img_x);title('虚数部分');xlabel('n');ylabel('imag');grid onsubplot(2,2,3)stem(n,mag_x);title('幅度');xlabel('n');ylabel('abs');grid onsubplot(2,2,4)stem(n,phase_x);title('相位');xlabel('n');ylabel('phase');grid on同学解题过程存在的问题:1注意括号应用,有同学如此定义指数序列:x=exp(-0.1+j*0.3)*n ,思考错在哪里?2绘图程序最好放在M 脚本中编辑,部分同学是在命令窗口中,这样不便于对程序的调试和修正。
matlab实验三答案详解
![matlab实验三答案详解](https://img.taocdn.com/s3/m/675b8e11ff00bed5b9f31d8e.png)
实验三选择结构程序设计实验内容1.程序代码function f(x)if x<0&x~=-3y=x*x+x-6;elseif x>0&x<5&x~=2&x~=3y=x*x-5.*x+6;elsey=x*x-x-1endy运行结果>> f(-7)y =36>> f(2.5)y =-0.2500>> f(90)y =80092.程序代码score=input('分数=');if score<=100&score>=0switch fix(score/10)case{9,10}grade='A';case{8}grade='B';case{7}grade='C';case{6}grade='D';otherwisegrade='E';endgradeelsedisp('输入有误,请输入正确的百分制成绩');end运行结果>> 分数=56grade =E>> f分数=345输入有误,请输入正确的百分制成绩3.程序代码time=input('工作时间=');switch timecase{time>120}wage=(time-120)*(1+15/100)*84+120*84;case{time<60}wage=time*84-700;otherwisewage=time*84end运行结果>> f工作时间=150wage =12600>> f工作时间=50wage =4200>> f工作时间=110wage =92404.程序代码A=10+90*rand(5);a=fix(A(3,3))b=fix(A(2,4))x=input('输入一个运算符号','s');switch xcase {'+'}c=a+b;case{'-'}c=a-b;case{'*'}c=a*b;case{'/'}c=a/b;otherwisec='无效'endc运行结果>> fa =55b =59输入一个运算符号+c =114>> fa =77b =14输入一个运算符号/c =5.5000>> fa =30b =58输入一个运算符号#c =无效c =无效5.程序代码A=input('输入一个5行6列的矩阵A=');n=input('输入一个正整数n=');if n<5&n>0disp(A(n,:));elsedisp(A(5,:));endlasterr运行结果>> f输入一个5行6列的矩阵A=[1 2 3 4 5 5;2 3 4 5 7 6;2 9 2 2 2 3;11 2 11 9 7 3;2 3 4 5 6 7] 输入一个正整数n=62 3 4 5 6 7ans =Undefined function or variable 'clcl'.>> f输入一个5行6列的矩阵A=[1 2 3 4 5 5;2 3 4 5 7 6;2 2 2 2 2 3;11 2 3 9 7 3;2 3 4 5 6 7]输入一个正整数n=32 2 2 2 2 3ans =Undefined function or variable 'clcl'.。
(完整版)MATLAB)课后实验答案[1]
![(完整版)MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/306ff29b690203d8ce2f0066f5335a8102d266cb.png)
(完整版)MATLAB)课后实验答案[1]实验⼀ MATLAB 运算基础1. 先求下列表达式的值,然后显⽰MATLAB ⼯作空间的使⽤情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =,其中2120.455i x +??=?- (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--L (4) 2242011122123t t z t t t t t ?≤=-≤,其中t =0:0.5:2.5 解:4. 完成下列操作:(1) 求[100,999]之间能被21整除的数的个数。
(2) 建⽴⼀个字符串向量,删除其中的⼤写字母。
解:(1) 结果:(2). 建⽴⼀个字符串向量例如:ch='ABC123d4e56Fg9';则要求结果是:实验⼆ MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S=?,其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对⾓阵,试通过数值计算验证2 2E R RS A O S +??=。
解: M ⽂件如下;5. 下⾯是⼀个线性⽅程组:1231112340.951110.673450.52111456x x x ??=???(1) 求⽅程的解。
(2) 将⽅程右边向量元素b 3改为0.53再求解,并⽐较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解: M ⽂件如下:123d4e56g9实验三选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ?+-<≠-?=-+≤<≠≠??--?且且及其他⽤if 语句实现,分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
解:M ⽂件如下:2. 输⼊⼀个百分制成绩,要求输出成绩等级A、B、C、D、E。
matlab-第三次实验报告
![matlab-第三次实验报告](https://img.taocdn.com/s3/m/8712168d8762caaedd33d4a3.png)
Matlab 第三次实验报告 2012029010010 尹康1. 对高阶多项式()()()()()2011220k p x x x x x k ==---=-∏编程求下面方程的解()190p x x ε+=并绘图演示方程的解与扰动量ε的关系。
解:程序代码:p=[1,-1];for i=2:20n=[1,-i];p=conv(p,n); % 求多项式乘积 endm=zeros(1,21); % m 的最高次幂为20,有21项 hold onx=1:20;d=[-1,0,0.1,0.5,1];for i=1:5delt=d(i);m(2)=delt;y=(roots(p+m))'; % 求多项式的根 plot(x,y,'-o','color',[i/5,i/20,i/10]); endtitle('方程p(x)=0的解与扰动量delt 的关系')legend('delt=-1','delt=0','delt=0.1','delt=0.5','delt=1') 运行结果:心得体会:1.当delt=0时,程序运行结果仍然不能精确得到理论解,进一步说明计算机的表达能力是有限的,计算机中的运算总是存在舍入误差; 2当delt ≠0时,总的来说,即使|delt|很小,它对解的扰动也是巨大的;考察它对每个解的扰动程度,发现解的理论值越大,那么它受到的扰动越大; 3.当|delt|≠0的绝对值一定时,它是负值所造成的扰动大于它是正值时的。
2. 对220n = ,生成对应的Hilbert 矩阵,计算矩阵的条件数;通过先确定解获得常向量b 的方法,确定方程组()n H x b =最后,用矩阵分解方法求解方程组,并分析计算结果。
解:程序代码:clear allclcfor n=2:20h=hilb(n);fprintf('\n\nn=%-10dcond(Hn)∞=%d',n,cond(h,'inf'))%cond :求矩阵范数X=1:n;b=h*X';[l u]=lu(h); %lu分解x=u\(l\b); %利用lu分解求线性其次方程组的根 x=x';fprintf('\n X’=')for i=1:nfprintf('%-8.2f',X(i))endfprintf('\n x’=')for i=1:nfprintf('%-8.2f',x(i))endend运行结果:心得体会:1.采用lu 分解求解(x=u\(l\b ))的效率大大高于直接求解(x=A\b ); 2.当n ≤10时,尽管系数矩阵Hn 的条件数很大,但该算法所得解仍然能与理论解很好地吻合,说明病态问题并非一定不能求解;当n ≥11时(截图中红线所指),系统提示:Warning: Matrix is close to singular or badly scaled. Results may be inaccurate ,所得解的误差已经不能忽略。
MATLAB实验及答案详解
![MATLAB实验及答案详解](https://img.taocdn.com/s3/m/610d39e7910ef12d2af9e7ac.png)
>> C3
C3 =
9 22 2
8 19 -6
36 88 3
>> D1
D1 =
2 -2 3
655
16 -10 3
>> D2
D2 =
0.5000 -2.0000 0.3333
0.6667 5.0000 0.2000
1.0000 -10.0000 3.0000
>> D3
D3 =
141
4 25 1
16 100 9
5.0000 - 5.0000i 6.0000 - 6.0000i
②求矩阵的逆矩阵
>> B=[1 2;3 4];
>> B1=inv(B) 运行后显示:
B1 = -2.0000 1.0000 1.5000 -0.5000
>> B2=B^(-1) 运行后显示: B2 =
-2.0000 1.0000 1.5000 -0.5000 ③关于矩阵求幂 >>A=[1 2;3 4]; B=[ 2 1;3 2]; >>A.^B 运行后显示: ans =
>>A=[1,2,3;4,5,6] 运行后显示: A= 123 456
在命令窗口输入:
>>b=A(1,2)
运行后显示:b =2
在命令窗口输入:
>>A(2,3)=-3
运行后显示:A =
123
4 5 -3
矩阵的操作
>>A=[1,2,3;4,5,6;7,8,9]
>>B=diag(A) %X 为矩阵时,V=diag(X,k)得到列向量 V,它取自 X 的第 K 个对角
matlab2009实验3参考解答
![matlab2009实验3参考解答](https://img.taocdn.com/s3/m/64582e80d4d8d15abe234e5e.png)
MATLAB 实验三MATLAB 计算可视化试验报告说明:1 做试验前请先预习,并独立完成试验和试验报告。
2 报告解答方式:将MATLAB 执行命令和最后运行结果从命令窗口拷贝到每题的题目下面,请将报告解答部分的底纹设置为灰色,以便于批阅。
3 在页眉上写清报告名称,学生姓名,学号,专业以及班级。
3 报告以Word 文档书写。
文档命名方式: 学号+姓名+_(下划线)+试验几.doc 如:110400220张三_试验1.doc 4 试验报告doc 文档以附件形式发送到maya_email@ 。
凡文档命名不符合规范,或者发送方式不正确,不予登记。
5 每次试验报告的最后提交期限:下周试验课之前。
6 欢迎访问我的教学论坛:/nbbs ,由此进入Matlab 专栏。
论坛注册时请在个人签名档写上自己的姓名和学号,以便于权限设置,可在论坛上交流学习,并可下载发布的试验报告和解答参考等相关资料。
一 目的和要求(1)熟练掌握MATLAB 二维曲线的绘制(2)熟练掌握图形的修饰(3)掌握三维图形的绘制(4)熟练掌握各种特殊图形的绘制(5)掌握句柄图形的概念(6)了解GUI 设计二 内容和步骤(说明:绘图题需在报告中粘贴图形结果。
)1 绘制二维曲线(必做)查阅stem 绘图指令,利用stem 绘制复指数序列n j e n x )3..01.0()(+-=(其中1010<<-n )的相关图形,要求:(1)把窗口分成2行2列。
(2)左上窗口绘制复指数序列的实部波形,窗口右上角绘制复指数序列的虚部波形,左下角绘制复指数序列的幅度波形,右下角绘制窗口的相位波形。
要求每个子图标注标题,X,Y 轴说明。
要求显示网格。
% 复指数序列程序n=-10:10;alpha=-0.1+0.3*j;x=exp(alpha*n);real_x=real(x); % 实部img_x=imag(x); % 虚部mag_x=abs(x); % 幅值phase_x=(180/pi)*angle(x); % 相角,转为角度subplot(2,2,1)stem(n,real_x);title('实数部分');xlabel('n');ylabel('real');grid onsubplot(2,2,2)stem(n,img_x);title('虚数部分');xlabel('n');ylabel('imag');grid onsubplot(2,2,3)stem(n,mag_x);title('幅度');xlabel('n');ylabel('abs');grid onsubplot(2,2,4)stem(n,phase_x);title('相位');xlabel('n');ylabel('phase');grid on同学解题过程存在的问题:1注意括号应用,有同学如此定义指数序列:x=exp(-0.1+j*0.3)*n,思考错在哪里?2绘图程序最好放在M脚本中编辑,部分同学是在命令窗口中,这样不便于对程序的调试和修正。
MATLAB全部实验及答案
![MATLAB全部实验及答案](https://img.taocdn.com/s3/m/0a4032f10242a8956bece473.png)
MATLAB全部实验及答案实验一、MATLAB基本操作实验内容及步骤1、命令窗口的简单使用(1)简单矩阵的输入(2)求[12+2×(7-4)]÷32的算术运算结果2、有关向量、矩阵或数组的一些运算(1)设A=15;B=20;求C=A+B与c=a+b?(2)设A=[1 2 3;4 5 6;7 8 9],B=[9 8 7;6 5 4;3 2 1];求A*B与A.*B?A*B就是线代里面的矩阵相乘 A.*B是对应位置的元素相乘(3)设a=10,b=20;求i=a/b=0.5与j=a\b=2?(4)设a=[1 -2 3;4 5 -4;5 -6 7]请设计出程序,分别找出小于0的矩阵元素及其位置(单下标、全下标的形式),并将其单下标转换成全下标。
clear,clca=[1 -2 3;4 5 -4;5 -6 7];[x,y]=find(a<0);c=[];for i=1:length(x)c(i,1)=a(x(i),y(i));c(i,2)=x(i);c(i,3)=y(i);c(i,4)=(y(i)-1)*size(a,2)+x(i);endc(5)在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8];看结果如何?如果改成运行A=[1,2;3,4]+i[5,6;7,8],结果又如何?前面那个是虚数矩阵,后面那个出错(6)请写出完成下列计算的指令:a=[1 2 3;3 4 2;5 2 3],求a^2=?,a.^2=?a^2= 22 16 1625 26 2326 24 28a.^2=1 4 99 16 425 4 9(7)有一段指令如下,请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X( : ) 转化为列向量(8)使用三元组方法,创建下列稀疏矩阵2 0 8 00 0 0 10 4 0 06 0 0 0方法一:clear,clcdata=[2 8 1 4 6];ir=[1 1 2 3 4 ];jc=[1 3 4 2 1];s=sparse(ir,jc,data,4,4);full(s)方法二:不用三元组法clear,clca=zeros(4,4);a(1,[1,3])=[2,8];a(2,4)=1;a(3,2)=4;a(4,1)=6;a(9) 写出下列指令的运行结果>> A = [ 1 2 3 ]; B = [ 4 5 6 ];>> C = 3.^A>> D = A.^B3、 已知⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t 若需要计算t ∈[-1,1],取间隔为0.01,试计算出相对应的y 值。
matlab实验指导答案详解(非常详细正确)
![matlab实验指导答案详解(非常详细正确)](https://img.taocdn.com/s3/m/e79999f8162ded630b1c59eef8c75fbfc77d9420.png)
matlab实验指导答案详解(⾮常详细正确)实验⼀ MATLAB ⼯作环境熟悉及简单命令的执⾏⼀、实验⽬的:熟悉MATLAB 的⼯作环境,学会使⽤MATLAB 进⾏⼀些简单的运算。
⼆、实验内容:MATLAB 的启动和退出,熟悉MATLAB 的桌⾯(Desktop ),包括菜单(Menu )、⼯具条(Toolbar )、命令窗⼝(Command Window)、历史命令窗⼝、⼯作空间(Workspace)等;完成⼀些基本的矩阵操作;学习使⽤在线帮助系统。
三、实验步骤:1、启动MATLAB ,熟悉MATLAB 的桌⾯。
2、在命令窗⼝执⾏命令完成以下运算,观察workspace 的变化,记录运算结果。
(1)(365-52?2-70)÷3 >>(365-52*2-70)/3 ans = 63.6667(2)>>area=pi*2.5^2 area = 19.6350(3)已知x=3,y=4,在MATLAB 中求z :()232y x y x z -= >>x=3 >>y=4>>z = x ^2 * y ^3 / (x - y) ^2 z = 576(4)将下⾯的矩阵赋值给变量m1,在workspace 中察看m1在内存中占⽤的字节数。
m1=11514412679810115133216 执⾏以下命令>>m1 =[16 2 3 13 ; 5 11 10 8 ; 9 7 6 12 ; 4 14 15 1 ] >>m1( 2 , 3 ) ans = 10 >>m1( 11 ) ans = 6>>m1( : , 3 ) ans =3 10 6 15>>m1( 2 : 3 , 1 : 3 ) ans =5 11 10 9 7 6>>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) ans = 34(5)执⾏命令>>help abs查看函数abs 的⽤法及⽤途,计算abs( 3 + 4i ) (6)执⾏命令>>x=0:0.1:6*pi; >>y=5*sin(x); >>plot(x,y)(6)运⾏MATLAB 的演⽰程序,>>demo ,以便对MATLAB 有⼀个总体了解。
MATLAB教程课后实验报告题目及解答[第一至第五章]
![MATLAB教程课后实验报告题目及解答[第一至第五章]](https://img.taocdn.com/s3/m/512b4436f111f18583d05a96.png)
MATLAB教程实验报告实验项目名称实验一 Matlab基本操作学生姓名汪德旺专业班级 09数教(1)班学号 0301090131实验成绩日期一. 实验目的和要求1、了解MATLAB 的开发环境。
2、熟悉Matlab的基本操作。
3、掌握建立矩阵的方法。
4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。
5、填写实验报告,实验报告文件取名为report1.doc。
6、于邮件附件形式将实验报告文件report1.doc 发到邮箱*******************,邮件主题为班级学号姓名,如:09数教1班15号张三。
二、实验内容1、先求下列表达式的值,然后显示MATLAB工作空间的使用情况并保存全部变量。
(1)0 122sin851ze =+(2)2212 1ln(0.4552i z x x+⎡⎤==⎢⎥-⎣⎦其中(3)0.30.330.3sin(0.3)ln,22a ae e az a--+=++a=-3.0,-2.9,-2.8,…, 2.8, 2.9,3.0(4)2242,011,12,0:0.5:2.521,23t tz t tt t t⎧≤<⎪=-≤<=⎨⎪-+≤<⎩其中t2.已知:1234413134787,2033657327 A B--⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值:(1)A+6*B和A-B+I(其中I为单位矩阵) (2)A*B和A.*B(3)A^3和A.^3(4)A/B和B\A(5)[A,B]和[A([1,3],:);B^2]3、设有矩阵A 和B123453016678910A=,B=17-691112131415023-41617181920970212223242541311⎡⎤⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦(1)求它们的乘积C 。
(2)将矩阵C 的右下角3*2子矩阵赋给D 。
(3)查看MATLAB 工作空间的使用情况。
915230-MATLAB第三版实验答案
![915230-MATLAB第三版实验答案](https://img.taocdn.com/s3/m/fc69726b6f1aff00bfd51e2e.png)
MATLAB程序设计与应用(第三版)实验参考答案实验1 MATLAB系统环境与运算基础%第1题(略)%第2题(略)%第3题%(1)z1=2*sin(85*pi/180)/(1+exp(2))%(2)x=12;y=1e-5;z2=pi/(x+y)%(3)x=[2,1+2i;-0.45,5];z2=0.5*log(x+sqrt(1+x.^2))%(3)a=-3.0:0.1:3.0;z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)%第4题A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7];A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]%第5题A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25] B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]C=A*B%求D方法1F=size(C)D=C(F(1)-2:F(1),F(2)-1:F(2))%求D方法2D=C(end-2:end,end-1:end)whos%第6题%(1)A=100:999;B=rem(A,21);C=length(find(B==0))%(2)A='lsdhKSDLKklsdkl';k=find(A>='A'&A<='Z');A(k)=[]实验2 MATLAB矩阵处理%第1题E=eye(3);R=rand(3,2);O=zeros(2,3);S=diag([2,3]);A=[E,R;O,S];A^2B=[E,(R+R*S);O,S^2]%第2题A=fix(10*rand(5))H=det(A)Trace=trace(A)Rank=rank(A)Norm=norm(A)%第3题H=hilb(5)P=pascal(5)Hh=det(H)Hp=det(P)Th=cond(H)Tp=cond(P)%第4题A=[-29,6,18;20,5,12;-8,8,5][V,D]=eig(A)%数学意义略%第5题方法一%(1):A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; b=[0.95,0.67,0.52]';x=inv(A)*b%(2):B=[0.95,0.67,0.53]';x=inv(A)*B%(3):cond(A)%第5题方法二A=hilb(4)A(:,1)=[]A(4,:)=[]B=[0.95,0.67,0.52]';X=inv(A)*BB1=[0.95,0.67,0.53]';X1=inv(A)*B1N=cond(B)N1=cond(B1)Na=cond(A) %矩阵A为病态矩阵实验3 顺序结构程序设计%程序(略)实验4 选择结构程序设计%第1题程序一x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0];y=[]; %建立存放所有y值的矩阵for x0=xif x0<0&x0~=-3y=[y,x0*x0+x0-6];elseif x0>=0&x0<5&x0~=2&x0~=3y=[y,x0*x0-5*x0+6];elsey=[y,x0*x0-x0-1];endendx %输出所有xy %输出所有y%第1题程序二x=[-5,-3,1,2,2.5,3,5];y=[];for a=1:7if x(a)<0&x(a)~=-3y=[y,(x(a))^2+x(a)-6];elseif x(a)>=0&x(a)<5&x(a)~=2&x(a)~=3y=[y,(x(a))^2-5*x(a)+6];elsey=[y,x(a)*x(a)-x(a)-1];endend%第2题程序一x=input('请输入一个百分制成绩:');if x>100|x<0disp('您输入的成绩不是百分制成绩,请重新输入。
MATLAB实验报告3
![MATLAB实验报告3](https://img.taocdn.com/s3/m/c08d06fe770bf78a652954ab.png)
3.6 基于Simulink 控制系统的稳态误差分析1.实验目的1)掌握使用Simulink 仿真环境进行控制系统稳态误差分析的方法。
2)了解稳态误差分析的前提条件是系统处于稳定状态。
3)研究系统在不同典型信号输入作用下,稳态误差的变化。
2.实验内容(1)研究系统在不同典型输入信号作用下,稳态误差的变化。
【例3-11】 已知一个单位负反馈系统开环传递函数为G(s)=)11.0(10 s s K ,分别K=10和K=1时,系统单位阶跃响应曲线并求出单位阶跃响应稳态误差。
【解】 首先对闭环系统判稳。
该系统为零极点模型,用函数roots()命令判断系统闭环全部特征根的实部都是负值,说明闭环系统稳定。
这样进行稳态误差分析才有意义。
K=10时的判稳程序如下:>> n1=100;d1=conv([1,0],[0.1,1]);G=tf(n1,d1);sys=feedback(G ,1);roots(sys.den{1})ans =-5.0000 +31.2250i-5.0000 -31.2250i然后在Simulink 环境下,建立系统数学模型,如图3-21所示。
设置仿真参数并运行,观察示波器Scope 中系统的单位阶跃响应曲线,如图3-22所示,并读出单位阶跃响应稳态误差。
图3-21 基于Simulink Ⅰ型控制系统单位阶跃响应(K=10)结构图图3-22 基于Simulink Ⅰ型控制系统单位阶跃响应稳态误差曲线【分析】实验曲线表明,Ⅰ型单位反馈系统在单位阶跃输入作用下,稳态误差e ssr=0,即Ⅰ型单位反馈系统稳态时能完全跟踪阶跃输入,是一阶无静差系统。
K=1时的判稳程序如下:>> n1=10;d1=conv([1,0],[0.1,1]);G=tf(n1,d1);sys=feedback(G,1);roots(sys.den{1})ans =-5.0000 + 8.6603i-5.0000 - 8.6603i在Simulink环境下建立的数学模型及仿真参数运行后示波器Scope中系统的单位阶跃响应曲线如下图所示。
金融MATLAB实验报告三解析
![金融MATLAB实验报告三解析](https://img.taocdn.com/s3/m/a342cf14d5bbfd0a7856735e.png)
安徽财经大学金融证券实验室实验报告实验课程名称《金融MATLAB 》开课系部_______________ 金融学院____________班级______________________________学号_____________________________姓名_____________________________指导教师____________________________________2015年**月**日实验分析过程一、收益、风险和有效前沿的计算从Wind咨询金融终端分别下载三只股票(美好集团、石油化服和首开股份)从2013年年初至今的日收盘价价格,经过相关处理得出三只股票的收益率均值、标准差以及协方差矩阵等数据,如下表。
现根据表格数据进行关于收益、风险和有效前沿计算。
PortRetur n =0.00060.00080.00090.00100.00120.00130.00140.00150.00170.0018PortWts =0 0.0439 0.95610 0.1724 0.82760.0275 0.2678 0.70470.1032 0.3054 0.59140.1789 0.3430 0.47810.2546 0.3806 0.36470.3303 0.4183 0.25140.4060 0.4559 0.13810.4817 0.4935 0.02481.0000 0 0画图:E3罡山frontcon (ExpRetur n,ExpCovaria nee, NumPorts)0.013 0.014 O.O'fS 0.016 0.017 0.018 0.019 0 02 0.021 0 022 0.C23-'' "■ i ■ I 4,;i::i : . (2)例2中如果各个资产投资上限为50%,求解有效前沿?>> ExpReturn=[0.0018,0.0016,0.0006];ExpCovariia nce=[0.0010,0.0004,0.0005;0.0004,0.0017,0.0003;0.0005,0.0003,0.0013];NumPorts =10;AssetBo un ds=[0,0,0;0.5,0.5,0.5];[PortRisk, PortReturn, PortWts] = frontcon (ExpReturn,ExpCovaria nee, NumPorts,[],AssetBou nds) PortRisk =0.01560.01600.01640.01680.01730.01770.01820.01870.01920.0197PortRetur n =0.00110.00120.00130.00130.00140.00140.00150.00160.00160.0017PortWts =0.1196 0.3804 0.50000.1967 0.3519 0.45140.2344 0.3706 0.39490.2722 0.3894 0.33850.3099 0.4081 0.28200.3476 0.4268 0.22560.3853 0.4456 0.16910.4230 0.4643 0.11260.4608 0.4831 0.05620.5000 0.5000 0frontcon (ExpRetur n,ExpCovaria nee, NumPorts,[],AssetBo unds)Risk (Standard Deviation)3•约束条件下有效前沿例3.例配置美好集团、石油化服、首开股份三个资产,美好集团最大配置60%,石油化服最大配置70% , 首开股份最大配置50%,美好集团为资产集合A,石油化服、首开股份组成资产计划B,集合A的最大配置为70%,集合B的最大配置为50%,集合A的配置不能超过集合B的3倍,则如何配置?解:>> NumAssets = 3;ExpReturn=[0.0018,0.0016,0.0006];ExpCovariia nce=[0.0010,0.0004,0.0005;0.0004,0.0017,0.0003;0.0005,0.0003,0.0013];NumPorts =5;PVal = 1;AssetMin = 0;AssetMax=[0.6,0.7,0.5];GroupA = [1 0 0]; GroupB = [0 1 1];GroupMax =[0.7,0.5];AtoBmax = 3;Co nSet = portco ns('PortValue', PV al, NumAssets,'AssetLims',...AssetMin, AssetMax, NumAssets, 'GroupCompariso n',GroupA, NaN,...AtoBmax, GroupB,GroupMax );[PortRisk, PortReturn, PortWts] = portopt(ExpReturn, ExpCovaria nee, NumPorts, [], Co nSet) PortRisk =0.01560.01650.01760.01870.0199PortRetur n =0.00110.00130.00140.00160.0017PortWts =0.1196 0.3804 0.50000.2468 0.3768 0.37640.3346 0.4204 0.24500.4225 0.4640 0.11350.6000 0.4000 0画图>> portopt(ExpRetur n, ExpCovaria nce,N umPorts, [], Con Set)Mean-Variance-Efficient Frontierx 10'31.8171.6151.41.31.21.10.0155 0.016 00165 0.017 0.0175 0018 0.0185 0G19 C.0195 0.02Risk (Standard Deviation)解:1.读取数据>> [typ, desc, fmt] = xlsfinfo('HHHZ.xlsX'); >> [data,textdate]= xlsread('HHHZ.xlsX'); >> HHHZ=data; >> save HHHZ HHHZ 2.三只股票和一只指数的精致曲线>> load HHHZ >> figure;>> hold on >> plot(HHHZ(:,1)/HHHZ(1,1),'k') >> plot(HHHZ(:,2)/HHHZ(1,2),'r-.') >> plot(HHHZ(:,3)/HHHZ(1,3),'bo')>> plot(HHHZ(:,4)/HHHZ(1,4),'g--') >> xlabel('time') >> ylabel('price')>> legend('中证800','华丽家族','华远地产','合肥城建')8---- 中证日00 7 .------- 华丽家族3. Beta 与 Alpha 计算、投资组合绩效分析从Wind 咨询金融终端分别下载三只股票(华丽家族、华远地产、合肥城建)和一只指数(中证800 )从2014年年初至今的日收盘价数据,加以处理整理。
Matlab实验三答案
![Matlab实验三答案](https://img.taocdn.com/s3/m/6413beee09a1284ac850ad02de80d4d8d15a01bc.png)
Matlab实验三答案实验三 MATLAB语言的程序设计一、实验目的:1、熟悉MATLAB程序编辑与设计环境2、掌握各种编程语句语法规则及程序设计方法3、函数文件的编写和设计4、了解和熟悉跨空间变量传递和赋值二、实验基本知识:1、程序流程控制语句for循环结构语法:for i=初值:增量:终值语句1……语句nend说明:1.i=初值:终值,则增量为1。
2.初值、增量、终值可正可负,可以是整数,也可以是小数,只须符合数学逻辑。
while 循环结构语法:while 逻辑表达式循环体语句end说明:1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。
若表达式的值为真,执行循环体语句一次、在反复执行时,每次都要进行判断。
若表达式的值为假,则程序执行end之后的语句。
2、为了避免因逻辑上的失误,而陷入死循环,建议在循环体语句的适当位置加break语句、以便程序能正常执行。
(执行循环体的次数不确定;每一次执行循环体后,一定会改变while 后面所跟关系式的值。
)3、while循环也可以嵌套、其结构如下:while逻辑表达式1循环体语句1while逻辑表达式2循环体语句2end循环体语句3endif-else-end分支结构if 表达式1语句1else if 表达式2(可选)语句2else(可选)语句3endend说明:1.if结构是一个条件分支语句,若满足表达式的条件,则往下执行;若不满足,则跳出if结构。
2.else if表达式2与else为可选项,这两条语句可依据具体情况取舍。
3.注意:每一个if都对应一个end,即有几个if,记就应有几个end。
switch-case结构语法:switch表达式case常量表达式1语句组1case常量表达式2语句组2……otherwise语句组nend说明:1.switch后面的表达式可以是任何类型,如数字、字符串等。
2.当表达式的值与case后面常量表达式的值相等时,就执行这个case后面的语句组如果所有的常量表达式的值都与这个表达式的值不相等时,则执行otherwise后的执行语句。
参考答案Matlab实验报告
![参考答案Matlab实验报告](https://img.taocdn.com/s3/m/51a2dcc2aef8941ea76e05da.png)
»k=input('’,’s’);Eie48458DHUEI4778»f=find(k>=’A’&k<=’Z’);f=9 10 11 12 13»k(f)=[ ]K=eie484584778●输入矩阵,并找出其中大于或等于5的元素。
(find)1.输入命令:»h=[4 8 10;3 6 9; 5 7 3];»[i,j]=find(h>=5)i=3 j=11 22 23 21 32 3●不采用循环的形式求出和式的数值解。
(sum)1.输入命令:»w=1:63;»q=sum(2.^w)q=1.8447e+0191.打开Matlab,新建M文件2.输入命令:y=0;i=1while 1f=1/(2i-1);y=y+fif y>5break;endi=i+1endn=i-1y=y-f3.保存,命名为eg1,并在命令窗口输入eg1:»eg1n=3094y=4.9999编程完成,对输入的函数的百分制成绩进行等绩转换,90~100为优,80~89为良,70~79为中,60~69为及格。
1.打开Matlab,新建M文件2.输入命令:1.打开Matlab ,新建M 函数文件2.输入命令:function y=hs(x)y=2*sqrt(x.^2+10)-53.保存,命名为hs,并在命令窗口输入»hs(12)y=19.8193ans=19.8193 编程求分段函数,并求2226;0356;05321;x x x x y x x x x x x x +-<≠=-+≤<≠≠-+且且及其它输入x=[-5.0,-3.0,1.0,2.0,2.5,3.0,3.5]时的输出y 。
1.打开Matlab,新建M 函数文件2.输入命令:function y=zhi(x)for i=1:length(x)if x(i)<0&x(i)~=-3y(i)=x(i).^2+x(i)-6;elseif x(i)>=0&x(i)<5&x(i)~=3&x(i)~2y(i)=x(i).^2-x(i)-5*x(i)+6;else y(i)=x(i).^2-x(i)+1;endendy3.保存,保存名为zhi,并在命令窗口输入»x=[-5,-3,1,2,2.5,3,3.5];»zhi(x)y=Columns 1 through 614.0000 13.0000 2.0000 3.0000 -0.2500 7.0000Column 70.7500ans=Columns 1 through 614.0000 13.0000 2.0000 3.0000 -0.2500 7.0000Column 70.7500实验三 Matlab的基本绘图命令一、实验目的:1.熟悉二维、三维绘图函数的基本命令。
金融MATLAB实验报告三解析
![金融MATLAB实验报告三解析](https://img.taocdn.com/s3/m/4658670f5a8102d276a22f7f.png)
目
的
凸度计算、利率的期限结构、技术指标计算等。
使用 MATLAB金融工具箱对下述 6 个主题进行数量分析,数据来源自行在网上搜
寻,要求是 2014 年之后的数据。(可参照各章的例题)
1. 期权定价分析 (第 10 章)
2. 收益、风险和有效前沿的计算 (第 12 章)
3. 投资组合绩效分析 (第 13 章)
2013-2014 年的数据,对华丽家族、华远地
日期 2013-01-04 2013-01-07 2013-01-08 2013-01-09 2013-01-10 2013-01-11 2013-01-14 2013-01-15
中证 800 收盘价 ( 元)
2643.83 2658.72 2657.14 2659.23 2666.74 2612.77
( 2)Alpha 计算 以中证 800 作为市场收益,使用 portalpha 函数分别计算华丽家族、华远地产、合肥城建的
beta 值。过程
如下:
>> load HHHZ
Rate=price2ret(HHHZ);
zz800=Rate(:,1);
hljz=Rate(:,2);
hydc=Rate(:,3);
0.0004,0.0017,0.0003; 0.0005,0.0003,0.0013]; NumPorts =10; AssetBounds=[0,0,0;0.5,0.5,0.5]; [PortRisk, PortReturn, PortWts] = frontcon(ExpReturn,ExpCovariance, NumPorts,[],AssetBounds) PortRisk =
0.2450 0.1135
金融MATLAB实验报告三复习资料详解
![金融MATLAB实验报告三复习资料详解](https://img.taocdn.com/s3/m/56fdbfdd2e3f5727a4e96258.png)
安徽财经大学金融证券实验室实验报告实验课程名称《金融MATLAB 》开课系部金融学院班级学号姓名指导教师年月日3.计算期权Delta。
例2.假设欧式股票期权,六个月后到期,执行价格90元,现价为102元,无股利支付,股价年化波动率为55%,无风险利率为8%,计算期权Delta。
解:clearPrice=102;>> Strike=90;>> Rate=0.08;>> Time=6/12;>> V olatility=0.55;[CallDelta, PutDelta] = blsdelta(Price, Strike, Rate, Time, Volatility)计算结果:CallDelta =0.7321PutDelta =-0.26794.利用不同的Price与Time计算Detla三维关系。
>> Price=60:1:102;>> Strike=90;Rate=0.08;>> Time=(1:1:12)/12;>> Volatility=0.55;>> [Price,Time]=meshgrid(Price,Time);[Calldelta, Putdelta] = blsdelta(Price, Strike, Rate, Time, Volatility);>> mesh(Price, Time, Putdelta);xlabel('Stock Price ');ylabel('Time (year)');zlabel('Delta');>>5.B-S公式隐含波动率计算例3:假设欧式股票期权,一年后,执行价格99元,现价为105元,无股利支付,股价年化波动率为40%,无风险利率为10%,则期权价格为:解:clear>> Price=105;>> Strike=99;>> Rate=0.1;>> Time=1;>> CallValue=15;>> CallVolatility = blsimpv(Price, Strike, Rate, Time, CallValue, [], [], [],。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
安徽财经大学金融证券实验室实验报告
实验课程名称《金融》TLABMA
金融学院部课系开
级班
学号
姓名
师导指教日年月
1.
2
一、期权定价分析
1.black-scholes方程求解
例1:假设欧式股票期权,六个月后到期,执行价格90元,现价为102元,无股利支付,
股价年化波动率为55%,无风险利率为8%,计算期权价格。
解:clear
Price=102;
>>Strike=90;
>>Rate=0.08;
>>Time=6/12;
>>V olatility=0.55;
[CallDelta,PutDelta]=blsprice(Price,Strike,Rate,Time,V olatility)
计算结果:
CallDelta=
23.5648
PutDelta=
8.0358
2.期权价格与波动率关系分析
Price=102;
>>Strike=90;
>>Rate=0.08;
>>Time=6/12;
V olatility=0.08:0.01:0.5;
>>N=length(V olatility)
Call=zeros(1,N);
Put=zeros(1,N);
for i=1:N
[Call(i),Put(i)]=blsprice(Price,Strike,Rate,Time,V olatility(i));
N=
43
end
plot(Call,'b--');
hold on
plot(Put,'b');
xlabel('V olatility')
ylabel('price')
legend('Call','Put')
3.计算期权Delta。
例2.假设欧式股票期权,六个月后到期,执行价格90元,现价为102元,无股利支付,
股价年化波动率为55%,无风险利率为8%,计算期权Delta。
解:clear
Price=102;
>>Strike=90;
>>Rate=0.08;
>>Time=6/12;
>>V olatility=0.55;
[CallDelta,PutDelta]=blsdelta(Price,Strike,Rate,Time,Volatility)计算结果:CallDelta=
0.7321
PutDelta=
-0.2679
4.利用不同的Price与Time计算Detla三维关系。
>>Price=60:1:102;
>>Strike=90;
Rate=0.08;
4.
5.B-S公式隐含波动率计算
例3:假设欧式股票期权,一年后,执行价格99元,现价为105元,无股利支付,股价年
化波动率为40%,无风险利率为10%,则期权价格为:
解:clear
>>Price=105;
>>Strike=99;
>>Rate=0.1;
>>Time=1;
>>CallValue=15;
>> CallVolatility = blsimpv(Price, Strike, Rate, Time, CallValue, [], [], [],
5.
{'Call'})
计算结果:
CallVolatility=
NaN
>>PutValue=7;
>>PutVolatility=blsimpv(Price,Strike,Rate,Time,PutValue,[],[], [],
{'Put'})
PutVolatility=
0.3455
6.期权二叉树模型的计算
例:假设欧式股票期权,三个月后到期,执行价格85元,现价为95元,无股利支付,股价
年化波动率为60%,无风险利率为10%。
解:clear
>>Price=95;
>>Strike=85;
>>Rate=0.1;
>>Time=4/12;
>>flag=1;
>>Increment=1/12;
>>Volatility=0.6;
>>[AssetPrice,OptionValue]=binprice(Price,Strike,Rate,Time, Increment,
Volatility,flag)
计算结果:
AssetPrice=
95.0000112.9654134.3283159.7312189.9379
079.891795.0000112.9654134.3283
67.18610079.891795.0000
6.
00056.501267.1861
47.51550000
OptionValue=
20.057432.495250.733375.4365104.9379 28.67088.860849.3283016.1268
10.00004.7685002.2739
00000
00000
7.
成绩:
指导教师签字:8.。