新浪博客

刘卫国 Matlab 例题 1-3章

2011-03-14 12:03阅读：

http://blog.sina.cn/dpool/blog/u/1864621744

第一章
例1.1 分别绘制函数 刘卫国 <wbr>Matlab <wbr>例题 <wbr>1-3章和 刘卫国 <wbr>Matlab <wbr>例题 <wbr>1-3章的曲线。
x=-2*pi:pi/180:2*pi;
plot(x,2.^(-abs(x)),':',x,sin(x));
例1.2 求方程
2x5-3x3 +71x2-9x+13=0的全部根。
p=[2,0,-3,71,-9,13];
x=roots(p)
例1.3 求解线性方程组。
刘卫国 <wbr>Matlab <wbr>例题 <wbr>1-3章
a=[2,3,-1;8,2,3;45,3,9];
b=[2;4;23];
x=inv(a)*b
例1.4 求积分 刘卫国 <wbr>Matlab <wbr>例题 <wbr>1-3章
quad('x.*log(1+x)',0,1)
第二章

例2.1 计算表达式 刘卫国 <wbr>Matlab <wbr>例题 <wbr>1-3章 的值，并将结果赋给变量x，然后显示出结果。
x=(5+cos(47*pi/180))/(1+sqrt(7)-2*i) %计算表达式的值
例2.2 利用M文件建立MYMAT矩阵。
(1)启动有关编辑程序或MATLAB文本编辑器(见第4章)，并输入待建矩阵：
MYMAT=[101,102,103,104,105,106,107,108,109;
201,202,203,204,205,206,207,208,209;
301,302,303,304,305,306,307,308,309];
(2)把输入的内容存盘(设文件名为mymatrix.m)。
(3)在MATLAB命令窗口中输入mymatrix，即运行该M文件，就会自动建立一个名为MYMAT的矩阵，可供以后使用。
例2.3 建立5阶方阵A，判断A的元素是否能被3整除。
A =[24,35,13,22,63;23,39,47,80,80; ...
90,41,80,29,10;45,57,85,62,21;37,19,31,88,76]
P=rem(A,3)==0 %判断A的元素是否可以被3整除
例2.4 在[0,3π]区间，求y=sin(x)的值。要求：
(1)消去负半波，即（π,2π）区间内的函数值置0。
(2) （ 刘卫国 <wbr>Matlab <wbr>例题 <wbr>1-3章, 刘卫国 <wbr>Matlab <wbr>例题 <wbr>1-3章）和（ 刘卫国 <wbr>Matlab <wbr>例题 <wbr>1-3章, 刘卫国 <wbr>Matlab <wbr>例题 <wbr>1-3章）区间内取值均为sin 刘卫国 <wbr>Matlab <wbr>例题 <wbr>1-3章。
方法1：
x=0:pi/100:3*pi;
y=sin(x);
y1=(x<pi|x>2*pi).*y; %消去负半波
q=(x>pi/3&x<2*pi/3)|(x>7*pi/3&x<8*pi/3);
qn=~q;
y2=q*sin(pi/3)+qn.*y1; %按要求处理第(2)步
方法2：
x=0:pi/100:3*pi;
y=sin(x);
y1=(y>=0).*y; %消去负半波
p=sin(pi/3);
y2=(y>=p)*p+(y<p).*y1; %按要求处理第(2)步
例2.5 建立矩阵A，然后找出在[10，20]区间的元素的位置。
(1) 建立矩阵A。
A=[4,15,-45,10,6;56,0,17,-45,0]
A =
4 15 -45 10 6
56 0 17 -45 0
(2) 找出大于4的元素的位置。
find(A>=10 & A<=20)
ans =
3
6

我的更多文章

下载客户端阅读体验更佳

APP专享