Problema funzione di trasferimento MatLab/Simulink

[Tidus]

Salve a tutti!
Ecco il mio problema: sto analizzando un modello elettrico di un'interfaccia microelettrodo-cellula e, per prima cosa, dovrei studiarne il comportamento in frequenza, tramite un classico diagramma di Bode.
Il circuito è questo:

Sorgente FidoCadJ | Ingrandisci | Cos'è Fidocad

dove Z è un'impedenza complessa, la cui formula è:
.
Devo ammettere che la Z, così scritta, con quell'esponente, mi sembra alquanto strana; nell''articolo accademico dal quale sto prendendo questo modello, gli autori sostengono che un' impedenza ad angolo di fase costante (piuttosto che un condensatore) risulta in miglior accordo con le misurazioni ottenute ed è un fattore empirico relativo alle irregolarità della superficie.

Ho allora inizialmente provato a disegnare il circuito su Simulink (primo utilizzo per me), ma mi sono bloccato quasi subito, perché ho notato che nella libreria Simscape non è presente il componente impedenza, il quale compare invece nella libreria SimRF. L'ho comunque selezionato, ma, quando faccio partire la simulazione del circuito, ottengo questo messaggio di errore: "Transfer function is not allowed in core Simscape.". Non ho trovato ancora nessuna spiegazione del perché ottenga questo errore, forse potrebbe essere dovuto al fatto che non stia usando tutti componenti Simscape e al software questo non piaccia.

Qualcuno potrebbe aiutarmi su questo fronte?

Ad ogni modo, ho allora provato ad agire sul fronte MatLab: nello stesso paper, è scritta la funzione di trasferimento:

dove:



Ma anche qui trovo degli ostacoli: documentandomi sulle varie guide, ho notato che in MatLab, il comando Bode(sys) ha bisogno di una funzione di trasferimento "canonica", del tipo .
C'è un altro modo per ottenere un diagramma di Bode, partendo da una funzione di trasferimento "non canonica" come la mia?
Ho anche pensato di fare la sostituzione e poi carta e penna sostituire e in in modo tale da ottenere una funzione di trasferimento nel formato che piace al comando Bode(sys), ma, a parte il casino che ne uscirebbe fuori, il problema principale sarebbe di nuovo quel fattore che non saprei come maneggiare.

C'è qualcuno che mi saprebbe dare una mano? Mi piacerebbe di più risolvere con Simulink, perché poi mi piacerebbe direttamente lavorare con il circuito e fare delle modifiche per vedere come varia il comportamento in frequenza, ma visto come stanno le cose, mi basterebbe anche sistemare dal lato MatLab!

Grazie mille!

[g.schgor]

Se mi dai i valori, te lo simulo con Microcap.

[DanteCpp]

Qua c'è una panoramica.

A me piace definire prima la variabile complessa

Codice: Seleziona tutto

s=tf([1 0],[1])

e poi costruire le funzioni di trasferimento più complesse come composizione di questa e funzioni elementari.

[MarkyMark]

Come può essere un'impedenza se ? Non mi tornano le dimensioni

[SandroCalligaro]

A me piace definire prima la variabile complessa

Codice: Seleziona tutto

s=tf([1 0],[1])

La si può definire anche in modo più "comodo" come

Codice: Seleziona tutto

s=tf('s')

[IsidoroKZ]

nell''articolo accademico dal quale sto prendendo questo modello, gli autori sostengono ...

Riferimento bibliografico all'articolo?

[dimaios]

Le fractional order transfer functions non sono simulabili nativamente da Matlab ma esistono varie soluzioni per farlo con risorse esterne.
Te ne propongo una, vedi se riesci ad utilizzarla con successo nella tua applicazione.

FOMCON Toolbox

[Tidus]

Se mi dai i valori, te lo simulo con Microcap.

Grazie per la disponibilità, ma mi piacerebbe avere qualcosa sottomano con cui "giocare", piuttosto che avere un semplice risultato. Non conscendo Microcap, mi sono informato e ho visto che è un software a pagamento, ma che esiste anche la versione di prova. Sapresti dirmi se questa DEMO è sufficiente per i miei scopi?

Qua c'è una panoramica.

A me piace definire prima la variabile complessa

Codice: Seleziona tutto

s=tf([1 0],[1])

e poi costruire le funzioni di trasferimento più complesse come composizione di questa e funzioni elementari.

Grazie mille, ora provo a vedere se così facendo riesco a risolvere il problema.

Come può essere un'impedenza se ? Non mi tornano le dimensioni

Già, non tornano i conti nemmeno a me!

Riferimento bibliografico all'articolo?

Hai ragione, mi sono dimenticato di inserirlo. L'articolo è questo qui e ciò di cui vi sto parlando riguarda il capitolo 2: ARTICOLO

Le fractional order transfer functions non sono simulabili nativamente da Matlab ma esistono varie soluzioni per farlo con risorse esterne.
Te ne propongo una, vedi se riesci ad utilizzarla con successo nella tua applicazione.

FOMCON Toolbox

Grazie mille, non sapevo dell'esistenza di questo Tool! Provo a vedere un po' e ti faccio sapere.

[Tidus]

Qua c'è una panoramica.

A me piace definire prima la variabile complessa

Codice: Seleziona tutto

s=tf([1 0],[1])

e poi costruire le funzioni di trasferimento più complesse come composizione di questa e funzioni elementari.

Niente da fare: il fatto che debba fare un (all'interno di ), con non scalare intero, fa bloccare MatLab, dandomi appunto questo errore: "In the expression "M^K", the exponent K must be a scalar integer."
Grazie comunque per l'idea!

[g.schgor]

Sì, la versione Demo (gratuita) di Microcap
è normalmente sufficiiccper ente per
la simulazione di piccoli circuiti.
Vantaggi sono la semplicità d'uso e
l'immediatezza dei risultati.