Inner matrix dimension < Matlab < Mathe-Software < Mathe < Vorhilfe
|
Hallo zusammen.
Eigentlich weiss ich was der Fehler bedeutet: inner matrix dimension must agree!!
Aber ich weiss bei meinem Programm wirklich nicht was falsch ist. Mein Programm habe ich unten angegeben. Den Radius r habe ich in eine Matrix geschrieben und möchte nun die Kurve plotten. Bei einem ähnlichen Programm von mir hat es wunderbar geklappt, aber hier sehe ich nicht durch. Ich hoffe Jemand kann mir weiter helfen.
[mm] E=16.7*10^6;
[/mm]
[mm] p=1*10^5;
[/mm]
a=0.005;
h=0.0006;
r=0:0.0001:a;
y=0.5;
[mm] K=E*h^3/(12*(1-y^2));
[/mm]
%% Processing
% Auslenkung
[mm] w=(p*(r.^2-a^2)*(r.^2-a^2-6.4*h^2))/(64*K)
[/mm]
plot(r,w);
|
|
| |
|
moin,
r ist doch ein Vektor!
Im letzten Term tritt das r zweimal auf. Einmal müsstest du transponieren.
Grob hast du doch
[mm] $w=v_1(r)*v_2(r)/Rest$
[/mm]
Entweder benutzt du den Operator ".*" oder transponierst einen Faktor.
Als Tipp:
Zeile auskommentieren und statt einer Zeile lieber Zwischenergebnisse abspeichern für Debuggen.
|
|
|
| |
|
Danke für die Antwort. Habe es jetzt halt so gelöst:
[mm] c=(p*(r.^2-a^2));
[/mm]
[mm] d=(r.^2-a^2-6.4*h^2)/(64*K);
[/mm]
w=c.*d;
Nun hätte ich noch eine weitere Frage. Diese Funktion möchte ich gerne im Kreis drehen und dann das Volumen bestimmen. Das Integral wird ja so aufgestellt:
[mm] 2*pi*\integral_{0}^{a}{w(r)dr}
[/mm]
Ich habe noch nie im Matlab integriert, gibt es da einen Befehl?
|
|
|
| |
|
moin,
> Danke für die Antwort. Habe es jetzt halt so gelöst:
> [mm]c=(p*(r.^2-a^2));[/mm]
> [mm]d=(r.^2-a^2-6.4*h^2)/(64*K);[/mm]
> w=c.*d;
Warum schreibst du dann nicht
| 1: |
| | 2: | w=(p*(r.^2-a^2).*(r.^2-a^2-6.4*h^2))/(64*K);
|
gleich?
>
>
> Nun hätte ich noch eine weitere Frage. Diese Funktion
> möchte ich gerne im Kreis drehen und dann das Volumen
> bestimmen. Das Integral wird ja so aufgestellt:
>
> [mm]2*pi*\integral_{0}^{a}{w(r)dr}[/mm]
>
> Ich habe noch nie im Matlab integriert, gibt es da einen
> Befehl?
Ich auch nicht. Aber eine ![[]](/images/popup.gif) Dokumentation von trapz gibt es!
| 1: |
| | 2: | E=16.7*1e6;
| | 3: | p=1e5;
| | 4: | a=0.005;
| | 5: | h=0.0006;
| | 6: | r=0:.0001:a;
| | 7: | y=0.5;
| | 8: | K=E*h^3/(12*(1-y^2));
| | 9: | w=(p*(r.^2-a^2).*(r.^2-a^2-6.4*h^2))/(64*K);
| | 10: | 2*pi*trapz(r,w)
|
wieschoo
|
|
|
|
