。

#xff08;采用单选数#xff0c;以及计算机科学技术中的递推回溯法#xff09;

以上的数独是图14-2的两个矩阵的和#xff0c;左侧的矩阵可以由kron和magic函…目录

四独游戏

采用蛮力试凑法来解决数独问题。

采用单选数以及计算机科学技术中的递推回溯法

以上的数独是图14-2的两个矩阵的和左侧的矩阵可以由kron和magic函数建立起来前一个函数用来求Kronecker乘积后者生成幻方矩阵。

kron(eye(3),

使用4×4网格可能的备选项用小数标注出来。

将单选数填在网格内如果没有单选数则采用递归回溯法来求解。

左有单选数

右无单选数

数独网格的一些运算可能改变其在图形用户界面中的显示但不能改变其基本特性。

所有变化基本上都是相同的数字谜。

这些等效的运算可以用matlab的数组运算表示。

randperm(9),

------------------------------function

[C,N]

0;z(nonzeros(X(tri(i),tri(j))))

0;C{i,j}

0;gui_init(X);endsudoku_gui(X,steps);%

Fill

candidates(X);endsudoku_gui(X,steps,C);%

Recursive

solution.breakendsudoku_gui(X,steps,C,-s);

Revert

------------------------------function

[C,N]

0;z(nonzeros(X(tri(i),tri(j))))

0;C{i,j}

------------------------------function

gui_init(X)%

get(0,defaulttextfontsize);fname

Lucida

Typewriter;clfshgset(gcf,color,white)axis

squareaxis

10]/11,color,dkmagenta,linewidth,4)line([1

10]/11,[m

m]/11,color,dkmagenta,linewidth,4)endH.a

zeros(9,9);for

text((j1/2)/11,(10.5-i)/11,string,

...units,normal,fontsize,fsize6,fontweight,bold,

...fontname,fname,color,color,horizont,center);endendstrings

zeros(1,4);for

uicontrol(style,toggle,string,strings{k},

...units,normal,position,[(k3)*0.125,0.05,0.10,0.05],

...callback,

set(gcf,user,H));endset(H.b(1),style,pushbutton)H.s

1;H.t

title(0,fontweight,bold);set(gcf,userdata,H)drawnowend

gui_init%

------------------------------function

get(gcf,userdata);if

0set(H.t,string,int2str(steps))drawnowendreturnelseset(H.t,string,int2str(steps))endk

[1:H.s-1

H.s1:4];set(H.b(k),value,0);dkblue

2/3];dkred

~isequal(get(H.a(i,j),color),dkblue)

...~isequal(get(H.a(i,j),color),cyan)if

X(i,j)

0set(H.a(i,j),string,int2str(X(i,j)),fontsize,fsize6,

nargin

1set(H.a(i,j),string,char3x3(C{i,j}),fontsize,fsize-4,

...color,dkred)elseset(H.a(i,j),string,char3x3(C{i,j}),fontsize,fsize-4,

...color,dkgreen)endendendendif

nargin

0set(H.a(z),color,cyan)elseset(H.a(-z),color,dkgreen)returnendend%

Gui

0set(H.a(z),color,cyan)elseset(H.a(-z),color,dkgreen)returnendend%

------------------------------function

char3x3(c)%

------------------------------function

gui_finish(X,steps)H

2;set(H.b(1:3),vis,off)set(gcf,userdata,H)set(H.b(4),string,close,value,0,

...callback,close(gcf))sudoku_gui(X,steps)end

gui_finishend

0;z(nonzeros(X(tri(i),tri(j))))

0;C{i,j}

matlab提供了很多求给定常微分方程数组近似解的函数这一常微分方程组数值解的函数包括ode23、ode45、ode113、ode23s、ode15s、ode23t、ode23tb。

函数名中的数字表示所用算法的阶次阶次和算法的复杂程度和精度有关。

所有这些函数都会自动选择近似步长来保证预先选择的精度要求。

所选择的阶次越高每一步计算量越大但是所用的步长也越大。

比如ode23算法比较二阶和三阶算法来估计计算步长而ode45比较的是四阶和五阶算法。

s表示为stiff

ode23(acircle,tspan,y0)plot(y(:,1),y(:,2),-o)axis

squareaxis([-1.1

odeset(outputfcn,odephas2)ode23(acircle,tspan,y0,opts)axis

squareaxis([-1.1