MATLAB程序设计常微分方程编程例题答案数学.docVIP

matlab程序设计-常微分方程 1、功能：计算函数在指定点的值 实验：求，，， M 文件 function f=fd(x) if x=-1 f=sqrt(-x); elseif x0 f=int(((sin(t))^(10)),0,x); else x=0 f=3*x+1; end 结果 fd(-2) ans = 1.4142 fd(-1) ans = 1 fd(0) ans = 1 ans = 1 fd(pi) ans = 63/256*pi 2、实验：分别计算椭圆与的焦距与离心率 M文件： function [c,e]=fun(a,b) c=sqrt(a^2-b^2); e=c/a; end 结果：  [c,e]=fun2(5,3) c = 4 e = 0.8000 [c,e]=fun2(1,sqrt(2)) c = 0 + 1.0000i e = 0 + 1.0000i 3、实验：分别绘制函数与的图像（要求：图像光滑） M文件： function ea(a,b,c) k=2*pi/b; t=0:0.0001:2*k; end y=a*sin(b*t+c) plot(t,y) 结果： ea(1,2,0) Columns 62793 through 62800 -0.0080 -0.0078 -0.0076 -0.0074 -0.0072 -0.0070 -0.0068 -0.0066 Columns 62801 through 62808 -0.0064 -0.0062 -0.0060 -0.0058 -0.0056 -0.0054 -0.0052 -0.0050 Columns 62809 through 62816 -0.0048 -0.0046 -0.0044 -0.0042 -0.0040 -0.0038 -0.0036 -0.0034 Columns 62817 through 62824 -0.0032 -0.0030 -0.0028 -0.0026 -0.0024 -0.0022 -0.0020 -0.0018 Columns 62825 through 62832 -0.0016 -0.0014 -0.0012 -0.0010 -0.0008 -0.0006 -0.0004 -0.0002 ea(10,500*pi,pi/6) y = Columns 1 through 8 5.0000 6.2932 7.4314 8.3867 9.1355 9.6593 9.9452 9.9863 Columns 9 through 16 9.7815 9.3358 8.6603 7.7715 6.6913 5.4464 4.0674 2.5882 Columns 17 through 24 1.0453 -0.5234 -2.0791 -3.5837 -5.0000 -6.2932 -7.4314 -8.3867 Columns 25 through 32 4、实验：求与 M文件 ： function s=qh(n) s=1; for i=2:n s=s*i; end 结果 ： JC(18) ans = 6.4024e+015 JC(0) ans = 1 5、功能：求最大的，使得 实验：分别取，，，求 M文件 function n=maxn(c) s=1 i=2 while sc s=s+(1/i); i=i+1; end n=i-2; 结果： maxn(1.5) s = 1 i = 2 ans = 1 maxn(1.50001) s = 1 i = 2 ans = 2 maxn(12) s = 1 i = 2 ans = 91379 6、功能：找出数组中最小元素所在的位置 例如，数组，其中最小元素的位置为 实验：找出数组中最小元素所在的位置 找出数组中最小元素所在的位置 M文件 function?[t,i]=zb(a)n=l