练习:一):先求下列表达式的值,然后显示在matlab工作空间的使用情况并全部保存到变量名为: xiti
z1=sind(85)/(1+exp(2))=0.1187
x=(5+cosd(47))/1+sqrt(7)-2i)=1.1980 + 0.6572i
3)已知:x=0:0.5:8 ,求y=sinx2
x=0:0.5:8
y=sin(x.^2)
plot(x,y)
二)保存到mb_01, mb_02 ……
1.写出完成下列操作的命令。
1)建立3阶单位矩阵a; a=eye(3)
2)建立5x6阶随机矩阵a1; a1=rand(5,6)
3)将矩阵a1第2~5行中的第1,3,5列元素赋给矩阵b; b=[a1(2:5,1),a1(2:5,3),a1(2:5,5)]
4)删除矩阵a的第7号元素;a(7)=[
三).设:1) 求:f(x)+g(x) ,f(x)-g(x);
2) f(x)*g(x) ,f(x)/g(x)
syms x
f1=3*x^5-5*x^4+2*x^3-7*x^2+5*x+6
f2=3*x^2+5*x-3
1) f1+f2= 3*x^5-5*x^4+2*x^3-4*x^2+10*x+3
f1-f2= 3*x^5-5*x^4+2*x^3-10*x^2+9
2)f1*f2= (3*x^5-5*x^4+2*x^3-7*x^2+5*x+6)*(3*x^2+5*x-3)
f1/f2= (3*x^5-5*x^4+2*x^3-7*x^2+5*x+6)/(3*x^2+5*x-3)
四)求多项式x4+8x3-10的根。
clear all;clc;
p=[1 8 0 0 10]
poly2sym(p)
x=roots(p)
m=poly(x)p =
ans =x^4+8*x^3+10x =
0.5569 + 0.8918i
0.5569 - 0.8918im =
3. 计算与的数组乘积。
a=[6,9,3;2,7,5]
b=[2,4,1;4,6,8]
a.*bans =
3. 对于,如果,,求解x。
a=[4,9,2;7,6,4;3,5,7]
b=[37;26;28]
x=a\bx =
4, 求解多项式x3-7x2+2x+40的根。
clear all;clc;
p=[1 -7 2 40]
poly2sym(p)
x=roots(p)
m=poly(x)
ans =x^3-7*x^2+2*x+40x =
m =5.求下列极限。
syms a x
f=1+a/x
limit(f,x,inf,'left') ans =
syms a x m
f=(x^(1/m)-a^(1/m))/x-a)
limit(f,x,a)
ans =a^(1/m)/a/m
6.已知:x=0:0.5:8 ,求y=sinx2,并绘制图形图形,要求1)采用红色虚线。
2)作坐标轴,图形标注,加图例。
3)加网格线,4)对点(0,0.7)处添加说明‘曲线’,1)
x=0:0.5:8
y=sin (x.*x)
plot (x,y,'r--'
xlabel(‘x=0:0.5:8’)
ylabel('y=sin (x.*x)')
legend('y=sin (x.*x)')
grid on
text(0,0.7,'曲线')
7.求积分和。
syms x
f=sqrt(x^3+x^4)
int(f,x)
ans =
1/48*(x^3+x^4)^(1/2)*(16*(x^2+x)^(3/2)+12*(x^2+x)^(1/2)*x+6*(x^2+x)^(1/2)-3*log(1/2+x+(x^2+x)^(1/2)))x/((x+1)*x)^(1/2)
syms x
f=x^3/ (x-1)^100
int(f,x,2,3)ans =
8.已知:f(x)=3x5+4x3-5x2-7.2x+5
(1)计算f(x)=0的全部根。
(2)计算f(5), f(7.8), f(9.6), f(12.3)的值。
1)p=[3 0 4 -5 -7.2 5]
poly2sym(p)
x=roots(p)
x =-0.3046 + 1.6217i
-0.3046 - 1.6217i
2) p=[3 0 4 -5 -7.2 5]
poly2sym(p,5)
poly2sym(p,7.8)
poly2sym(p,9.6)
poly2sym(p,12.3)ans =
ans =ans =
ans =
9.已知:f(x)=3x5-5x4+2x3-7x2-5x+6
g(x)=3x2+5x-3
(1) 求f(x)+g(x)和f(x)-g(x)。
p=[3 -5 2 -7 -5 6]
q=[3 5 -3]
f=poly2sym(p)
g=poly2sym(q)
f+gf-g
ans =3*x^5-5*x^4+2*x^3-4*x^2+3
ans =3*x^5-5*x^4+2*x^3-10*x^2-10*x+9
2)求f(x) g(x)和f(x)/g(x)。
simplif(f*g)
simplify(f/g)
ans =9*x^7-28*x^5+4*x^4-56*x^3+14*x^2+45*x-18
ans =3*x^5-5*x^4+2*x^3-7*x^2-5*x+6)/(3*x^2+5*x-3)
10. 在一个图形窗口中以子图形式同时绘制正弦、余弦、正切、余切曲线。
x=linspace(0,6)
subplot(221)
plot(x,sin(x))
xlabel('x')
ylabel('sin(x)')
title('sin(x)')
subplot(222)
plot(x,cos(x))
xlabel('x')
ylabel('cos(x)')
title('cos(x)')
y1=tan(x);y2=cot(x)
subplot(223)
plot(x,y1)
xlabel('x')
title('tan(x)')
subplot(224)
plot(x,y2)
xlabel('x')
title('cot(x)')
ylabel('cot(x)')
ylabel('tan(x)')
11.化简:s=(x2+y2)2+(x2-y2)2
syms x y
s=(x^2+y^2)^2+(x^2-y^2)^2
simplify(s)
ans =
2*x^4+2*y^4
12.求级数之和:
syms n x
s=symsum(n*x^n,n,x,inf)
s =1/(x-1)^2*(-x^(2+x)+2*x*x^x)
13.已知表达式为试编写。
能对任意一组输入x值求相应y值的程序,令函数名为fun13。
function y=fun13(x)
x=input('x=?'
if x>=1
y=10;elseif x>1&x<1
y=0;else
y=-10end
14.编写一个函数文件,用于求两个矩阵的乘积和点乘,然后在命令文件中调用该函数。假设给定矩阵:
1)a为一个三阶单位矩阵,a1为2行3列随机矩阵。
2)b1为由第一小题的第(3)问求得的结果矩阵b而生成的新矩阵,即抽取b的1~3行所有列元素得到b1。
试计算:a*b1与a.*b1,a1*b1,a1.*b1(可以吗?)并比较结果,体会乘与点乘的区别。
建立m文件。
1)function z=chengji(x,y)
z=x*y2)function z=diancheng(x,y)
z=x.*y
a=eye(3)
a1=rand(2,3)
a1=rand(5,6)
b=[a1(2:5,1),a1(2:5,3),a1(2:5,5)]
b1=[b(1:3,1),b(1:3,2),b(1:3,3)]
q1]=chengji(a,b1)
q2]=diancheng(a,b1)
q3]=chengji(a1,b1)
q4]=diancheng(a1,b1)
15以子图的形式绘制曲线:
x=linspace(0,6)
t=linspace(0,6)
subplot(221)
y1=100./(1+x.^2)
plot(x,y1)
xlabel('x')
ylabel('y')
title('y=100/(1+x^2)')
subplot(222)
y2=sqrt(1-x.^2)
plot(x,y2)
xlabel('x')
ylabel('y')
title('x^2+y^2=1')
subplot(223)
y3=1/(2*pi)*exp(-x.^2/2)
plot(x,y3)
xlabel('x')
ylabel('y')
title(' y=1/(2*pi)*exp(-x^2/2)')
subplot(224)
MATLAB习题
实验一 matlab运算基础。班级 01 学号 201105810118 姓名 王宝龙 一 实验目的 1.熟悉启动和退出matlab的方法。2.熟悉matlab命令窗口的组成。3.掌握建立矩阵的方法。4.掌握matlab各种表达式的书写规则及常用函数的使用。二 实验内容 1.求下列表达式的值,观察m...
Matlab1习题
第一章matlab 入门1 习题11.执行下列指令,观察其运算结果,理解其意义 1 12 34 10 2i 5 exp 12 34 6 log 110100 7 prod 12 34 8 a,b min 1020 3040 9 abs 12 34 pi 11 find 1020 3040 40,30...
matlab练习题
向量和矩阵。1.向量a b 计算向量a,b的数量积,向量积和对应元素乘积。解 a 1,2,3 a b 2,1,5 b dot a,b ans cross a,b ans a.bans 2.用matlab生成下列矩阵 解 a 3 4 9 1 0 7 2 1 4 a 解 b zeros 2,2 b 解 ...