上机1：
1. 课后习题2，3，8
习题2：1.5e3
1.5e3=1.5×103）
习题3：合法的变量名是：xyz_3,abcdefgh
变量名由英文、数字和下连符三种字符组成；变量名的第一个字符必须是英文字母；变量名最多可包含63个字符。）
习题8：第一个指令能实现目的，第二个和第三个不能。
2. a=-8;p=[1,0,0,-a];r=roots(p)
r= -2.0000
1.0000+1.7321i
1.0000-1.7321i
书p10例1.3-6
令，，则，。求解方程。
指令r=roots(p)，求多项式的根。先构造一个多项式，p=[1,0,0,-a]是多项式的系数分量，r=roots(p)求多项式的根。）
给出数值型符号结果（可参考课本2.2.2节，）
a=sqrt(3)+pi
sa=sym('sqrt(3)+pi')
ca=class(a)
csa=class(sa)
vpa(sa-a)a =
sa =sqrt(3)+pi
ca =double
csa =sym
ans =6669047319828670e-15 %例2-1-1
上机2：习题2-9
clear all
syms x y
f2=int(int((x^2+y^2),y,1,x^2),x,1,2)
vf2=vpa(f2)f2 =
vf2 =
习题10clear all
syms t x
f=sin(t)/t;
y=int(f,t,0,x);
ezplot(y,[0,2*pi])
subs(y,x,4.5)ans =
习题21dsolve('dx=a*t^2+b*t','x(0)=2')或dsolve('dx=a*t^2+b*t','x(0)=2','t')
ans =1/3*a*t^3+1/2*b*t^2+2
习题22s=dsolve('df-3*f-4*g=0,dg-3*g+4*f=0','f(0)=0,g(0)=1','x')
s = f: [1x1 sym]
g: [1x1 sym]
disp([
exp(3*x)*sin(4*x), exp(3*x)*cos(4*x)]
上机3：1 指令1：x=linspace(0,2*pi,10)
指令2：x=0:2*pi/9:2*pi
指令3：x=[0 2*pi/9 4*pi/9 6*pi/9 8*pi/9 10*pi/9 12*pi/9 14*pi/9 16*pi/9 18*pi/9]
指令4：x=[0,2*pi/9,4*pi/9,6*pi/9,8*pi/9,10*pi/9,12*pi/9,14*pi/9,16*pi/9,18*pi/9]
x =columns 1 through 10
rand('state',0)
a=rand(3,5);
r=(a>0.5);
ri,ci]=find(r) %全下标。
rj=find(r单下标。ri =
ci =
rj =
数组运算法：
cleart=0:0.1:10;
y=1-exp(-0.5.*t).*cos(2.*t);
plot(t,y)
xlabel('t'),ylabel('y')
标量循环运算法：
cleart=0:0.1:10;
l=length(t);
for k=1:l
y(k)=1-exp(-0.5*t(k))*cos(2*t(k));
endplot(t,y)
xlabel('t'),ylabel('y')
上机4：课后习题17：
1）求“商”及“余”多项式。
clear all
format rat %为避免浮点显示。
n=conv([3,0,1,0],[1,0,0,0.5]);计算分子多项式。
d=conv([1,2,-2],[5,2,0,1]);计算分母多项式。
q,r]=deconv(n,d求商多项式和余多项式。
cq='商多项式q(x)为 ';cr='余多项式r(x)为 ';
disp([cq,poly2str(q,'x')]disp([cr,poly2str(r,'s')]
商多项式q(x)为 0.6 x - 1.44
余多项式r(x)为 -4.4409e-016 s^6 + 8.8818e-016 s^5 + 21.
88 s^4 - 5.34 s^3 - 5.52 s^2 + 4.
58 s - 2.88
2）利用计算所得“商”和“余”验算分子多项式。
qd=conv(q,d计算“商”与“分母”的乘积。
nn=qd+r重复得到的分子多项式。
norm(nn-n) %采用计算矩阵范数的方式比较两个矩阵之间的“差”ans =
作业19：clear all
randn('state',1);
u=2*(randn(1,100)>0.5)-1;
ku=0:99;%u有100个元素。
h=[0.05,0.24,0.40,0.24,0.15,-0.1,0.1];
s=conv(u,h);
ns1=0+1;ns2=99+size(h,2);%ns1为s的左端点，ns2为s的右端点。
ns=ns1:ns2;
subplot(2,1,1),stem(ku,u),text(50,0.5,'系统输入信号')
subplot(2,1,2),stem(ns,s),text(50,1,'系统输出信号')
xlabel('n')
上机5：习题1：
clear all
t=0:2*pi/200:2*pi;
x=4*cos(t);y=2*sin(t);
plot(x,y,'r.',markersize',20)
axis([-4,4,-pi,pi])
xlabel('x'),ylabel('y')
习题2：clear all
thi=0:pi/50:2*pi;
r=1-cos(thi);
polar(thi,r)
title('p=1-cos\theta')
说明：polar currently does not support code generation.
可以在figure窗口中点中曲线在单击鼠标右键选择line width 进行线宽设置。目前polar
指令不支持线宽设置，它采用的是默认线宽。
习题5：clear all
t=0:0.01:13;
x=sin(t);y=cos(t);z=t;
plot3(x,y,z,'g-')
view([-33,77]),box on
上机6：1.课后习题1
for循环语句。
sum=0;
for n=0:1000000
sum=sum+(0.2)^n;
endk=sum;kk =
while循环语句。
n=0;sum=0;
while n<=1000000
sum=sum+(0.2)^n;
n=n+1;
endk=sum;kk =
数值计算程序：
clear all
n=0:1000000;
k=sum((0.2).^n)k =
符号计算程序：
clear all
syms n
k=symsum((0.2)^n,n,0,1000000);
kk=vpa(k)
k=single(kk)kk =
k =function cirline(n,r,str)
cirline(n) without input,plot the unit circle; when input is number n which is greater than 2,plot regular polygon with mumbers of n sides.
cirline(n)
n the number of sides
cirline 画单位圆。
cirline(n) n为大于2的数时，画正n边形。
tian tian编写于2008-11-18
switch nargin
case 0
n=100;r=1;str='-b';
case 1
r=1;str='-b';
case 2
str='-b';
case 3
otherwise
error('输入量太多。')
end;t=0:2*pi/n:2*pi;
x=r*sin(t);y=r*cos(t);
if nargin==0
plot(x,y,str,'linewidth',3);
title('circle')
axis off
elseif n>2
plot(x,y,str,'linewidth',3);
title('poly gon nit6 edges')
axis off
elseerror('出错提示')
endcirline(6)
matlab上机实验答案
三 假设已知矩阵，试给出相应的matlab命令，将其全部偶数行提取出来，赋给矩阵，用命令生成矩阵，用上述命令检验一下结果是不是正确。a magic 8 b a 2 2 end,五 选择合适的步距绘制出下面的图形。1 其中 2 其中。1 t 1 0.0001 1 y sin 1.t plot t,y ...
Matlab上机实习
matlab上机实习报告。一 实验内容。实验四。1 使用函数，实现方阵左旋90 或右旋90 的功能。例如，原矩阵为a,a左旋后得到b,右旋后得到c。b rot90 a b 10 11 12 c rot90 a,3 c 3 2 1 2 建立一个方阵a，求a的逆矩阵和a的行列式的值，并验证a与a 1是互...
MATLAB上机作业
高等代数机算与应用作业题。一 机算题。1 利用函数rand和函数round构造一个5 5的随机正整数矩阵a和b。解 在command window中输入如下内容 a round rand 5 10 结果如下 a 再键入 b round rand 5 10 结果如下 b 1 计算a b，a b和6a ...