fangchaneijit20.m

来自「数学建模竞赛：DNA序列的分析」· M 代码 · 共 108 行

a=[1.92054829335959
1.87420481222328
2.99877425939725
1.99108685834404
2.58714772159097
2.22265998103721
2.77339135441486
2.12248112085229
4.98789744071797
2.19089023002066
1.68477343831791
5.30876695061553
1.79078098932397
2.56996947122022
2.64540892978706
2.91717010524325
2.61033576649687
2.43241991988774
2.32879377001924
2.75122955590116
]
t=[2.26172160430602
3.43019442105455
5.87672933802624
2.26803028036598
3.67069829515121
2.46256580486697
4.21272305915337
1.60391677452661
5.46079311436613
2.38198365944147
2.45526782265635
2.40963673742521
1.96396101212393
9.14239210126832
3.48341937273074
4.07839837844874
4.97863857949773
1.95134699683086
3.28120470983236
2.22673153598123
]
c=[2.70801280154532
2.67435485922184
4.09756024970958
3.8328157000404
3.35395465945912
2.34731366095903
4.0275367945655
4.86441157798145
4.90238037419918
2.60320315508926
3.54668117659609
3.66238221356822
3.37950479087971
3.95028392435328
4.64347981172577
2.71277344887962
3.16227766016838
4.38431547932197
2.14269293295285
5.05210941328566
]
g=[3.71756122181041
1.84842275106824
2.1901109697328
2.78601951929716
2.40787239328948
2.78860463125802
2.46779614743843
4.0629241859991
1.95442253411601
1.49240501448927
3.09159683893174
3.08327970781338
2.30799476317426
1.98306320689475
1.92182325789369
1.90424626720229
2.00567377026456
4.29946277073765
2.4181850730701
4.52543171020652]


s=[a,t,c,g]
as=[3.21472881345033 7.73920234155702 5.353333556465 2.55596503346531]
bs=[3.52246934793863 1.85980014068493 12.5705612573442 10.8411699404179]
sas=s*as'

ss=sum((s.*s)')
ass=sum(as.*as)
bss=sum(bs.*bs)
A=sqrt((ss*ass))

neijiA=sas./A'

sbs=s*bs'
B=sqrt(ss*bss)
neijiB=sbs./B'
biji=[neijiA,neijiB]

neijiA-neijiB