Ciao a tutti!
Mi scuso in anticipo per la domanda forse banale ma sono nuova nel mondo di Mathlab. Ho da risolvere un sistema di due equazioni in due incognite. Una equazione è alla quarta e inoltre ho sia le incognite che i termini noti che sono dei vettori di 24 componenti. Ho provato ad utilizzare vpasolve per ogni singola componente e ottengo dei valori ma quando voglio mettere il vettore insieme mi da come risultato empty sym. Ho provato anche con un ciclo for ma niente. Cosa sbaglio secondo voi? vi ringrazio in anticipo
il codice l'ho scritto così:
% Definizione dei coefficienti delle equazioni
sigma=5.76E-8;
epsilon=0.85;
a=0.65; %assorbimento
h=10.6; %coefficiente di convezione aria/estradosso
Rc4=0.596;%Resistenze
Rci=0.381;
data = readtable('Tabelladati.csv');
% Dividing training and test dataset
G= table2array(data(1:24,"G_invroma"));
Tint_inv= table2array(data(1:24,"T_int"));
Test_inv=table2array(data(1:24,"T_e"));
k1=a.*G;
k2=(epsilon*sigma).*Test_inv.^4;
k3=h.*Test_inv(1);
n=size(data);
z=omes(24,1);
k=k2+k3+k1;
a1=(-(sigma*epsilon).*z);
a2=(-(h+1/Rc4).*z);
a3=((1/Rc4).*z);
cond1=273.*z;
cond2=273.*z;
b2=((1/Rc4).*z);
b3=((-1/Rci-1/Rc4).*z);
b4=((1/Rci).*Tint_inv;
% Definizione delle variabili
syms ('x','real');
syms ('y','real');
eqn1=a1.*x^4+a2.*x+a3.*y==-k;
eqn2=b2.*x+b3.*y==-b4;
S = vpasolve([eqn1 eqn2],[x,y]);
% Visualizzazione delle soluzioni
disp('Le soluzioni del sistema di equazioni sono:');
disp(S);