Maximale positiven UND negativen Wert einer Matrix

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

authrey

Forum-Newbie


	Beiträge: 3

	Anmeldedatum: 25.02.18

	Wohnort: ---

	Version: ---

Verfasst am: 25.02.2018, 16:11 Titel: Maximale positiven UND negativen Wert einer Matrix

Servus,

ich würde gerne aus mehreren Matrizen die maximalen Eintrage ausgeben lassen. Mein
Problem ist jedoch, dass in diesen Matrizen auch negative Werte vorkommen. Anstatt den größten negativen Wert in die Matrix zu schreiben bekomme ich eine null mit dem Befehl max(A) (ist ja auch logisch weil 0 in diesem Fall der höchste Wert ist). Für meinen Anwendungsfall benötige ich aber unabhängig vom Vorzeichen den höchsten Wert.

Hat jemand eine Idee?

Nutze zurzeit den unteren Befehl.

Code:

f_max(k,j)=max(f(1,:))

Funktion ohne Link?

huhu0815

Forum-Century


	Beiträge: 141

	Anmeldedatum: 20.02.14

	Wohnort: Dresden

	Version: R2013b, R2014b, R2016b, R2018a

Verfasst am: 25.02.2018, 16:55 Titel:

Hi,

unabhängig vom vorzeichen, kannst du abs verwenden.

Code:

A = -10 * rand(3, 2);
mitVorzeichen= max(A);
ohneVorzeichen = max(abs(A));

Funktion ohne Link?

_________________

Viele Grüße

authrey

Themenstarter

Forum-Newbie


	Beiträge: 3

	Anmeldedatum: 25.02.18

	Wohnort: ---

	Version: ---

Verfasst am: 25.02.2018, 22:57 Titel:

Vielen Dank für deinen Tipp!
Mir werden jetzt alle Werte Ausgespuckt.

Gibt es eine Möglichkeit mir aus die maximalen negativen Werte in die Matrix zu schreiben? Ich möchte die Matrix plotten, und wenn die zuvor negativen Werte jetzt positiv sind habe ich Sinngemäß ein Problem.
Den kompletten Code habe ich angehängt, vielleicht hilft es weiter.

Danke nochmal für deine Hilfe!

Code:

clear all
clc

alpha=18.5;
b=2.8387;
a=1.8316;
l=6.8711;
Fp=200;
q=linspace(10,160,16);
E=210888;
h=linspace(0.4,1,7);
c=0.402;

c_1=-4*(Fp/(b^2));
c_2=4*((Fp*(2*a+b))/(b^2));
c_3=-4*((Fp*a*(a+b))/(b^2));

for j=1:length(q)
for k=1:length(h)
A=[-(a^2)/2 -a 0 0 0 0 0 0 (a^2)/2 a 1 0; -(a^3)/6 -(a^2)/2 0 0 0 0 0 0 (a^3)/6 (a^2)/2 a 1; 0 0 0 0 -((a+b)^2)/2 -(a+b) -1 0 ((a+b)^2)/2 (a+b) 1 0; 0 0 0 0 -((a+b)^3)/6 -((a+b)^2)/2 -(a+b) -1 ((a+b)^3)/6 ((a+b)^2)/2 (a+b) 1; -1 0 0 0 0 0 0 0 1 0 0 0; -a -1 0 0 0 0 0 0 a 1 0 0; 0 0 0 0 -1 0 0 0 1 0 0 0; 0 0 0 0 -(a+b) -1 0 0 (a+b) 1 0 0; 0 0 0 0 (l^3)/6 (l^2)/2 l 1 0 0 0 0; 0 0 0 0 (l^2)/2 l 1 0 0 0 0 0; 0 0 1 0 0 0 0 0 0 0 0 0; 0 0 0 1 0 0 0 0 0 0 0 0];
A_inv=inv(A);
B=[-(1/60)*c_1*(a^5)-(1/24)*c_2*(a^4)-(1/6)*c_3*(a^3) ; -(1/360)*c_1*(a^6)-(1/120)*c_2*(a^5)-(1/24)*c_3*(a^4) ; -(1/60)*c_1*((a+b)^5)-(1/24)*c_2*((a+b)^4)-(1/6)*c_3*((a+b)^3) ; -(1/360)*c_1*((a+b)^6)-(1/120)*c_2*((a+b)^5)-(1/24)*c_3*((a+b)^4) ; -(1/3)*c_1*(a^3)-(1/2)*c_2*(a^2)-c_3*a ; -(1/12)*c_1*(a^4)-1/6*c_2*(a^3)-(1/2)*c_3*(a^2) ; -(1/3)*c_1*((a+b)^3)-(1/2)*c_2*((a+b)^2)-c_3*(a+b) ; -(1/12)*c_1*((a+b)^4)-(1/6)*c_2*((a+b)^3)-(1/2)*c_3*((a+b)^2) ; (1/24)*q(j)*(l^4) ; (1/6)*q(j)*(l^3) ; 0 ; 0];
C=A_inv*B;
C_1=C(1,1);
C_2=C(2,1);
C_3=C(3,1);
C_4=C(4,1);
C_5=C(5,1);
C_6=C(6,1);
C_7=C(7,1);
C_8=C(8,1);
C_9=C(9,1);
C_10=C(10,1);
C_11=C(11,1);
C_12=C(12,1);

I=(c*2)*(h(k)^3)/12;
s=linspace(0,l,10000);

for i=1:length(s)
x=s(1,i);
if x <= a
Q(1,i)=q(j)*x-C_1;
elseif x > a && x <= (a+b)
Q(1,i)=q(j)*x-(1/3)*c_1*(x^3)-(1/2)*c_2*(x^2)-c_3*x-C_9;
else
Q(1,i)=q(j)*x-C_5;
end
i=i+1;
end

for i=1:length(s)
x=s(1,i);
if x <= a
M(1,i)=(1/2)*q(j)*(x^2)-C_1*x-C_2;
elseif x > a && x <= (a+b)
M(1,i)=(1/2)*q(j)*(x^2)-(1/12)*c_1*(x^4)-(1/6)*c_2*(x^3)-(1/2)*c_3*(x^2)-C_9*x-C_10;
else
M(1,i)=(1/2)*q(j)*(x^2)-C_5*x-C_6;
end
i=i+1;
end

for i=1:length(s)
x=s(1,i);
if x <= a
f(1,i)=(1/(E*I))*(-(1/24)*q(j)*(x^4)+(1/6)*C_1*(x^3)+(1/2)*C_2*(x^2)+C_3*x+C_4);
f(2,i)=x;
elseif x > a && x <= (a+b)
f(1,i)=(1/(E*I))*(-(1/24)*q(j)*(x^4)+(1/360)*c_1*(x^6)+(1/120)*c_2*(x^5)+(1/24)*c_3*(x^4)+(1/6)*C_9*(x^3)+(1/2)*C_10*(x^2)+C_11*x+C_12);
f(2,i)=x;
else
f(1,i)=(1/(E*I))*(-(1/24)*q(j)*(x^4)+(1/6)*C_5*(x^3)+(1/2)*C_6*(x^2)+C_7*x+C_8);
f(2,i)=x;
end
i=i+1;
end

sigma_vergleich_max(k,j)=sqrt((h(k)/(I*2)*max(M))^2+3*(c*h(k))*max(Q))/1.168;

f_max(k,j)=max(abs(f(1,:)));

k=k+1;
end
j=j+1;
end

for l=1:length(h)
subplot(2,1,1)
plot(q(1,:),sigma_vergleich_max(l,:))
hold on
text(161,max(sigma_vergleich_max(l,:)), num2str(h(1,l)))

subplot(2,1,2)
plot(q(1,:),f_max(l,:))
hold on
text(161,max(f_max(l,:)), num2str(h(1,l)))
end

Funktion ohne Link?

authrey

Themenstarter

Forum-Newbie


	Beiträge: 3

	Anmeldedatum: 25.02.18

	Wohnort: ---

	Version: ---

Verfasst am: 26.02.2018, 11:39 Titel:

Habs gelöst, danke!
Für die die es Interessiert, Code ist unten.

Code:

if max(f(1,:)) ~= 0
f_max(k,j)=-max(f(1,:));
else
f_max(k,j)=-min(f(1,:));
end

Funktion ohne Link?


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 - 2024 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.