Polynomauswertung

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

WiMa2014

Forum-Newbie


	Beiträge: 4

	Anmeldedatum: 02.11.14

	Wohnort: ---

	Version: ---

Verfasst am: 02.11.2014, 15:13 Titel: Polynomauswertung

Hallo Leute, ich bin neu hier und bin gerade dabei ein Programm zur Berechnung der reellen Lösungen der allgemeinen quadratischen Gleichung p(x)=ax^2+bx^2+c = 0 zu schreiben.

Dazu habe ich einige Beispiele gemacht. Aber mein Programm läuft nicht. Jetzt wollte ich mal fragen ob ihr einen Fehler erkennt.

Liebe Grüße

Hier mein kleines Programm:

Code:

av = [2 0 0 1 0 1 4 10]
bv = [10 0 1 0 0 -2 -12 10]
cv = [1 1 0 0 0 8 9 0.01]
k= 1 %Index

% Diese Vektoren haben wir erstellt, damit wir schonmal die Probe im Überblick
% haben und nur noch jedes mal einen Index weiter gehen müssen um die Probe
% durchzuführen. Jetzt wollen wir für p(x) = ax^2 + bx + c = 0 die reellen
% Lösungen.

while k <= length(av)

a= av(k)
b= bv(k)
c= cv(k)

if a~=0

p = b/a
q = c/a

if q < p^2/4

u = p^2/4
v = u - q
w = sqrt(v)

x2 = - p/2 - w
x1 = q/x2

disp('Es gibt zwei mögliche Nullstellen x1 = ',num2str(x1),'und x2 = ',num2str(x2))
% Damit lassen wir uns die möglichen Nullstellen ausgeben
end

else if b~=0

x=-c/b

disp(' Es gibt nur eine mögliche Nullstelle x = ' ,num2str(x))

else

disp (' Da a=0 und b=0 muss auch c=0 sein! Und es ist keine quadratische Gleichung vorhanden')

k = k+1

end
end

Funktion ohne Link?

[EDITED, Jan, Bitte Code-Umgebung verwenden - Danke!]

Seban

Forum-Meister


	Beiträge: 600

	Anmeldedatum: 19.01.12

	Wohnort: ---

	Version: ab R2014b

Verfasst am: 02.11.2014, 20:14 Titel:

Hallo,

Vorweg: Wenn du Befehle mit ; beendest, werden diese nicht im Command Window ausgegeben.
Und der Befehl heißt elseif, nicht else if. Dann stimmt auch die Anzahl der end (gab bei mir eine Fehlermeldung, ein end zu wenig).

k=k+1 ist falsch platziert, die while Schleife läuft unendlich lange weil k konstant 1 bleibt.

Hast du keine Fehlermeldungen erhalten?

Bei mir erscheint folgende Meldung (ebenfalls für Line 40):

Code:

Error using disp
Too many input arguments.

Error in TEMP (line 32)
disp('Es gibt zwei mögliche Nullstellen x1 = ',num2str(x1),'und x2 = ',num2str(x2))

Funktion ohne Link?

Code:

fprintf

Funktion ohne Link?

bietet sich als Alternative an an.

Ist es nich sinnvoll x, x1, x2 auch zu speichern anstatt sie nur auszugeben?

Code:

av = [2 0 0 1 0 1 4 10];
bv = [10 0 1 0 0 -2 -12 10];
cv = [1 1 0 0 0 8 9 0.01];
k= 1; %Index

% Diese Vektoren haben wir erstellt, damit wir schonmal die Probe im Überblick
% haben und nur noch jedes mal einen Index weiter gehen müssen um die Probe
% durchzuführen. Jetzt wollen wir für p(x) = ax^2 + bx + c = 0 die reellen
% Lösungen.

x = NaN(size(av)); % preallocate
x1 = NaN(size(av));
x2 = NaN(size(av));

while k <= length(av)

a= av(k);
b= bv(k);
c= cv(k);

if a~=0

p = b/a;
q = c/a;

if q < p^2/4

u = p^2/4;
v = u - q;
w = sqrt(v);

x2(k) = - p/2 - w;
x1(k) = q/x2(k);

% disp('Es gibt zwei mögliche Nullstellen x1 = ',num2str(x1),'und x2 = ',num2str(x2))
% Damit lassen wir uns die möglichen Nullstellen ausgeben
end

elseif b~=0

x(k)=-c/b;

% disp(' Es gibt nur eine mögliche Nullstelle x = ' ,num2str(x))

else

disp (' Da a=0 und b=0 muss auch c=0 sein! Und es ist keine quadratische Gleichung vorhanden')

end

k = k+1;
end

Funktion ohne Link?

Grüße,
Seban
_________________

Richtig fragen
Debugging

Nras

Forum-Meister


	Beiträge: 608

	Anmeldedatum: 21.02.12

	Wohnort: ---

	Version: 7.12.0.635 (R2011a)

Verfasst am: 03.11.2014, 00:22 Titel:

Hallo,
falls das nicht nur eine Fingerübung ist, und du wirklich Nullstellen von Polynomen suchst, dann würde man das vermutlich anders lösen. Die Polynome werden in Matlab durch einen Vektor mit den Koeffizienten repräsentiert, darauf kann man dann bestimmte Funktionen anwenden, wie zum beispiel polyval() oder roots().

Code:

f = @(x) 1*x.^2 - 1*x - 2; % Ein Polynom
p = [1, 1, -2]; % Polynomrepräsentation in Matlab
r = roots(p); % Nullstellen des Polynoms (-2 und 1)

ezplot(f) % Plot
hold on
plot(r, f(r), 'ro', 'markerfacecolor', 'red') % plot Nullstellen
% alternativ:
% plot(r, polyval(p, r), 'ks', 'markerfacecolor', 'red')

Funktion ohne Link?

Viele Grüße,
Nras.

WiMa2014

Themenstarter

Forum-Newbie


	Beiträge: 4

	Anmeldedatum: 02.11.14

	Wohnort: ---

	Version: ---

Verfasst am: 05.11.2014, 13:20 Titel:

Danke für die schnelle Hilfe. Hat mir einiges geholfen.
Liebe grüße


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.