while Schleife

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

gydoh

Forum-Anfänger


	Beiträge: 33

	Anmeldedatum: 16.11.15

	Wohnort: ---

	Version: ---

Verfasst am: 21.09.2016, 09:35 Titel: while Schleife

Hallo liebe Leute,

ich habe ein Problem mit einer while Schleife (oder vielleicht kann man es auch anders lösen).
Hier mein Code:

Code:

%% Clear
clc
clear all
close all

%% Data for the Tjaerborg WT
omega = 2.3; % [s^-1] rotational speed
rho_air = 1.225; % [kg/m^3] density of air
R = 30.56; % [m] Rotor radius
B = 3; % [-] Number of blades
P_rate = 2000; % [kW] rated Power
U_in = 5; % [m/s] Cut in speed
U_out = 25; % [m/s] Cut out speed
r_i = [6.46 9.46 12.46 15.46 18.46 21.46 24.46 27.46 28.96 29.86 30.55]; % Point on the blade
c = [3.3 3.0 2.7 2.4 2.1 1.8 1.5 1.2 1.05 0.96 0.95]; % Chord length
beta = [8 7 6 5 4 3 2 1 0.5 0.2 0.16]; % Twisting angle
theta_p = 0; % Pitch angle

%% importing data
x1 = importdata('tjaere4.dat');[m,n]=size(x1);
alfatab = x1(:,1);
cltab = x1(:,2);
cdtab = x1(:,3);

%% BEM
a = 0;
a_prime = 0;

flow_angle = atan((1-a)*U_in)/((1+a_prime)*omega*r_i(1));
alpha_local = flow_angle*(180/pi)-(beta(1)+theta_p);
Cl = interp1(alfatab,cltab,alpha_local);
Cd = interp1(alfatab,cdtab,alpha_local);
Cn = Cl*cos(flow_angle)+Cd*sin(flow_angle);
Ct = Cl*sin(flow_angle)+Cd*cos(flow_angle);
sigma = (c(1)*B)/(2*pi*r_i(1));
a1 = 1/(((4*sin(flow_angle)^2)/(sigma*Cn))+1);
a1_prime = 1/(((4*sin(flow_angle)*cos(flow_angle))/(sigma*Ct))-1);

eps1 = abs(a1-a);
eps2 = abs(a1_prime-a_prime);

Funktion ohne Link?

Mein problem ist jetzt, dass ich a1 und a1_prime so lange berechnen muss, bis eps1 und eps2 < 10^-4 sind. Das heißt ich muss immer wieder bei dem Schritt anfangen in dem der flow_angle berechnet wird. Zu Beginn setze ich a und a_prime = 0 und berechne dann neue a und a_prime die ich a1 und a1_prime genannt habe. Und das so lange bis eps1 und eps2 < 10^-4 sind.
Habe es mit einer while schleife versucht, aber dann kann ich eps nicht mehr berechnen, da ich kein vorheriges a und a_prime mehr habe.

Ich hoffe jemand versteht mein Problem und kann mir helfen. Es ist etwas schwierig zu erkläre. Vielen Dank schonmal!

LG Gydoh

Harald

Forum-Meister


	Beiträge: 24.501

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 21.09.2016, 09:54 Titel:

Hallo,

ich sehe nicht ganz, wo das Problem liegt

Code:

% Initialisierung
eps1 = 1;
eps2 = 1;

while eps1 > 1e-4 | eps2 > 1e-4

flow_angle = atan((1-a)*U_in)/((1+a_prime)*omega*r_i(1));
alpha_local = flow_angle*(180/pi)-(beta(1)+theta_p);
Cl = interp1(alfatab,cltab,alpha_local);
Cd = interp1(alfatab,cdtab,alpha_local);
Cn = Cl*cos(flow_angle)+Cd*sin(flow_angle);
Ct = Cl*sin(flow_angle)+Cd*cos(flow_angle);
sigma = (c(1)*B)/(2*pi*r_i(1));
a1 = 1/(((4*sin(flow_angle)^2)/(sigma*Cn))+1);
a1_prime = 1/(((4*sin(flow_angle)*cos(flow_angle))/(sigma*Ct))-1);

eps1 = abs(a1-a);
eps2 = abs(a1_prime-a_prime);
a = a1;
a_prime = a1_prime;
end

Funktion ohne Link?

Testen konnte ich das aufgrund fehlender Daten nicht.

Grüße,
Harald


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.