推 enoz:恩~~~~try是一種胺基酸= =" 01/26 20:36
我有三個變數 n1 n2 n3 (RANGE都在0~30之間)
帶入23條方程式可以得到23個值 d1 d2 d3.......d23
再把它分類成BC兩組
B=[d1 d2 d3 d4 d5 d6 ]
C=[d7 d8 ............d23]
D=max(B)
E=min(C)
要Try 出一組n1 n2 n3 使 E>D
我產生200000組亂數
都TRY不來
我是這樣寫的
=====================================================
clear
clc
N1=[];N2=[];N3=[];
for n=1:200000
x=0+30*rand(200000,3);
n1=x(n,1);
n2=x(n,2);
n3=x(n,3);
a1=(4*(15.3-n1)^2+(0-n2)^2+(0-n3)^2)^0.5 ;
a2=(4*(17.8-n1)^2+(1-n2)^2+(3.1-n3)^2)^0.5;
a3=(4*(17.8-n1)^2+(0-n2)^2+(2.7-n3)^2)^0.5;
a4=(4*(15.8-n1)^2+(5.3-n2)^2+(7.2-n3)^2)^0.5;
a5=(4*(18-n1)^2+(1.4-n2)^2+(2-n3)^2)^0.5;
a6=(4*(15.2-n1)^2+(5.6-n2)^2+(9.1-n3)^2)^0.5;
a7=(4*(18.4-n1)^2+(0-n2)^2+(2-n3)^2)^0.5;
a8=(4*(16-n1)^2+(9-n2)^2+(5.1-n3)^2)^0.5;
a9=(4*(17.8-n1)^2+(3.1-n2)^2+(5.7-n3)^2)^0.5;
a10=(4*(16.8-n1)^2+(5.7-n2)^2+(8-n3)^2)^0.5;
a11=(4*(16.6-n1)^2+(7.9-n2)^2+(6.5-n3)^2)^0.5;
a12=(4*(15.5-n1)^2+(10.4-n2)^2+(7-n3)^2)^0.5;
a13=(4*(19-n1)^2+(1.8-n2)^2+(7.4-n3)^2)^0.5;
a14=(4*(20-n1)^2+(8.6-n2)^2+(4.1-n3)^2)^0.5;
a15=(4*(16-n1)^2+(5.7-n2)^2+(15.8-n3)^2)^0.5;
a16=(4*(15.8-n1)^2+(6.1-n2)^2+(16.4-n3)^2)^0.5;
a17=(4*(18.4-n1)^2+(6.3-n2)^2+(13.7-n3)^2)^0.5;
a18=(4*(15.3-n1)^2+(18-n2)^2+(6.1-n3)^2)^0.5;
a19=(4*(17.4-n1)^2+(13.7-n2)^2+(11.3-n3)^2)^0.5;
a20=(4*(15.8-n1)^2+(8.8-n2)^2+(19.4-n3)^2)^0.5;
a21=(4*(18.4-n1)^2+(16.4-n2)^2+(10.2-n3)^2)^0.5;
a22=(4*(15.1-n1)^2+(12.3-n2)^2+(22.3-n3)^2)^0.5;
a23=(4*(15.6-n1)^2+(16-n2)^2+(42.3-n3)^2)^0.5;
B=[a23,a10,a15,a16,a17,a19,a20,a21,a22]; %bad solvent
C=[a1,a2,a3,a4,a5,a6,a7,a8,a9,a11,a12,a13,a14,a18]; % good solvet
D=max(C);
E=min(B);
if E>D
n1=n1
n2=n2
n3=n3
N1=[N1 n1];
N2=[N2 n2];
N3=[N3 n3];
disp('one ans')
disp('=======================================')
end
end
N=[N1' N2' N3'];
===============================================================
可是跑了半天還是沒有答案
有沒有辦法把 0~30中所有 n1 n2 n3 的組合產生來TRY
有辦法嗎??
謝謝了
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.115.55.210