-
Notifications
You must be signed in to change notification settings - Fork 1
/
h_q.m
32 lines (29 loc) · 3.04 KB
/
h_q.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
function [y_q] = h_q(Y,alpha,thm,thp)
y_q(1,1)=Y(2) - alpha(1,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^6 - (alpha(1,7)*(Y(1) +...
Y(2) + Y(4)/2 - thp)^6)/(thm - thp)^6 + (6*alpha(1,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5*...
(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp) + (6*alpha(1,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*...
(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^5 - (15*alpha(1,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) -...
1)^4*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm - thp)^2 + (20*alpha(1,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm -...
thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^3)/(thm - thp)^3 - (15*alpha(1,5)*((Y(1) + Y(2) + Y(4)/2 -...
thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^4;
y_q(2,1)=Y(3) - alpha(2,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^6 - (alpha(2,7)*(Y(1) +...
Y(2) + Y(4)/2 - thp)^6)/(thm - thp)^6 + (6*alpha(2,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5*...
(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp) + (6*alpha(2,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*...
(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^5 - (15*alpha(2,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) -...
1)^4*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm - thp)^2 + (20*alpha(2,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm -...
thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^3)/(thm - thp)^3 - (15*alpha(2,5)*((Y(1) + Y(2) + Y(4)/2 -...
thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^4;
y_q(3,1)=Y(4) - alpha(3,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^6 - (alpha(3,7)*(Y(1) +...
Y(2) + Y(4)/2 - thp)^6)/(thm - thp)^6 + (6*alpha(3,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5*...
(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp) + (6*alpha(3,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*...
(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^5 - (15*alpha(3,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) -...
1)^4*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm - thp)^2 + (20*alpha(3,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm -...
thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^3)/(thm - thp)^3 - (15*alpha(3,5)*((Y(1) + Y(2) + Y(4)/2 -...
thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^4;
y_q(4,1)=Y(5) - alpha(4,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^6 - (alpha(4,7)*(Y(1) +...
Y(2) + Y(4)/2 - thp)^6)/(thm - thp)^6 + (6*alpha(4,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5*...
(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp) + (6*alpha(4,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*...
(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^5 - (15*alpha(4,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) -...
1)^4*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm - thp)^2 + (20*alpha(4,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm -...
thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^3)/(thm - thp)^3 - (15*alpha(4,5)*((Y(1) + Y(2) + Y(4)/2 -...
thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^4;