standardtmhmmcreator.java

马尔科夫模型的c语言实现

	modelmaker.setTransition("in18", "e");
	tot = setChiSquareDistrib("in18", 17, 18, f_i, q_i, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("in18");
	
	modelmaker.setTransition("in19", "in20");
	modelmaker.setTransition("in19", "in21");
	modelmaker.setTransition("in19", "inglob1");
	modelmaker.setTransition("in19", "e");
	tot = setChiSquareDistrib("in19", 19, 20, f_i, q_i, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("in19");
	
	modelmaker.setTransition("in20", "in21");
	modelmaker.setTransition("in20", "e");
	modelmaker.initializeTransitionProbabilities("in20");

	modelmaker.setTransition("in21", "in22");
	modelmaker.setTransition("in21", "e");
	modelmaker.initializeTransitionProbabilities("in21");
	
	modelmaker.setTransition("in22", "in23");
	modelmaker.setTransition("in22", "e");
	modelmaker.initializeTransitionProbabilities("in22");

	modelmaker.setTransition("in23", "in24");
	modelmaker.setTransition("in23", "e");
	modelmaker.initializeTransitionProbabilities("in23");
	
	modelmaker.setTransition("in24", "in25");
	modelmaker.setTransition("in24", "e");
	modelmaker.initializeTransitionProbabilities("in24");

	modelmaker.setTransition("in25", "in26");
	modelmaker.setTransition("in25", "e");
	modelmaker.initializeTransitionProbabilities("in25");

	modelmaker.setTransition("in26", "in27");
	modelmaker.setTransition("in26", "e");
	modelmaker.initializeTransitionProbabilities("in26");

	modelmaker.setTransition("in27", "in28");
	modelmaker.setTransition("in27", "e");
	modelmaker.initializeTransitionProbabilities("in27");

	modelmaker.setTransition("in28", "in29");
	modelmaker.setTransition("in28", "e");
	modelmaker.initializeTransitionProbabilities("in28");

	modelmaker.setTransition("in29", "ohelixi1");
	modelmaker.setTransition("in29", "e");
	modelmaker.initializeTransitionProbabilities("in29");

	modelmaker.setTransition("inglob1", "in20");
	modelmaker.setTransition("inglob1", "e");
	Module m1 = modelmaker.getModule("inglob1");
	LinkedList l = m1.getVertices();
	Vertex v = ((Vertex)l.get(0));
	v.setTransitionProbability(v.getNumber(), ((double)(lengthI - 1))/(double)(lengthI));
	modelmaker.initializeTransitionProbabilities("inglob1");
	
	modelmaker.setTransition("outglob10", "outglob11");
	//modelmaker.setTransition("outglob10", "outglob29");
	//modelmaker.setTransition("outglob10", "ihelixo1");
	modelmaker.setTransition("outglob10", "e");
	modelmaker.initializeTransitionProbabilities("outglob10");

	modelmaker.setTransition("outglob11", "outglob12");
	//modelmaker.setTransition("outglob11", "outglob28");
	//modelmaker.setTransition("outglob11", "outglob29");
	modelmaker.setTransition("outglob11", "e");
	modelmaker.initializeTransitionProbabilities("outglob11");

	modelmaker.setTransition("outglob12", "outglob13");
	//modelmaker.setTransition("outglob12", "outglob27");
	//modelmaker.setTransition("outglob12", "outglob28");
	modelmaker.setTransition("outglob12", "e");
	modelmaker.initializeTransitionProbabilities("outglob12");

	modelmaker.setTransition("outglob13", "outglob14");
	//modelmaker.setTransition("outglob13", "outglob26");
	//modelmaker.setTransition("outglob13", "outglob27");
	modelmaker.setTransition("outglob13", "e");
	modelmaker.initializeTransitionProbabilities("outglob13");

	modelmaker.setTransition("outglob14", "outglob15");
	//modelmaker.setTransition("outglob14", "outglob25");
	//modelmaker.setTransition("outglob14", "outglob26");
	modelmaker.setTransition("outglob14", "e");
	modelmaker.initializeTransitionProbabilities("outglob14");	
	
	modelmaker.setTransition("outglob15", "outglob16");
	//modelmaker.setTransition("outglob15", "outglob24");
	//modelmaker.setTransition("outglob15", "outglob25");
	modelmaker.setTransition("outglob15", "e");
	modelmaker.initializeTransitionProbabilities("outglob15");
	
	modelmaker.setTransition("outglob16", "outglob17");
	//modelmaker.setTransition("outglob16", "outglob23");
	//modelmaker.setTransition("outglob16", "outglob24");
	modelmaker.setTransition("outglob16", "e");
	modelmaker.initializeTransitionProbabilities("outglob16");

	modelmaker.setTransition("outglob17", "outglob18");
	//modelmaker.setTransition("outglob17", "outglob22");
	//modelmaker.setTransition("outglob17", "outglob23");
	modelmaker.setTransition("outglob17", "e");
	modelmaker.initializeTransitionProbabilities("outglob17");

	modelmaker.setTransition("outglob18", "outglob19");
	//modelmaker.setTransition("outglob18", "outglob21");
	//modelmaker.setTransition("outglob18", "outglob22");
	modelmaker.setTransition("outglob18", "e");
	modelmaker.initializeTransitionProbabilities("outglob18");
	
	//modelmaker.setTransition("outglob19", "outglob20");
	//modelmaker.setTransition("outglob19", "outglob21");
	modelmaker.setTransition("outglob19", "outglobLong");
	modelmaker.setTransition("outglob19", "e");
	modelmaker.initializeTransitionProbabilities("outglob19");
	
	modelmaker.setTransition("outglob20", "outglob21");
	modelmaker.setTransition("outglob20", "e");
	modelmaker.initializeTransitionProbabilities("outglob20");

	modelmaker.setTransition("outglob21", "outglob22");
	modelmaker.setTransition("outglob21", "e");
	modelmaker.initializeTransitionProbabilities("outglob21");
	
	modelmaker.setTransition("outglob22", "outglob23");
	modelmaker.setTransition("outglob22", "e");
	modelmaker.initializeTransitionProbabilities("outglob22");

	modelmaker.setTransition("outglob23", "outglob24");
	modelmaker.setTransition("outglob23", "e");
	modelmaker.initializeTransitionProbabilities("outglob23");
	
	modelmaker.setTransition("outglob24", "outglob25");
	modelmaker.setTransition("outglob24", "e");
	modelmaker.initializeTransitionProbabilities("outglob24");

	modelmaker.setTransition("outglob25", "outglob26");
	modelmaker.setTransition("outglob25", "e");
	modelmaker.initializeTransitionProbabilities("outglob25");

	modelmaker.setTransition("outglob26", "outglob27");
	modelmaker.setTransition("outglob26", "e");
	modelmaker.initializeTransitionProbabilities("outglob26");

	modelmaker.setTransition("outglob27", "outglob28");
	modelmaker.setTransition("outglob27", "e");
	modelmaker.initializeTransitionProbabilities("outglob27");

	modelmaker.setTransition("outglob28", "outglob29");
	modelmaker.setTransition("outglob28", "e");
	modelmaker.initializeTransitionProbabilities("outglob28");

	modelmaker.setTransition("outglob29", "ihelixo1");
	modelmaker.setTransition("outglob29", "e");
	modelmaker.initializeTransitionProbabilities("outglob29");
	
	modelmaker.setTransition("outglobLong", "outglob20");
	modelmaker.setTransition("outglobLong", "e");
	m1 = modelmaker.getModule("outglobLong");
	l = m1.getVertices();
	v = ((Vertex)l.get(0));
	v.setTransitionProbability(v.getNumber(), ((double)(lengthOL - 1))/(double)(lengthOL));
	modelmaker.initializeTransitionProbabilities("outglobLong");
	
	distribSum = 0.0;
	for(int i = 1; i <= 20; i++) {
	    distribSum += Math.pow(i, (f_o/2 - 1)) * Math.exp((0-i)/q_o);
	}
	tot = 1.0;
	
	modelmaker.setTransition("out10", "out11");
	modelmaker.setTransition("out10", "out29");
	modelmaker.setTransition("out10", "ihelixo1");
	modelmaker.setTransition("out10", "e");
	tot = setChiSquareDistrib("out10", 1, 2, f_o, q_o, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("out10");

	modelmaker.setTransition("out11", "out12");
	modelmaker.setTransition("out11", "out28");
	modelmaker.setTransition("out11", "out29");
	modelmaker.setTransition("out11", "e");
	tot = setChiSquareDistrib("out11", 3, 4, f_o, q_o, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("out11");

	modelmaker.setTransition("out12", "out13");
	modelmaker.setTransition("out12", "out27");
	modelmaker.setTransition("out12", "out28");
	modelmaker.setTransition("out12", "e");
	tot = setChiSquareDistrib("out12", 5, 6, f_o, q_o, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("out12");

	modelmaker.setTransition("out13", "out14");
	modelmaker.setTransition("out13", "out26");
	modelmaker.setTransition("out13", "out27");
	modelmaker.setTransition("out13", "e");
	tot = setChiSquareDistrib("out13", 7, 8, f_o, q_o, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("out13");

	modelmaker.setTransition("out14", "out15");
	modelmaker.setTransition("out14", "out25");
	modelmaker.setTransition("out14", "out26");
	modelmaker.setTransition("out14", "e");
	tot = setChiSquareDistrib("out14", 9, 10, f_o, q_o, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("out14");	
	
	modelmaker.setTransition("out15", "out16");
	modelmaker.setTransition("out15", "out24");
	modelmaker.setTransition("out15", "out25");
	modelmaker.setTransition("out15", "e");
	tot = setChiSquareDistrib("out15", 11, 12, f_o, q_o, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("out15");
	
	modelmaker.setTransition("out16", "out17");
	modelmaker.setTransition("out16", "out23");
	modelmaker.setTransition("out16", "out24");
	modelmaker.setTransition("out16", "e");
	tot = setChiSquareDistrib("out16", 13, 14, f_o, q_o, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("out16");

	modelmaker.setTransition("out17", "out18");
	modelmaker.setTransition("out17", "out22");
	modelmaker.setTransition("out17", "out23");
	modelmaker.setTransition("out17", "e");
	tot = setChiSquareDistrib("out17", 15, 16, f_o, q_o, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("out17");

	modelmaker.setTransition("out18", "out19");
	modelmaker.setTransition("out18", "out21");
	modelmaker.setTransition("out18", "out22");
	modelmaker.setTransition("out18", "e");
	tot = setChiSquareDistrib("out18", 17, 18, f_o, q_o, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("out18");
	
	modelmaker.setTransition("out19", "out20");
	modelmaker.setTransition("out19", "out21");
	modelmaker.setTransition("out19", "outglobShort");
	modelmaker.setTransition("out19", "e");
	tot = setChiSquareDistrib("out19", 19, 20, f_o, q_o, distribSum, tot);
	modelmaker.initializeTransitionProbabilities("out19");
	
	modelmaker.setTransition("out20", "out21");
	modelmaker.setTransition("out20", "e");
	modelmaker.initializeTransitionProbabilities("out20");

	modelmaker.setTransition("out21", "out22");
	modelmaker.setTransition("out21", "e");
	modelmaker.initializeTransitionProbabilities("out21");
	
	modelmaker.setTransition("out22", "out23");
	modelmaker.setTransition("out22", "e");
	modelmaker.initializeTransitionProbabilities("out22");

	modelmaker.setTransition("out23", "out24");
	modelmaker.setTransition("out23", "e");
	modelmaker.initializeTransitionProbabilities("out23");
	
	modelmaker.setTransition("out24", "out25");
	modelmaker.setTransition("out24", "e");
	modelmaker.initializeTransitionProbabilities("out24");

	modelmaker.setTransition("out25", "out26");
	modelmaker.setTransition("out25", "e");
	modelmaker.initializeTransitionProbabilities("out25");

	modelmaker.setTransition("out26", "out27");
	modelmaker.setTransition("out26", "e");
	modelmaker.initializeTransitionProbabilities("out26");

	modelmaker.setTransition("out27", "out28");
	modelmaker.setTransition("out27", "e");
	modelmaker.initializeTransitionProbabilities("out27");

	modelmaker.setTransition("out28", "out29");
	modelmaker.setTransition("out28", "e");
	modelmaker.initializeTransitionProbabilities("out28");

	modelmaker.setTransition("out29", "ihelixo1");
	modelmaker.setTransition("out29", "e");
	modelmaker.initializeTransitionProbabilities("out29");

	modelmaker.setTransition("outglobShort", "out20");
	modelmaker.setTransition("outglobShort", "e");
	m1 = modelmaker.getModule("outglobShort");
	l = m1.getVertices();
	v = ((Vertex)l.get(0));
	v.setTransitionProbability(v.getNumber(), ((double)(lengthOS - 1))/(double)(lengthOS));
	modelmaker.initializeTransitionProbabilities("outglobShort");

	modelmaker.setTransition("ohelixi1", "ohelixi2");
	modelmaker.initializeTransitionProbabilities("ohelixi1");
	modelmaker.setTransition("ohelixi2", "ohelixi3");
	modelmaker.initializeTransitionProbabilities("ohelixi2");
	modelmaker.setTransition("ohelixi3", "ohelixi4");
	modelmaker.initializeTransitionProbabilities("ohelixi3");
	modelmaker.setTransition("ohelixi4", "ohelixi5");
	modelmaker.initializeTransitionProbabilities("ohelixi4");
	modelmaker.setTransition("ohelixi5", "ohelixi6");
	modelmaker.initializeTransitionProbabilities("ohelixi5");	
	modelmaker.setTransition("ohelixi6", "ohelixi7");
	modelmaker.initializeTransitionProbabilities("ohelixi6");
	modelmaker.setTransition("ohelixi7", "ohelix");
	modelmaker.initializeTransitionProbabilities("ohelixi7");	

	modelmaker.setTransition("ihelixo1", "ihelixo2");
	modelmaker.initializeTransitionProbabilities("ihelixo1");
	modelmaker.setTransition("ihelixo2", "ihelixo3");
	modelmaker.initializeTransitionProbabilities("ihelixo2");
	modelmaker.setTransition("ihelixo3", "ihelixo4");