高手帮忙~~~帮忙分析一段MATLAB程序~~~~高加分!!!

作者&投稿:蓬荔 (若有异议请与网页底部的电邮联系)
关于Matlab的一个小程序,希望高手帮忙啊。。。~

用perms函数就可以实现了...

你透露的信息太少了,起码应该告知一下这个程序是用来干什么的,大概用的是什么算法啊。
看来只能猜了,从最后画出的图来看,感觉这个程序的目的用逼近的方法得到一条分界线的样子,具体是什么算法,看不出来。你看看我的注解,自己研究一下吧。
在网页里可能显得比较乱,你复制粘贴到matlab的编辑器了再看吧。

w1=[0.1 6.8 -3.5 2.0 4.1 3.1 -0.8 0.9 5.0 3.9; 1.1 7.1 -4.1 2.7 2.8 5.0 -1.3 1.2 6.4 4.0]; %输入第一组十个点的平面坐标w1,对应图中的红点
w2=[7.1 -1.4 4.5 6.3 4.2 1.4 2.4 2.5 8.4 4.1;4.2 -4.3 0.0 1.6 1.9 -3.2 -4.0 -6.1 3.7 -2.2];%输入第二组十个点的平面坐标w2,对应图中的蓝星
ww1=[ones(1,size(w1,2)); w1]; %将2行10列的w1前面增加一个全为1的行,将其扩展为3行10列的ww1
ww2=[ones(1,size(w2,2)); w2]; %将2行10列的w2前面增加一个全为1的行,将其扩展为3行10列的ww2
X=[ww1 -ww2]; %将3行10列的ww1和ww2,合并成3行20列的X
k=1; % k在这里应该是定义了一个逼近次数的计数器
W=ones(size(X,1),1); % 感觉W应该是初始化了一个3行1列的向量,里面存得三个数用于得到最后的那条直线
ok=0; %给ok赋初值0, ok在这里应该是定义了一个循环结束的标志,
while(ok==0) %ok初始值为0,当ok的值变为1的时候,下面的循环语句就跳出了
for i=1:size(X,2) %i从1增加到20,相当于一个数组指针,依次对X中的各个列进行操作
if (W'*X(:,i)<0) %这是if语句整个程序的关键,真正的算法就在这里实现,
k=k+1;
W=W+X(:,i); %如果W与X的第i列的点积小于零,则将W与X的第i列进行加和,并赋值给W
break; %跳出for循环体
else
if (i==size(X,2)) %如果W与X的第i列的点积不小于零,则需要看一下是否已经对X中的20组数都进行了判断,
ok=1; %如果已经对X中的20组数都进行了判断,则结束while循环,结束逼近过程
end
end
end
end
figure(1) %新建一个图像窗口
plot(w1(1,:),w1(2,:),'r.') %绘制第一组的十个点,用红点表示
hold on %保持刚才绘制的十个点不被清除,继续绘图
plot(w2(1,:),w2(2,:),'*') %绘制第二组的十个点,用蓝星表示
xmin=min(min(w1(1,:)),min(w2(1,:))); %得到20个点的最小x值
xmax=max(max(w1(1,:)),max(w2(1,:))); %得到20个点的最大x值
ymin=min(min(w1(2,:)),min(w2(2,:))); %得到20个点的最小y值
ymax=max(max(w1(2,:)),max(w2(2,:))); %得到20个点的最大y值
xindex=xmin-1:(xmax-xmin)/100:xmax+1; %生成一个x序列,用于绘制直线
yindex=-W(2)*xindex/W(3)-W(1)/W(3); %计算直线与上面的x序列相对应的y坐标
plot(xindex,yindex) %绘制直线

常微分(ODE)方程的数值求解器有:非刚性求解器(计算的精度从低到高)ode23,ode45,ode113,刚性方程求解器(适用的刚性从弱到强) ode15s,ode23s,ode23t,ode23tb,隐方程求解器ode15i
所谓刚性方程,就是指它的解的曲线有剧烈的或缓慢的变化。如van der Pol方程(教材p144例4)就是一个刚性方程.
Ode求解器默认的相对误差是1e-3,绝对误差是1e-6,要改变默认的精度设置,可以用odeset来设定Options。具体设置方法,求助于help功能。Ode求解器中可以求解带有参数的微分方程。
常用的精度设置如
Options=odeset(‘RelTol’, 1e-5,’AbsTol’, [1e-8, 1e-7]);
其中绝对误差可以对每个未知函数的分量分别规定,写成一个向量,维数等于方程的维数。如对各分量的绝对误差设置相同则只须写一个标量误差。
每一积分步第i个分量的误差满足e(i) <= max(RelTol*abs(y(i)),AbsTol(i)).
如果只要对解的范数作误差控制,而不需对解的每个分量作误差控制,则在Options中可以加上选项’NormControl’, ‘on’. 这时每一积分步误差的范数满足norm(e) <= max(RelTol*norm(y),AbsTol). 这个选项对那些解的范数会等于零的方程特别有用,不用此选项时,为了要达到苛刻的误差要求,步长会取得很小,将大大减慢求解过程以致求解失败.
时滞是常数的时滞常微分方程DDE的数值求解器有dde23,要改变默认的设置,可以用ddeset来设定。

常微分方程求解器的options还可以设置一个有用的功能,语法是
Options=odeset( 'Events',@EVENTS);(Options可以是自定义名,其中各种设置如精度设置可以写在一起,用逗号分开),用ode45求解时,可用格式
[TOUT,YOUT,TE,YE,IE] = ode45(@ODEFUN,TSPAN,Y0,Options, P)
在odeset中设置了一个事件函数@EVENTS,是自编m函数,函数名自定。函数的格式是
[VALUE,ISTERMINAL,DIRECTION] = EVENTS(T,Y,P).
事件函数的输入是和微分方程的函数输入相同,顺序相同. 事件函数的输出是3个列向量. 例如 [VALUE,ISTERMINAL,DIRECTION],名称自定,列向量的维数是事件的个数,VALUE(I) 是事件函数的第I个分量表达式的值, ISTERMINAL(I)=1表示事件函数的第I个分量的值等于零时积分终止,不然等于0 , DIRECTION(I)=0 表示要计算事件函数的第I个分量的所有的零点 (默认), +1 仅计算事件函数第I个分量在零点递增的零点 -1仅计算事件函数的第I个分量在零点递减的零点.
TSPAN是微分方程求解区间,例如可用 [0 pi]表示,也可以指定TOUT为求解区间中的一些点列,如TSPAN=0: 0.1: 3.2; Y0是微分方程初值问题中的初值,如Y0=[0 0];
TOUT是微分方程求解区间TSPAN中的时间点列,YOUT是时间点列TOUT上方程解的值,是矩阵形式,列数等于方程的维数,TOUT(:,J)是解的第J个分量在时间点列TOUT上的向量值。 TE是列向量,是事件发生的时刻序列. YE的各行向量是事件发生各时刻方程的解的向量值, IE表示在TE时刻发生的事件在事件函数中的序号.

求高手帮忙分析一下这句话,谢谢
答:第一句结构很简单 主语:We 谓语动词:can call 宾语:these people 宾语补足语:"generalists"第二句是一个较长的复杂句 1、从大结构上分析,它是一个多重复合句:第一重 主句:these "generalists" are particularly needed for positions in administration,非限制性定语从句:where it is their job...

请高手帮忙分析下,此段代码的功能是什么??
答:后台设置的一些内容,如类目是不是要显示,类目文本颜色,或者其它等内容,都会存入数据库。而上面这个函数就是从数据库读出数据,根据管理员设置的数据,生成符合管理员要求的菜单部分的源代码。此部分代码分三部分,菜单 + 订阅 + 广告。

请高手帮忙翻译一下下面的一段文字!万分感谢了。
答:security management is an important part in the enterprise management.there are many specifical functional sub-system among the enterprise manament system.Security management is one of them. This sub-system is made up of corresponding person in the corresponding enterprise persons. the Sub-...

高手帮忙分析几个英语句子
答:是个简单句 a)主语:One new trend 一种新的趋势 b)谓语动词:involve 包括 c)宾语:动名词短语 sourcing portions of a research program to another country, 一个对另一个国家的研究计划的一些段落中所获取得来源。丢掉还有写错必要的关键词造成了分析的困难,不过都已经一一改正了。

请高手帮忙分析一下,为什么无线连接电脑只有20m宽带,连接手机却有40m网...
答:第一,那个不是宽带,是带宽。先搞清楚概念 第二,这个20M 40M的,是指定手机到路由器的连接速率,而不是你家宽带的速率 第三,这个改不了,它是路由器与手机/电脑的无线网卡自动匹配的。

麻烦高手帮忙分析一下,谢谢.
答:麻烦高手帮忙分析一下,谢谢.  我来答 1个回答 #热议# 婚姻并不幸福的父母,为什么也会催婚?...基本运势在这段时间的好运主要是和「居家」以及「手足」较有关系。只要能维持居家和手足情谊的稳定,

请高手帮忙翻译一下一段心得体会!
答:对与公司产品有关的资料、说明书、广告等,均必须努力研讨、熟记,同时要收集竞争对手的广告、宣传资料、说明书等,加以研究、分析,以便做到“知己知彼”,只有如此才能在金融危机中采取相应对策生存下来。同时要想作为一名成功的业务员,我们要做前期的市场调查,根据不同年龄段客户的特点采取不同的营销...

恋爱高手帮忙分析一下.
答:恋爱高手帮忙分析一下. 200 我是警察,她是机关的会计.我长的一般,她很漂亮.我家庭一般,她家境挺好的,我情况她也都知道了.我05年大学毕业考的公务员,她06年刚刚毕业的.我比她大几个月.我以前见过不少人家介绍的,其... 我是警察,她是机关的会计.我长的一般,她很漂亮.我家庭一般,她家境挺好的,我情况她...

求高手帮忙翻译一段文字
答:Chapter 1 analysis curved face work piece super voice examination system and it be related technique of research present condition, give thesis of main research contents and it importance meaning.Chapter 2 set up a mold for mathematics for measure of the curved face super voice.First ...

网速慢,请高手帮忙分析下!
答:只有这样,才能最大限度地避免串扰,保证数据传输。本人在实践中发现不按正确标准(T586A、T586B)制作的网线,存在很大的隐患。表现为:一种情况是刚开始使用时网速就很慢;另一种情况则是开始网速正常,但过了一段时间后,网速变慢。后一种情况在台式电脑上表现非常明显,但用笔记本电脑检查时网速却...