一。用克莱姆(cramer)方法求解线性方程组。
d=[2,2,-1,1;4,3,-1,2;8,5,-3,4;3,3,-2,2]; 定义系数矩阵。
b=[4;6;12;6]; 定义常数项向量。
d1=[b,d(:,2:4)];用方程组的右端向量置换d的第1列。
d2=[d(:,1:1),b,d(:,3:4)];用方程组的右端向量置换d的第2列。
d3=[d(:,1:2),b,d(:,4:4)];用方程组的右端向量置换d的第3列。
d4=[d(:,1:3),b]; 用方程组的右端向量置换d的第4列。
dd=det(d);
x1=det(d1)/dd;
x2=det(d2)/dd;
x3=det(d3)/dd;
x4=det(d4)/dd;
x1,x2,x3,x4]
二。x(0:pi/100:2*pi)';
y1=2*exp(-0.5*x)*[1,-1];
y2=2*exp(-0.5*x).*sin(2*pi*x);
plot(x,y1,'b:')
axis([0,2*pi,-2,2]);
hold on;
plot(x,y2,'k');
grid on;
box on;
hold off;
三。e=eye(3);
r=rand(3,2);
o=zeros(2,3);
s=disp(2,3);
a=[e,r,o,s];
a^2b=[e,(r+r*s);o,s^2])
四。function [x_star,index,it]=erfen(fun,a,b,ep)
求解非线性方程的二分法,其中fun(x)为需要求根的函数;a,b为初始区间的端点;ep为精度,当(b-a)/2《ep时终止计算,缺省值为1e-5,x_star为迭代成功时,输出方程的根,迭代不成攻时,输出两端点的函数值;index为指标变量,当index=1时,表明。
迭代成功,index=0,表明初始区间不是有根区间;it为迭代次数。
if nargin<4 ep=1e-5;end
fa=feval(fun,a);fb=feval(fun,b);
if fa*fb>0
x_star=[fa,fb];index=0;it=0;
return;
endk=1;
while abs(b-a)/2>=ep
x=(a+b)/2;fx=feval(fun,x);
if fx*fa<0
b=x;fb=fx;
elsea=x;fa=fx;
endk=k+1;
endx_star=(a+b)/2;index=1;it=k;
五。x=0:(2*pi)/100:2*pi;y=cos(0.5+((3*sin(x)).1+x.^2)))
设 y=cos[0.5+((3sinx)/(1+x^2))]把x=0~2π间分为101点,画出以x为横坐标,y为纵坐标的曲线;
plot(x,y)
xlabel('x');
ylabel('y');
title('y=cos(0.5+((3*sin(x)).1+x.^2)))
六。x=-10:1/20:10;y=x; [x,y]=meshgrid(x,y);
%在[-10,10;-10,10]范围内画出函数的三维图形。
z=sin((x.^2+y.^2).^0.5)./x.^2+y.^2).^0.5); surf(x,y,z)
xlabel('x');
ylabel('y');
zlabel('z');
title('z=sin((x.^2+y.^2).^0.5)./x.^2+y.^2).^0.5)')
七。x = 0 1; 0 1];
clusters = 5;
points = 10;
std_dev = 0.05;
p = nngenc(x,clusters,points,std_dev);
plot(p(1,:)p(2,:)r');
title('输入样本向量');
xlabel('p(1)')
ylabel('p(2)')
建立网络。net=newc([0 1;0 1],5,0.1);
得到网络权值,并在图上绘出。
figure;
plot(p(1,:)p(2,:)r');
w=hold on;
plot(w(:,1),w(:,2),'ob');
hold off;
title('输入样本向量及初始权值');
xlabel('p(1)')
ylabel('p(2)')
figure;
plot(p(1,:)p(2,:)r');
hold on;
训练网络。net=init(net);
net=train(net,p);
得到训练后的网络权值,并在图上绘出。
w=plot(w(:,1),w(:,2),'ob');
hold off;
title('输入样本向量及更新后的权值');
xlabel('p(1)')
ylabel('p(2)')
a=0;p = 0.6 ;0.8];
a=sim(net,p)
八。i = round(100*rand);
go = 1;
promotion = 在1到100的数字中猜一个吧,goodluck!!'
x = input(promotion);
while go ==1
if x ==i
input('you won');
go = 0
elseif x > i
input('high');
elseinput('low');
endend
九。for i=1:10
a=89+i;
b=79+i;
c=69+i;
d=59+i;
endc=[d,c];
name=;
mark=;
rank=cell(1,5);
s=struct('name',name,'mark',mark.'rank',rank);
for i=1:5
switch s(i).mark
case 100
s(i).rank='满分';
case a
s(i).rank='优秀 ';
case b
s(i).rank='良好';
case c
s(i).rank='及格';
otherwise
s(i).rank='不及格';
endend
disp(['学生姓名 ',得分 ',等级 ']
for i=1:5
disp([s(i).name,blanks (6),num2str(s(i).mark),blanks (6),s(i).rank]);end十。
= round(100*rand);
go = 1;
promotion = 在1到100的数字中猜一个吧,goodluck!!'
x = input(promotion);
while go ==1
if x ==i
input('you won');
go = 0
elseif x > i
input('high');
elseinput('low');
endend
十一。function y=f(x)
if x<10
y=1;for t=1:x
y=y*t;
endt=0:0.001:x;
plot3(t,sin(t),cos(t));
grid on;
elseif x<100
y=3*x+4;
t=x:0.001:100;
plot3(t,sin(t),cos(t));
grid on;
end十二。
求a的特征值和特征向量。
a=[1,2,2;1,-1,1;4,-12,1];
e=eig(a)
v,d]=eig(a)
v,d]=eig(a,'nobalance')
十三。x=linspace(0,2*pi,101);
y=cos(0.5+((3*sin(x)).1+x.^2)))
plot(x,y)
十四。a=[1 -4 0 3 -2 6];
x=linspace(-2,8,100);
y=polyval(a,x);
plot(x,y,x,zeros(1,100))
十五。特征值的方法解方程。
p=[3,-7,0,5,2,-18];
a=compan(p); a的友矩阵。
x1=eig(a) %求a的特征值。
x2=roots(p) %直接多项式p的零点。
十六。设银行年利率为11.25%。将10000元钱存入银行,多长时间会连本带利翻一番?
money=100000
years=0
while money<20000
years=years+1
money=money*(1+11025/100)
end十七。
function f=funl(x)
if x>1
f=x^2+1
endif x<=1
f=2*xend
十八。function dicetest(result)
检测掷一次骰子的点数是单数还是双数:
switch result
case disp('odd number of eyes')
case disp('even number of eyes')
otherwise
sidp('what kind of dice do you h**e?')
a = for k = 1:5
for j = 1:5
if k ==j
matlab 1
6.4.3 对象句柄的获取方法。图像对象建立后,句柄值时图形对象的唯一标志。句柄值的获取除可以创建时直接赋值,还有如下几种获取方法 1.追溯法获取句柄。h pa get h known,parent 功能 获取h known的父对象的句柄值,并赋给h pa。h ch get h known,chil...
matlab实验 1
1 了解matlab 用户界面。2 掌握matlab 的路径搜索 3 了解matlab 的帮助系统。4 掌握matlab 中标点符号的使用方法。5 掌握matlab 的数据类型。6 掌握常量和变量的使用方法。7 了解浮点数和复数的使用方法。8 掌握matlab 的文件操作。1 熟悉matlab 7用...
matlab笔记 1
中setdiff 函数作用 判断2个数组中不同元素。c setdiff a,b 返回在a中有,而b中没有的值,结果向量将以升序排序返回。在集合论中,c a b。a和b也可以是字符串细胞数组。2.matlab之size 返回矩阵的行数和列数。a,b size x 返回矩阵的行数a,列数b size x...