Matlab Hogyan - Hidrodinamikai Rendszerek Tanszék

4. FEJEZET. NEMLINEÁRIS FELADATOK 20
ventillator.eps
4.3. ábra. Ventillátor- és fojtásjelleggörbe ill. a kiszámolt metszéspont.
1 function ventillator
ventillator.m
2
3 a=0.6;
4
5 q0=0.5; dp0=0.6;
6 [x,z] = fminsearch(@hiba,[q0,dp0],[],a);
7
8 xx=0:0.01:0.9;
9 yv=jelleggorbe(xx);
10 yf=1/a^2*xx.^2;
11 plot(xx,yv,xx,yf,x(1),x(2),’ro’)
12 axis([0 0.9 0 0.8]), grid on
13 xlabel(’Q’), ylabel(’\Delta p’)
14
15 function R=hiba(x,a)
16
17 Q =x(1); dp=x(2);
18 r1 = jelleggorbe(Q)-x(2);
19 r2 = Q-a*sqrt(dp);
20 R=r1^2+r2^2;
21
22 function dp=jelleggorbe(q)
23
24 dp=0.27+0.18*(1+1.5*(4*q-1)-0.5*(4*q-1).^3);
Külön magyarázatot igényel a 6. sorban az üres szögletes zárójel. Az a változó szeretnénk a főprogram 3.
sorában beállítani és aztán átadni az összes többi eljárásnak is. Ezt megtehetjük az fminsearch argumentumlistájában,
de ott a kezdeti értékek után az options változónak kell állnia és csak azután kerülhetnek
be az egyéb extra argumentumok. Ezért, bár semmilyen opciót nem kívánunk igénybe venni, ki kell tennünk
a szögletes zárójelet és csak azután kerülhet be az a változó.
Ártatlannak tűnhet még a 24. sor végén kitett pont a ∧3 előtt, azonban ez nagyon is fontos. Ugyanis ha a
jellegörbe eljárást csak skalár értékkel hívnánk meg, ez a pont nem kellene, de a 9. sorban az xx vektor