数学模型实验题

这个没有什么排版，题目我放在最后面了，内容大概吧都没什么大问题，以后发现会改的。

16题是一个重点哦！！！！要多看看

---

1.

f='exp(x^2)'    f=sym('exp(x^2)')
eqt='a*(x-1)^2+b*(x-1)+c=0'
difeq='sin(x)+x=0'

---

2.

>> syms x y
>> expand(sin(x-y))
 
ans =
 
sin(x)*cos(y)-cos(x)*sin(y)

---

3.

>> digits(3)
>> vpa(sqrt(2))
 
ans =
 
1.41

>> digits(5)
>> vpa(sqrt(2))
 
ans =
 
1.4142

>> digits(9)
>> vpa(sqrt(2))
 
ans =
 
1.41421356

---

4.

>> digits(32)
>> a=vpa(sqrt(2))
 
a =
 
1.4142135623730950488016887242097
 
 
>> b=vpa('sqrt(2)')
 
b =
 
1.4142135623730950488016887242097

此处显示运行结果都一样。

---

5.

>> syms x y
>> f=exp(-y)
 
f =
 
exp(-y)
 
 
>> g=x^(-2)
 
g =
 
1/x^2
 
 
>> t=compose(f,g)
 
t =
 
exp(-1/x^2)
 
 
>> diff(t)
 
ans =
 
2/x^3*exp(-1/x^2)

---

6.

>> syms x t
>> f=log((2+x)/(2-x))
 
f =
 
log((2+x)/(2-x))
 
 
>> f=subs(f,x,t^2)
 
f =
 
log((2+t^2)/(2-t^2))
 
 
>> subs(f,t,1)

ans =

    1.0986

---

7.

>> syms x
>> f=x^3+x
 
f =
 
x^3+x
 
 
>> finverse(f,x)

Warning: finverse(x^3+x) is not unique.
> In sym.finverse at 43
 
ans =
 
1/6*(108*x+12*(12+81*x^2)^(1/2))^(1/3)-2/(108*x+12*(12+81*x^2)^(1/2))^(1/3)

---

8.

>> factor(sym('32784'))
 
ans =
 
(2)^4*(3)*(683)
 
 
>> factor(sym('26788312456987'))
 
ans =
 
(53)*(277)*(1824692627)

---

9.

>> syms x
>> factor(x^3+2*x^2+2*x+1)
 
ans =
 
(x+1)*(x^2+x+1)
 
 
>> factor(x^6-1)
 
ans =
 
(x-1)*(x+1)*(x^2+x+1)*(x^2-x+1)

---

10.

>> syms x
>> limit((3^x+9^x)^(1/x),inf)
 
ans =
 
9

>> syms x y
>> limit(limit(x*y/(sqrt(x*y+1)-1),x,0),y,0)
 
ans =
 
2

---

11.

>> syms k n
>> symsum(k^2,1,n)
 
ans =
 
1/3*(n+1)^3-1/2*(n+1)^2+1/6*n+1/6

>> symsum(1/k^2,1,inf)
 
ans =
 
1/6*pi^2

---

12.

>> syms x y
>> f=(x-y)^3*sin(x+2*y)
 
f =
 
(x-y)^3*sin(x+2*y)

>> int(f,y,-x,x)
 
ans =
 
3/8*sin(3*x)+4*x^3*cos(x)-3*x^2*sin(x)-3/2*cos(x)*x+3/8*sin(x)

---

13.

>> syms x a
>> f=x^2/sqrt(a^2-x^2)
 
f =
 
x^2/(a^2-x^2)^(1/2)
 
 
>> int(f,x)
 
ans =
 
-1/2*x*(a^2-x^2)^(1/2)+1/2*a^2*atan(x/(a^2-x^2)^(1/2))

---

14.

>> syms x
>> f=x^(-x)
 
f =
 
x^(-x)
 
 
>> int(f,x,0,1)
Warning: Explicit integral could not be found.
> In sym.int at 58
 
ans =
 
int(x^(-x),x = 0 .. 1)

---

15.

符号求解

>> syms x1 x2 x3
>> [x1,x2,x3]=solve('2*x1+2*x2+x3=2','x2+2*x3=1','x1+x2+3*x3=3')
 
x1 =
 
6/5
 
x2 =
 
-3/5

x3 =
 
4/5

数值计算

>> A=[2 2 1;0 1 2;1 1 3]

A =

     2     2     1
     0     1     2
     1     1     3

>> b=[2;1;3]

b =

     2
     1
     3

>> A\b

ans =

    1.2000
   -0.6000
    0.8000

---

16.

不要使用solve，错误。。。

>> solve('cos(2*x)+sin(x)-1=0')
ans =
     pi
      0
 1/6*pi
 5/6*pi

>> solve('sin(x)=0')
ans =0

以上都是错的，看图像不只有4个解，这是一个对称函数。上面这个例子也是证明

正确的

这个代码写在一个myfun的m文件中，并保存，之后来调用。

function f=myfun(x)
f=cos(2*x)+sin(x)-1

这个写在主页面中，或者另外一个不重名的m文件中都可以运行出来。 

x=fsolve(@myfun,[-2 2])

 最后运行的结果贴在这里了，这个是最终的结果，之前还有很多解，这个是最好的。

Optimization terminated: first-order optimality is less than options.TolFun.

x =

   -3.1416    2.6180

---

17.

>> syms x y
>> dsolve('Dy+6*y=exp(x)','x')
 
ans =
 
1/7*exp(x)+exp(-6*x)*C1

---

18.

>> syms x t
>> dsolve('2*D2x-5*Dx-3*x=45*exp(2*t)','x(0)=2,Dx(0)=1','t')
 
ans =
 
4*exp(-1/2*t)+7*exp(3*t)-9*exp(2*t)

>> t=0:0.1:2;
>> plot(t,4*exp(-1/2*t)+7*exp(3*t)-9*exp(2*t))

---

19.

>> [x,y]=dsolve('Dx=y','Dy=-x')
 
x =
 
-C1*cos(t)+C2*sin(t)
 
y =
 
C1*sin(t)+C2*cos(t)

---

20.

>> syms x y
>> f=3*x^3-4*x^4-5
f =
3*x^3-4*x^4-5

>> diff(f)
ans =
9*x^2-16*x^3
>> solve(ans)
ans = 
    0
    0
 9/16

ans=9/16即为极值对应的x值

---

 

 

---

x=[-1:0.1:2];
f=3*x.^3-4*x.^4-5;
plot(x,f)

[x2,maxf]=fminbnd('-3*x^3+4*x^4+5',-1,2)%极小值
maxf=-maxf

[x2,maxf]=fminunc('-3*x^3+4*x^4+5',-1,2)%最小值
maxf=-maxf

---