Zaloguj się
Blog na Matlab.pl
Forum polskich użytkowników
 
UŻYTKOWNICY GRUPY PROFIL Zaloguj się, by sprawdzić wiadomości FAQ
 



Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Przepływ ciepła
Forum MATLAB Strona Główna-> MATLAB
Post Wysłany: 25 Grudnia 2016, Nie 12:43 pm Temat postu: Przepływ ciepła Odpowiedz z cytatem
 
AUTOR:
Fasol00005



Dołączył: 12 Paź 2015
Posty: 2


Ogląda profil użytkownika Wyślij prywatną wiadomość
Witam, przerabiam jeden z podręczników do nauku matlaba. Jestem na rodziale przepływ ciepła. Zapisany jest w podręczniku gotowy skrypt. Niestety nie działa. Prosze o pomoc.
k=1e-7;
n=10;
L=0.36;
T0 = 25;
TN = -15
z=ones(1,n-1)*(T0+TN)*0.25;
[t,T] = ode45(@rownanie,linspace(0,200000,100),z,[],n);
surf(T)
view([31 22]);
xlabel('warstwa - x')
ylabel('czas')
zlabel('temperatura')

function dy = rownanie(t,T,n)
dy=zeros(n-1,1);
global k L T0 TN
d=n*n*k/(L^2);
dy(1)=d*(T(2)-2*T(1)+T0);
for i=2:n-2
dy(i)=d*(T(i+1)-2*T(i)+T(i-1));
end
dy(n-1)=d*(TN-2*T(n-1)+T(n-2));


 

Post Wysłany: 26 Grudnia 2016, Pon 8:46 am Temat postu: Odpowiedz z cytatem
 
AUTOR:
alhambra
Może pisać książki


Dołączył: 05 Sie 2008
Posty: 675


Ogląda profil użytkownika Wyślij prywatną wiadomość
tak na szybko po obejrzeniu kodu ...

w linijce
[t,T] = ode45(@rownanie,linspace(0,200000,100),z,[],n);

usuń []
bo w sumie dajesz 4 parametry z czego jeden pusty w a funkcji
function dy = rownanie(t,T,n)
masz możliwość implementacji tylko 3 parametrów



_________________
Korepetycje z Matlaba, pomoc ....
 

Forum MATLAB Strona Główna-> MATLAB
Wyświetl posty z ostatnich:   

Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Wszystkie czasy w strefie CET (Europa)

Skocz do:  

Statystyki forum:



Od dnia 08.06.2006 forum odwiedzano 34068316
Najwięcej użytkowników 266 było obecnych 19 Lutego 2015, Czw 7:03 pm

Aktualnie online:




Najnowsze posty na forum:
(Scilab, Octave, ...) zmiana funkcji improfile  (23 Styczenia 2017, Pon 1:02 pm)
Zadanie : Obracanie obrazu  (23 Styczenia 2017, Pon 11:19 am)
Zadanie : Obracanie obrazu  (23 Styczenia 2017, Pon 11:18 am)
Zapis do pliku wektora po DCT obrazu  (22 Styczenia 2017, Nie 11:45 pm)
Zlece prosty projekt w Matlab  (22 Styczenia 2017, Nie 3:52 pm)
Wprowadzanie wlasnej daty  (22 Styczenia 2017, Nie 3:32 pm)
Zadanie  (22 Styczenia 2017, Nie 11:37 am)
Odwracanie macierzy, pomocy  (20 Styczenia 2017, Pią 9:20 pm)
[Matlab] odwracanie macierzy  (20 Styczenia 2017, Pią 9:17 pm)
Zamiana zmiennych globalnych na zmienne lokalne  (20 Styczenia 2017, Pią 4:57 pm)
Twoje prawa:
Nie możesz pisać nowych tematów
Nie możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach
Nie możesz załączać plików na tym forum
Nie możesz ściągać plików na tym forum