Implementatie formule in matlab
-
- Berichten: 2
Implementatie formule in matlab
Hallo,
Ik heb een vraag met betrekking tot het implementeren van een functie met 2 variabelen in matlab en dan een plot hieruit krijgen.
To en deltaT zijn constanten;
De exacte formule is T(x,t)=To+deltaT*exp(-t)*cos(pi*x).
Ik ben al een hele dag hiermee in de weer en krijg rare resultaten die nergens op slaan, ik weet dat ik geen held ben in Matlab. Heeft iemand een idee hoe ik deze functie kan implementeren in Matlab. De plot moet eigenlijk een verbranding voorstellen van een gas. Dus er zou binnen een kort tijdsdomein, dus t tussen 0 en 3 seconden, een hoge Temperatuur T moeten uitkomen.
Alvast bedankt..
Ik heb een vraag met betrekking tot het implementeren van een functie met 2 variabelen in matlab en dan een plot hieruit krijgen.
To en deltaT zijn constanten;
De exacte formule is T(x,t)=To+deltaT*exp(-t)*cos(pi*x).
Ik ben al een hele dag hiermee in de weer en krijg rare resultaten die nergens op slaan, ik weet dat ik geen held ben in Matlab. Heeft iemand een idee hoe ik deze functie kan implementeren in Matlab. De plot moet eigenlijk een verbranding voorstellen van een gas. Dus er zou binnen een kort tijdsdomein, dus t tussen 0 en 3 seconden, een hoge Temperatuur T moeten uitkomen.
Alvast bedankt..
-
- Berichten: 4.246
Re: Implementatie formule in matlab
Zoiets als dit:
T0=
deltaT=
tb=
te=
xb=
xe=
t=tb:1:te
x=xb:1:xe
plot(t, To+deltaT*exp(-t)*cos(pi*x))
Volgens mij moeten er nog een paar punten in de formule vanwege vectorvermenigvuldiging.
T0=
deltaT=
tb=
te=
xb=
xe=
t=tb:1:te
x=xb:1:xe
plot(t, To+deltaT*exp(-t)*cos(pi*x))
Volgens mij moeten er nog een paar punten in de formule vanwege vectorvermenigvuldiging.
Quitters never win and winners never quit.
-
- Berichten: 2
Re: Implementatie formule in matlab
Bedankt voor het antwoord, uiteindelijk ben ik erachter gekomen wat het moest zijn:
figure
for k = [1:51];
T0=0.25;
x=[0:0.001:1];
DeltaT=0.1;
t=(k-1)*3/50;
plot(x,T0+((DeltaT*exp(-t))*cos(x*pi)));
ylim([0.15 0.35]);
title('Temperatuur tov x');
xlabel('x');
ylabel('T');
M(:,k) = getframe(gcf);
end
figure
for k = [1:51];
T0=0.25;
x=[0:0.001:1];
DeltaT=0.1;
t=(k-1)*3/50;
plot(x,T0+((DeltaT*exp(-t))*cos(x*pi)));
ylim([0.15 0.35]);
title('Temperatuur tov x');
xlabel('x');
ylabel('T');
M(:,k) = getframe(gcf);
end