（3）c = 1:2:5 （4）c = [ 1, 2] [ 3, 4]
⑤已知符号表达式f?1?sin2x，g?2x?1，计算x＝0.5时，f的值；计算复合
函数f(g(x))。 源程序：
>> syms x; >> f=1-sin(x)^2 >> g=2*x+1;
>> f1=subs(f,'0.5') >> fg=compose(f,g)
运行结果：
f = 1-sin(x)^2 f1 =
1-sin((0.5))^2 fg =
1-sin(2*x+1)^2
⑥求 limx2?1x?2x2?3x?1。
源程序：
>> syms x;
>> ans=limit((x^2-1)/(x^2-3*x+1),x,2)
运行结果：
ans = -3
⑦求函数 f(x)= cos 2x -sin 2x的积分；求函数g(x)?ex?xsinx的导数。源程序：
（1）>> syms x y z a b >> S=cos(2*x)-sin(2*x);
16 / 24
>> int(S)
（2）>> syms x y t u v z a b >> S=sqrt(exp(x)+x*sin(x));
>> diff(S) 运行结果：
（1）ans =
1/2*sin(2*x)+1/2*cos(2*x) （2） ans =
1/2/(exp(x)+x*sin(x))^(1/2)*(exp(x)+sin(x)+x*cos(x))
?⑧计算定积分?60(sinx?2)dx
源程序：
>> syms x y z a b >> S=sin(x)+2; >> int(S,0,pi/6)
运行结果：
ans =
-1/2*3^(1/2)+1/3*pi+1
⑨求下列线性代数方程组的解。
??x?y?z?10?3x?2y?z?14 ??2x?3y?z?1源程序：
>> [x,y,z]=solve('x+y+z-10','3*x+2*y+z-14','2*x+3*y-z-1')
运行结果：
x = 1 y = 2
17 / 24
z = 7
⑩求解当y(0)=2，z(0)=7时，微分方程组的解。
?dy?dx?z?sinx ?dz??y?1?x?dx源程序：
>> [y,z]=dsolve('Dy-z=sin(x)','Dz+y=1+x','y(0)=2,x(0)=7','x'); Warning: Explicit solution could not be found.
> In C:\\MATLAB6p5\\toolbox\\symbolic\\dsolve.m at line 326
>> [y,z]=dsolve('Dy-z=sin(x)','Dz+y=1+x','y(0)=2,z(0)=7','x'); >> pretty(y) >> pretty(z)
运行结果：
cos(x) + 6 sin(x) + 1/2 x sin(x) + 1 + x - 3/2 sin(x) + 6 cos(x) + 1 + 1/2 cos(x) x
4 MATLAB程序设计
实验目的：掌握MATLAB程序设计的主要方法，熟练编写MATLAB函数。 实验内容：
① 用π /4≈1-1/3+1/5-1/7+…公式求π的近似值，直到最后一项的绝对值小于10-6 为止，试编写其M脚本文件。 源程序：
>> t=1; pi=0; n=1; s=1;
>> while abs(t)>1e-7 pi=pi+t; n=n+2;
18 / 24
s=-s; t=s/n; end >>
>> pi=pi*4;
>> fprintf('pi=%f/n',pi)
运行结果：
pi=3.141592/n
②分别用for和while结构计算1+22+33+…+100100的运行程序。源程序： >> clear sum=0; for i=1:100 sum=sum+i^i; end >> sum
>> clear all; sum=0; i=1;
while i<=100 sum=sum+i^i; i=i+1; end
>> sum
执行： sum =
1.0037e+200
19 / 24
sum =
1.0037e+200
5 MATLAB绘图
实验目的：掌握MATLAB二维图形绘制，掌握图形属性的设置和图形修饰；掌握图像文件的读取和显示。 实验内容：①绘制图形
图1
其中x初值为0，终值为12.5，步长为0.1，第一个曲线y1为sin(?x)，第二
6个曲线y2为cos(?x)。写出图1的绘制源程序。按照以下的步骤进行(1)产生曲线
6的数据(共有3组数据：x,y1,y2)；(2)选择合适的线形、标记、颜色(正弦曲线为红色，余弦曲线为紫色)；(3)添加图例及文字说明信息；(4)添加坐标轴说明与图标题。 源程序：
>> x=0:0.1:12.5; >> y1=sin(pi*x/6); >> y2=cos(pi*x/6);
>> plot(x,y1,'r--',x,y2,'p-');
20 / 24