17.04.2014, 13:06
Hallo
Ich versuche mit Labwindows CVI von meiner Funktion xd(t) das Integral zu ziehen, bekomm aber mit NumericIntegral leider keine richtige Lösung raus.
[/
double kpr = 1;
double tn = 1;
double kps = 1;
double t1 = 1;
double xd[10000];
double a = 1/(2*t1);
double b = sqrt(1/(t1*tn));
double wd = sqrt(pow(b,2)-pow(a,2));
double fak = ((kpr*kps)/(t1*tn));
double ISE;
for(int t = 0; t< 10000; t++)
{
xd[t] = pow(1-(fak*(1/pow(b,2))*(1-((cos(wd*t)+(a/wd)*sin(wd*t)))*exp(-a*t))),2);
}
NumericIntegration (xd, 10000, 0.1, TRAPEZOIDAL_RULE, &ISE);
]
Mathematisch soll das Integral so aussehen wie im Anhang zu sehen ist.
LG und danke für die Hilfe
Ich versuche mit Labwindows CVI von meiner Funktion xd(t) das Integral zu ziehen, bekomm aber mit NumericIntegral leider keine richtige Lösung raus.
[/
double kpr = 1;
double tn = 1;
double kps = 1;
double t1 = 1;
double xd[10000];
double a = 1/(2*t1);
double b = sqrt(1/(t1*tn));
double wd = sqrt(pow(b,2)-pow(a,2));
double fak = ((kpr*kps)/(t1*tn));
double ISE;
for(int t = 0; t< 10000; t++)
{
xd[t] = pow(1-(fak*(1/pow(b,2))*(1-((cos(wd*t)+(a/wd)*sin(wd*t)))*exp(-a*t))),2);
}
NumericIntegration (xd, 10000, 0.1, TRAPEZOIDAL_RULE, &ISE);
]
Mathematisch soll das Integral so aussehen wie im Anhang zu sehen ist.
LG und danke für die Hilfe