MATLAB三 1

发布 2022-09-20 20:17:28 阅读 8973

一。用克莱姆(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...