generate.m

高速公路路面质量改进方案 采用matlab编程 达到了稳定点

function M=generate(n)
syms p1 p2 
A=diag([sqrt((1-p2)^2-p1^2)*(1+p2)/(1-p2);sqrt(1-p2^2);ones(n-2,1)]);
B=diag([-p1*sqrt((1+p2)/(1-p2));-p1*ones(n-2,1)],-1);
C=diag([-p2*ones(n-2,1)],-2);
M=A+B+C 
M = inv(M)