Ansatz beim Plot

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Forum-Anfänger


	Beiträge: 20

	Anmeldedatum: 15.01.11

	Wohnort: ---

	Version: ---

Verfasst am: 15.01.2011, 02:55 Titel: Ansatz beim Plot

Moin, habe ein Problem. Ich bin gerade dabei mich mit Matlab anzufreunden. Und wollte das Matlab, ein paar Rechnungen für mich macht und diese dann Plotet jedoch habe ich ein kleines Problem.

Die Variablen
d_p_umg0, MesspDpU, dp60, F bestehen aus mehreren Werte die ich gemessen habe. Nun möchte ich mit denn Formeln alle für alle variablen einen wirkungsgrad , einen massenstrom ...
Die Ergebnisse würde ich dann gerne ploten.
nun stecke ich gerade bei
??? Error using ==> mpower
Matrix must be square.

Error in ==> labor at 72
expansionszahl=(kapa*((tau^(2/kapa))/(kapa-1))*(1-(tau^((kapa-1)/kapa)))/(1-tau))^0.5

fest. Help

Hier der m.code

Code:

clear all;
clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Variablen %

n=50; % Umdrehungen
p_umg=101330; % Umgebungsdruck In Pa

% Messp1pU in Pa CH03
d_p_umg0=[
2.6911
2.5937];

% Messwert in Pa (MesspDpU) CH02
MesspDpU=[
1.1207
1.0749];

% Messp3p1 CH04
dp60=[
1.0929
1.2715];

ra=0.15; % Radius außen in m
ri=0.11; % Radius innen in m
T_umg=294.3787; % Temperatur der Umgebung
Mv=0.22; % Störmoment
l=0.55; % Hebelarm in m

% Kraft
F=[
1.2866
1.3249];
alpha=0.992;
kapa=1.4;
R=287.058;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Rechnungen %

% Kalibrierung CH02
out=((15/2.986)+(20/3.989)+(25/4.994)+(30/6.001)+(35/7.007)+(40/8.015)+(45/9.034)+(50/10.044))/8;
fakCH02=out

% Kalibrierung CH03
out1=((15/3.015)+(20/4.021)+(25/5.028)+(30/6.036)+(35/7.044)+(40/8.051)+(45/9.060)+(50/10.068))/8;
fakCH03=out1

% Kalibrierung CH04
out2=((15/2.990)+(20/3.994)+(25/4.998)+(30/6.007)+(35/7.016)+(40/8.031)+(45/9.049)+(50/10.066))/8;
fakCH04=out2

% Kalibrierung CH05
out3=(-(1/1.027)-(2/2.052)-(3/3.078)-(4/4.106)-(5/5.132)-(6/6.159)-(7/7.183)-(8/8.210))/8;
fakCH05=out3*9.81

for iter=0:2

% Umrechnung
MesspDpU(:,1)=MesspDpU(:,1)*fakCH02
d_p_umg0(:,1)=d_p_umg0(:,1)*fakCH03
dp60(:,1)=dp60(:,1)*fakCH04
F(:,1)=F*fakCH05

% von mbar in Pa
MesspDpU(:,1)=MesspDpU(:,1)*100;
d_p_umg0(:,1)=d_p_umg0(:,1)*100;
dp60(:,1)=dp60(:,1)*100;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

tau(:,1)=1-MesspDpU(:,1)/p_umg
expansionszahl=(kapa*((tau^(2/kapa))/(kapa-1))*(1-(tau^((kapa-1)/kapa)))/(1-tau))^0.5
Ad=pi*(ra^2)
d_umg=p_umg/(R*T_umg)

volumenstromD=alpha*expansionszahl*Ad*sqrt((2*MesspDpU/d_umg))

p0=p_umg-d_p_umg0
d0=d_umg*(p0/p_umg)^(1/kapa)
volumenstrom0=(d_umg*volumenstromD)/d0
ua=2*pi*n*0.2
lieferzahl=volumenstrom0/(pi*ua*(0.2^2-0.11^2))

F=abs(F);
Mw=l*F
Pv=2*pi*n*Mv
Pw=2*pi*n*Mw
Pi=Pw-Pv

p6=dp60+p0
T0=p0/(R*d0)

Stuzenarbeit =(1.4/(1.4-1))*287.058*T0*((p6/p0)^((1.4-1)/1.4) -1)
massenstrom=volumenstromD*d_umg
Yi=Pi/massenstrom
Druckzahl=(2*Stuzenarbeit)/ua^2

wirkungsgrad=massenstrom*Stuzenarbeit/Pi
Leistungszahl=(lieferzahl*Druckzahl)/wirkungsgrad

end

Funktion ohne Link?

Sco

Forum-Meister


	Beiträge: 699

	Anmeldedatum: 15.08.10

	Wohnort: Dundee

	Version: 2008a, 2010a

Verfasst am: 15.01.2011, 03:45 Titel:

Hallo,

wenn du elementweise multiplizieren, potenzieren, oder im allgemeinen eine Rechenoperation ausführen möchtest musst du den Punktoperator verwenden. Siehe dieses Beispiel:

Code:

vec = [1 1; 2 2; 3 3];
vec^2 % = Matrixmultiplikation-->Fehler, da Matrixdimensionen nicht stimmen
vec.^2 % jedes Element wird potenziert-->Richtig

Funktion ohne Link?

Also müsste deine Zeile folgendermaßen aussehen:

Code:

expansionszahl=(kapa*((tau.^(2/kapa))/(kapa-1)).*(1-(tau.^((kapa-1)/kapa)))/(1-tau)).^0.5

Funktion ohne Link?

Dieser Fehler tritt auch noch in darauffolgenden Zeilen deines Codes auf.

MFG

Sco

Themenstarter

Forum-Anfänger


	Beiträge: 20

	Anmeldedatum: 15.01.11

	Wohnort: ---

	Version: ---

Verfasst am: 15.01.2011, 13:55 Titel:

vielen dank habe es geschafft, interessanter weise geht es ohne scheife

Code:

clear all;
clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Variablen %

n=50; % Umdrehungen
p_umg=101330; % Umgebungsdruck In Pa

% Messp1pU in Pa CH03
d_p_umg0=[
2.6911
2.5937];

% Messwert in Pa (MesspDpU) CH02
MesspDpU=[
1.1207
1.0749];

% Messp3p1 CH04
dp60=[
1.0929
1.2715];

ra=0.15; % Radius außen in m
ri=0.11; % Radius innen in m
T_umg=294.3787; % Temperatur der Umgebung
Mv=0.22; % Störmoment
l=0.55; % Hebelarm in m

% Kraft
F=[
1.2866
1.3249];
alpha=0.992;
kapa=1.4;
R=287.058;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Rechnungen %

% Kalibrierung CH02
out=((15/2.986)+(20/3.989)+(25/4.994)+(30/6.001)+(35/7.007)+(40/8.015)+(45/9.034)+(50/10.044))/8;
fakCH02=out

% Kalibrierung CH03
out1=((15/3.015)+(20/4.021)+(25/5.028)+(30/6.036)+(35/7.044)+(40/8.051)+(45/9.060)+(50/10.068))/8;
fakCH03=out1

% Kalibrierung CH04
out2=((15/2.990)+(20/3.994)+(25/4.998)+(30/6.007)+(35/7.016)+(40/8.031)+(45/9.049)+(50/10.066))/8;
fakCH04=out2

% Kalibrierung CH05
out3=(-(1/1.027)-(2/2.052)-(3/3.078)-(4/4.106)-(5/5.132)-(6/6.159)-(7/7.183)-(8/8.210))/8;
fakCH05=out3*9.81

%for iter=0:2

% Umrechnung
MesspDpU(:,1)=MesspDpU(:,1).*fakCH02
d_p_umg0(:,1)=d_p_umg0(:,1).*fakCH03
dp60(:,1)=dp60(:,1).*fakCH04
F(:,1)=F(:,1).*fakCH05

% von mbar in Pa
MesspDpU(:,1)=MesspDpU(:,1).*100;
d_p_umg0(:,1)=d_p_umg0(:,1).*100;
dp60(:,1)=dp60(:,1).*100;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

tau(:,1)=1-MesspDpU(:,1)./p_umg
expansionszahl(:,1)=(kapa*((tau(:,1).^(2/kapa))/(kapa-1)).*(1-(tau(:,1).^((kapa-1)/kapa)))./(1-tau(:,1))).^0.5

Ad=pi*(ra^2)
d_umg=p_umg/(R*T_umg)

volumenstromD(:,1)=alpha*expansionszahl(:,1)*Ad.*sqrt((2*MesspDpU(:,1)/d_umg))

p0(:,1)=p_umg(:,1)-d_p_umg0(:,1)
d0(:,1)=d_umg(:,1)*(p0(:,1)./p_umg(:,1)).^(1/kapa)
volumenstrom0(:,1)=(d_umg*volumenstromD(:,1))./d0(:,1)
ua=2*pi*n*0.2
lieferzahl(:,1)=volumenstrom0(:,1)./(pi*ua*(0.2^2-0.11^2))

F(:,1)=abs(F(:,1));
Mw(:,1)=l.*F(:,1)
Pv=2*pi*n*Mv
Pw(:,1)=2*pi*n*Mw(:,1)
Pi(:,1)=Pw(:,1)-Pv

p6(:,1)=dp60(:,1)+p0(:,1)
T0(:,1)=p0(:,1)./(R*d0(:,1))

Stuzenarbeit(:,1) =(1.4/(1.4-1))*287.058*T0(:,1).*((p6(:,1)./p0(:,1)).^((1.4-1)./1.4) -1)
massenstrom(:,1)=volumenstromD(:,1)*d_umg(:,1)
Yi(:,1)=Pi(:,1)./massenstrom(:,1)
Druckzahl(:,1)=(2*Stuzenarbeit(:,1))./ua^2

wirkungsgrad(:,1)=massenstrom(:,1).*Stuzenarbeit(:,1)./Pi
Leistungszahl(:,1)=(lieferzahl(:,1).*Druckzahl(:,1))./wirkungsgrad(:,1)

Funktion ohne Link?

jetzt muß ich die sachen nur noch ploten, also Stuzenarbeit zum Massenstrom.
ich glaube ich mu0 es so machen plot(Stuzenarbeit(:,1),Massenstrom(:,1))
gibt es eine funktion die gleich die Regression zeichnet ?

Sco

Forum-Meister


	Beiträge: 699

	Anmeldedatum: 15.08.10

	Wohnort: Dundee

	Version: 2008a, 2010a

Verfasst am: 15.01.2011, 15:55 Titel:

Hallo,

es gibt soweit ich weiss keine Funktion die die Regression automatisch, anhand der x- und y-Daten, mit in den plot zeichnet (zumal ja auch entscheident ist was für eine Art Regression du möchtest (Type,Grad)). Du kannst jedoch mit dem

Code:

hold on

Funktion ohne Link?

Befehl Daten im Plot zusätzlich mit darstellen (überlagern).

Für die Regression gibt es eine Fülle von Befehlen, z.B:

Code:

fit
polyfit

Funktion ohne Link?

Ein Beispiel wäre z.B. hier: http://www.gomatlab.de/gerade-durch.....07,highlight,sco+fit.html

Ansonsten bitte die Forumsuche verwenden, da dieses Thema schon zu genüge besprochen wurde.

MFG

Sco

Themenstarter

Forum-Anfänger


	Beiträge: 20

	Anmeldedatum: 15.01.11

	Wohnort: ---

	Version: ---

Verfasst am: 15.01.2011, 17:09 Titel:

ok vielen dank, läuft jedoch verstehe ich einfach nicht warum er beim zweiten Graphen keine Regression machen möchte... es ist bestimmt ein ganz dummer Fehler ......sitze schon seit 2h dran

Code:

clear all;
clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Variablen %

p_umg=[101330 101320]; % Umgebungsdruck In Pa
ra=0.15; % Radius außen in m
ri=0.11; % Radius innen in m
T_umg=[294.3787 294.2925]; % Temperatur der Umgebung
Mv=0.22; % Störmoment
l=0.55; % Hebelarm in m
alpha=0.992; % Durchflusszahl
kapa=1.4;
R=287.058; % Gaskonstante

% Umdrehungen Hz
n=[50 75];

% Messp1pU in Pa CH03
d_p_umg0=[
2.6911 6.1967
2.5937 5.1927
2.4820 4.9056
2.3676 4.6496
2.2937 4.3050
2.2151 4.0999
2.1566 3.8271
2.0001 3.5551
1.7946 3.3233
1.6438 3.1069
1.5558 2.8609
1.4653 2.5886
1.3395 2.4937
1.2722 2.3001
1.2002 2.0892
1.0726 1.9140
0.9456 1.7586
0.8828 1.5949
0.7586 1.5536
0.7526 1.4887
0.6331 0];

% Messwert in Pa (MesspDpU) CH02
MesspDpU=[
1.1207 2.5587
1.0749 2.1623
1.0324 2.0319
0.9841 1.9310
0.9523 1.7878
0.9210 1.7031
0.8995 1.5873
0.8324 1.4798
0.7412 1.3808
0.6802 1.2889
0.6452 1.1883
0.6076 1.0761
0.5510 1.0264
0.5196 0.9535
0.4907 0.8606
0.4389 0.7920
0.3882 0.6893
0.3623 0.5840
0.2874 0.5385
0.2808 0.5439
0.2302 0];

% Messp3p1 CH04
dp60=[
1.0929 2.0997
1.2715 3.6836
1.4229 4.0688
1.5817 4.3754
1.6755 4.7896
1.7786 4.9893
1.8446 5.2160
2.0651 5.4068
2.2616 5.5539
2.3587 5.6657
2.4233 5.7842
2.4667 5.8825
2.5323 5.9230
2.5662 5.9696
2.5934 5.9769
2.6435 5.9558
2.6619 4.6825
2.6615 4.6507
2.0832 4.5935
2.0957 4.5516
2.0782 0];

% Kraft
F=[
1.2866 2.7487
1.3249 3.1211
1.3651 3.1944
1.4003 3.2398
1.4173 3.2936
1.4341 3.3040
1.4452 3.3136
1.4939 3.3079
1.5057 3.2631
1.4925 3.2487
1.4903 3.2196
1.4792 3.1562
1.4556 3.1401
1.4486 3.0975
1.4324 3.0361
1.4046 2.9711
1.3698 2.6170
1.3525 2.5207
1.1655 2.5027
1.1636 2.4738
1.1101 0];

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Rechnungen %

% Kalibrierung CH02
out=((15/2.986)+(20/3.989)+(25/4.994)+(30/6.001)+(35/7.007)+(40/8.015)+(45/9.034)+(50/10.044))/8;
fakCH02=out

% Kalibrierung CH03
out1=((15/3.015)+(20/4.021)+(25/5.028)+(30/6.036)+(35/7.044)+(40/8.051)+(45/9.060)+(50/10.068))/8;
fakCH03=out1

% Kalibrierung CH04
out2=((15/2.990)+(20/3.994)+(25/4.998)+(30/6.007)+(35/7.016)+(40/8.031)+(45/9.049)+(50/10.066))/8;
fakCH04=out2

% Kalibrierung CH05
out3=(-(1/1.027)-(2/2.052)-(3/3.078)-(4/4.106)-(5/5.132)-(6/6.159)-(7/7.183)-(8/8.210))/8;
fakCH05=out3*9.81

for i=1:2

% Umrechnung
MesspDpU(:,i)=MesspDpU(:,i).*fakCH02;
d_p_umg0(:,i)=d_p_umg0(:,i).*fakCH03;
dp60(:,i)=dp60(:,i).*fakCH04;
F(:,i)=F(:,i).*fakCH05;

% von mbar in Pa
MesspDpU(:,i)=MesspDpU(:,i).*100;
d_p_umg0(:,i)=d_p_umg0(:,i).*100;
dp60(:,i)=dp60(:,i).*100;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

tau(:,i)=1-MesspDpU(:,i)./p_umg(:,i);
expansionszahl(:,i)=(kapa*((tau(:,i).^(2/kapa))/(kapa-1)).*(1-(tau(:,i).^((kapa-1)/kapa)))./(1-tau(:,i))).^0.5;

Ad=pi*(ra^2);
d_umg(:,i)=p_umg(:,i)./(R*T_umg(:,i));

volumenstromD(:,i)=alpha*expansionszahl(:,i)*Ad.*sqrt((2*MesspDpU(:,i)/d_umg(:,i)));

p0(:,i)=p_umg(:,i)-d_p_umg0(:,i);
d0(:,i)=d_umg(:,i)*(p0(:,i)./p_umg(:,i)).^(1/kapa);
volumenstrom0(:,i)=(d_umg(:,i).*volumenstromD(:,i))./d0(:,i);
ua(:,i)=2*pi*n(:,i)*0.2;
lieferzahl(:,i)=volumenstrom0(:,i)./(pi*ua(:,i).*(0.2^2-0.11^2));

F(:,i)=abs(F(:,i));
Mw(:,i)=l.*F(:,i);
Pv(:,i)=2*pi*n(:,i)*Mv;
Pw(:,i)=2*pi*n(:,i)*Mw(:,i);
Pi(:,i)=Pw(:,i)-Pv(:,i);

p6(:,i)=dp60(:,i)+p0(:,i);
T0(:,i)=p0(:,i)./(R*d0(:,i));

Stuzenarbeit(:,i) =(1.4/(1.4-1))*287.058*T0(:,i).*((p6(:,i)./p0(:,i)).^((1.4-1)./1.4) -1);
massenstrom(:,i)=volumenstromD(:,i)*d_umg(:,i);
Yi(:,i)=Pi(:,i)./massenstrom(:,i);
Druckzahl(:,i)=(2*Stuzenarbeit(:,i))./ua(:,i)^2;

wirkungsgrad(:,i)=massenstrom(:,i).*Stuzenarbeit(:,i)./Pi(:,i);
Leistungszahl(:,i)=(lieferzahl(:,i).*Druckzahl(:,i))./wirkungsgrad(:,i);

end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%% Ergebnisse %%%%%%%%%%%%%%%%%%%
MesspDpU
d_p_umg0
dp60
F
MesspDpU
d_p_umg0
dp60
tau
expansionszahl
Ad
d_umg
volumenstromD
p0
d0
volumenstrom0
ua
lieferzahl
F
Mw
Pv
Pw
Pi
p6
T0
Stuzenarbeit
massenstrom
Yi
Druckzahl
wirkungsgrad
Leistungszahl
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%% Plot %%%%%%%%%%%%%%%%%%%
%figure
grid
hold on
p=polyfit(volumenstrom0(:,1),Stuzenarbeit(:,1),4);
f=polyval(p,volumenstrom0(:,1));
plot(volumenstrom0(:,1),Stuzenarbeit(:,1),'o',volumenstrom0(:,1),f,'-','Color','red','MarkerFacecolor','green')

p=polyfit(volumenstrom0(:,2),Stuzenarbeit(:,2),4);
f=polyval(p,volumenstrom0(:,2));
plot(volumenstrom0(:,2),Stuzenarbeit(:,2),'x',volumenstrom0(:,2),f,'-','Color','blue','MarkerFacecolor','green')

xlabel('Volumenstrom [m^3/s]')
ylabel('Stuzenarbeit [m^2/s^2]')
legend('50Hz','Location','NorthEast');

Funktion ohne Link?

Sco

Forum-Meister


	Beiträge: 699

	Anmeldedatum: 15.08.10

	Wohnort: Dundee

	Version: 2008a, 2010a

Verfasst am: 15.01.2011, 17:51 Titel:

Hallo,

das Probem ist, dass dein Vektor "(volumenstrom02(:,2)" ein NaN enthält, so dass polyfit als Ergebnis nur NaN's ausgibt.

Wie du das lösen möchtest musst du wohl selber lösen. Du könntest z.B. diesen Wert durch die benachbarten interpolieren, oder für den fit diesen einen Wert löschen und mit den restlichen Werten die Regression berechnen. Letzteres würde z.B. so funktionieren:

Code:

hold on
p=polyfit(volumenstrom0(:,1),Stuzenarbeit(:,1),4);
f=polyval(p,volumenstrom0(:,1));
plot(volumenstrom0(:,1),Stuzenarbeit(:,1),'o',volumenstrom0(:,1),f,'-','Color','red','MarkerFacecolor','green')

volumenstrom02 = volumenstrom0;
Stuzenarbeit2 = Stuzenarbeit;
[row,col] = find(isnan(volumenstrom02));
volumenstrom02(row,:) = [];
Stuzenarbeit2(row,:) = [];

p=polyfit(volumenstrom02(:,2),Stuzenarbeit2(:,2),4);
f=polyval(p,volumenstrom02(:,2));
plot(volumenstrom0(:,2),Stuzenarbeit(:,2),'x',volumenstrom02(:,2),f,'-','Color','blue','MarkerFacecolor','green')

Funktion ohne Link?

MFG

Sco

Themenstarter

Forum-Anfänger


	Beiträge: 20

	Anmeldedatum: 15.01.11

	Wohnort: ---

	Version: ---

Verfasst am: 15.01.2011, 18:00 Titel:

aaah verstehe, könnte ich nicht auch einfach bei durchrechen die 0 wecklassen

das heist nicht Stuzenarbeit(:,i)=.......

sondern Stuzenarbeit(von 0 bis 16,i)...

wollte es auch so machen, habe es aber nicht hinbekommen. Dachte so an
Stuzenarbeit(0:16,i)... geht leider nicht.

Sco

Forum-Meister


	Beiträge: 699

	Anmeldedatum: 15.08.10

	Wohnort: Dundee

	Version: 2008a, 2010a

Verfasst am: 15.01.2011, 18:06 Titel:

Hallo,

aus deinem letzten Posting werde leider absolut nicht schlau. Ist deine Frage nun beantwortet oder verbleiben Fragen?

Zitat:

könnte ich nicht auch einfach bei durchrechen die 0 wecklassen

Welche 0 denn?

Zitat:

das heist nicht Stuzenarbeit(:,i)=.......

sondern Stuzenarbeit(von 0 bis 16,i)...

wollte es auch so machen, habe es aber nicht hinbekommen. Dachte so an
Stuzenarbeit(0:16,i)... geht leider nicht.

Verstehe ich nicht. Worauf bezieht sich das denn?

MFG

Sco

Themenstarter

Forum-Anfänger


	Beiträge: 20

	Anmeldedatum: 15.01.11

	Wohnort: ---

	Version: ---

Verfasst am: 15.01.2011, 18:17 Titel:

Es ist eigentlich nur noch eine Verständnis Frage.

Ganz am Anfang habe ich eine Matrix 21x2 Definiert
Die d_p_umg0
Im Grunde brauche ich denn Wert 0 ganz am Ende also bei d_p_umg0(21,2) nicht bei der Berechnung. Kann ich bei der Definition nicht etwas anderes Hinschreiben ( end ) oder so, wodurch dieser Wert bei der Berechnung keine Rolle mehr spielt.
lg

Themenstarter

Forum-Anfänger


	Beiträge: 20

	Anmeldedatum: 15.01.11

	Wohnort: ---

	Version: ---

Verfasst am: 15.01.2011, 19:10 Titel:

Edit: Ok ich komme echt nicht mehr weiter es kann doch nicht so schwer sein eine Legende zu ändern. Er macht es bei mir nacheinander. Zu erst Punkt 1 , dann Regression, dann Punkt zwei. Ich habe alles Probiert....
Help....

und eine kleinigkeit, was muß ich tun damit Regression nicht mehr in der Legende angezeigt wird ?

Code:

figure
grid
hold on
p=polyfit(volumenstrom0(:,1),Yi(:,1),4);
f=polyval(p,volumenstrom0(:,1));
ha=plot(volumenstrom0(:,1),Yi(:,1),'o',volumenstrom0(:,1),f,'-','Color','red','MarkerFacecolor','green')

volumenstrom02 = volumenstrom0;
Yi2 = Yi;
[row,col] = find(isnan(volumenstrom02));
volumenstrom02(row,:) = [];
Yi2(row,:) = [];

p=polyfit(volumenstrom02(:,2),Yi2(:,2),4);
f=polyval(p,volumenstrom02(:,2));
hb=plot(volumenstrom0(:,2),Yi(:,2),'x',volumenstrom02(:,2),f,'-','Color','blue','MarkerFacecolor','green')

xlabel('Volumenstrom [m^3/s]')
ylabel('innere Stuzenarbeit [m^2/s^2]')
hleg1 = legend('50Hz','Regresion 1','75 Hz');
set(hleg1,'Location','NorthEast')
set(hleg1,'Interpreter','none')
set(hleg1,'FontAngle','italic','TextColor',[.3 .2 .1])

Funktion ohne Link?

Sco

Forum-Meister


	Beiträge: 699

	Anmeldedatum: 15.08.10

	Wohnort: Dundee

	Version: 2008a, 2010a

Verfasst am: 15.01.2011, 20:22 Titel:

Hallo,

Zitat:

und eine kleinigkeit, was muß ich tun damit Regression nicht mehr in der Legende angezeigt wird ?

um dies zu realisieren musst du deine plot Struktur aufteilen und separat handles vergeben:

Code:

figure
hax = axes;
set(hax,'NextPlot','add'); % ersetzt den "hold on" Befehl und ist schneller
grid

p=polyfit(volumenstrom0(:,1),Yi(:,1),4);
f=polyval(p,volumenstrom0(:,1));
ha=plot(volumenstrom0(:,1),Yi(:,1),'o','Color','red','MarkerFacecolor','green');
plot(volumenstrom0(:,1),f,'-','Color','red','MarkerFacecolor','green');

volumenstrom02 = volumenstrom0;
Yi2 = Yi;
[row,col] = find(isnan(volumenstrom02));
volumenstrom02(row,:) = [];
Yi2(row,:) = [];

p=polyfit(volumenstrom02(:,2),Yi2(:,2),4);
f=polyval(p,volumenstrom02(:,2));
hb=plot(volumenstrom0(:,2),Yi(:,2),'x','Color','blue','MarkerFacecolor','green');
plot(volumenstrom02(:,2),f,'-','Color','blue','MarkerFacecolor','green');

xlabel('Volumenstrom [m^3/s]')
ylabel('innere Stuzenarbeit [m^2/s^2]')
hleg1 = legend([ha hb],'50Hz','75 Hz');
set(hleg1,'Location','NorthEast')
set(hleg1,'Interpreter','none')
set(hleg1,'FontAngle','italic','TextColor',[.3 .2 .1])

Funktion ohne Link?

Eine goolge Suche hätte dich auch schnell an das Ziel geführt!

Du könntest in deinem Code auch folgende Zeilen am ENDE einfügen um es etwas drekig zu basteln (jedoch ist die erste Lösung viel besser, da hier nur die Linie in der Legend Box gelöscht wird, aber die Legendbox gleich groß bleibt):

Code:

hleg1 = legend('50Hz','','75 Hz'); % Regressinsstring ist leer
set(hleg1,'Location','NorthEast')
set(hleg1,'Interpreter','none')
set(hleg1,'FontAngle','italic','TextColor',[.3 .2 .1])

hc = get(hleg1,'Child');
set(hc(5),'XData',[0 0]);

Funktion ohne Link?

MFG

Sco


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.