SIM/qucs_s
1
0
Fork 0
mirror of https://github.com/ra3xdh/qucs_s synced 2025-03-28 21:13:26 +00:00
Code Issues Packages Projects Releases Wiki Activity

Add initial support for S-parameter analysis with Xyce

Browse Source
This commit is contained in:
Vadim Kuznetsov 2022-03-12 18:04:35 +01:00
parent da6bc383e1
commit cb2eb53621
6 changed files with 127 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

43
qucs/components/source_ac.cpp
View File

@ -16,7 +16,10 @@
***************************************************************************/
#include "source_ac.h"
#include "node.h"
#include "extsimkernels/spicecompat.h"
#include <cmath>
Source_ac::Source_ac()
{
@ -52,6 +55,7 @@ Source_ac::Source_ac()
tx = x1+4;
ty = y2+4;
Model = "Pac";
SpiceModel = "P";
Name = "P";
// This property must be the first one !
@ -86,3 +90,42 @@ Element* Source_ac::info(QString& Name, char* &BitmapFile, bool getNewOne)
if(getNewOne) return new Source_ac();
return 0;
}
QString Source_ac::ngspice_netlist()
{
QString s = "*Stub! Will be added in Ngspice-37\n";
return s;
}
QString Source_ac::xyce_netlist()
{
QString s = spicecompat::check_refdes(Name,SpiceModel);
foreach(Port *p1, Ports) {
QString nam = p1->Connection->Name;
if (nam=="gnd") nam = "0";
s += " "+ nam; // node names
}
s += QString(" port=%1 ").arg(getProperty("Num")->Value);
QString s_z0 = spicecompat::normalize_value(getProperty("Z")->Value);
double z0 = s_z0.toDouble();
QString s_p = spicecompat::normalize_value(getProperty("P")->Value);
double p = s_p.toDouble();
double vrms = sqrt(z0/1000.0)*pow(10, p/20.0);
double vamp = vrms*sqrt(2.0);
s += QString(" z0=%1 ").arg(s_z0);
QString f = spicecompat::normalize_value(getProperty("f")->Value);
s += QString(" AC %1 ").arg(vamp);
s += QString(" SIN 0 %1 %2").arg(vamp).arg(f);
s += "\n";
return s;
}
QString Source_ac::spice_netlist(bool isXyce)
{
if (isXyce) {
return xyce_netlist();
} else {
return ngspice_netlist();
}
}

5
qucs/components/source_ac.h
View File

@ -20,13 +20,16 @@
#include "component.h"
class Source_ac : public Component {
private:
QString ngspice_netlist();
QString xyce_netlist();
public:
Source_ac();
~Source_ac();
Component* newOne();
static Element* info(QString&, char* &, bool getNewOne=false);
QString spice_netlist(bool isXyce);
};
#endif

66
qucs/components/sp_sim.cpp
View File

@ -16,10 +16,13 @@
***************************************************************************/
#include "sp_sim.h"
#include "main.h"
#include "misc.h"
#include "schematic.h"
SP_Sim::SP_Sim()
{
isSimulation = true;
Description = QObject::tr("S parameter simulation");
QString s = Description;
@ -41,6 +44,7 @@ SP_Sim::SP_Sim()
ty = y2+1;
Model = ".SP";
Name = "SP";
SpiceModel = ".SP";
// The index of the first 4 properties must not changed. Used in recreate().
Props.append(new Property("Type", "lin", true,
@ -103,3 +107,65 @@ void SP_Sim::recreate(Schematic*)
Props.next()->Name = "Points";
}
}
int SP_Sim::getSPortsNumber()
{
int p_num = 0;
if (containingSchematic != NULL) {
auto comps = containingSchematic->DocComps;
for(Component *pc = comps.first(); pc != 0; pc = comps.next()) {
if (pc->Model == "Pac") p_num++;
}
return p_num;
} else {
return p_num;
}
}
QString SP_Sim::ngspice_netlist()
{
QString s = "*Stub!\n";
return s;
}
QString SP_Sim::xyce_netlist()
{
QString s = ".AC ";
QString unit;
if (Props.at(0)->Value=="log") { // convert points number for spice compatibility
double Np,Fstart,Fstop,fac = 1.0;
misc::str2num(Props.at(3)->Value,Np,unit,fac); // Points number
Np *= fac;
misc::str2num(Props.at(1)->Value,Fstart,unit,fac);
Fstart *= fac;
misc::str2num(Props.at(2)->Value,Fstop,unit,fac);
Fstop *= fac;
double Nd = ceil(log10(Fstop/Fstart)); // number of decades
double Npd = ceil(Np/Nd); // points per decade
s += QString("DEC %1 ").arg(Npd);
} else { // no need conversion
s += QString("LIN %1 ").arg(Props.at(3)->Value);
}
QString fstart = spicecompat::normalize_value(Props.at(1)->Value); // Start freq.
QString fstop = spicecompat::normalize_value(Props.at(2)->Value); // Stop freq.
s += QString("%1 %2 \n").arg(fstart).arg(fstop);
s += ".LIN format=touchstone sparcalc=1\n"; // enable s-param
int ports_num = getSPortsNumber();
s += ".PRINT ac format=std file=spice4qucs_sparam.prn ";
for (int i = 0; i < ports_num; i++) {
for (int j = 0; j < ports_num; j++) {
s += QString(" sdb(%1,%2) s(%1,%2) sp(%1,%2) y(%1,%2) z(%1,%2) ").arg(i+1).arg(j+1);
}
}
s += "\n";
return s;
}
QString SP_Sim::spice_netlist(bool isXyce)
{
if (isXyce) {
return xyce_netlist();
} else {
return ngspice_netlist();
}
}

5
qucs/components/sp_sim.h
View File

@ -22,12 +22,17 @@
class SP_Sim : public Component {
private:
QString xyce_netlist();
QString ngspice_netlist();
int getSPortsNumber();
public:
SP_Sim();
~SP_Sim();
Component* newOne();
static Element* info(QString&, char* &, bool getNewOne=false);
void recreate(Schematic*);
QString spice_netlist(bool isXyce);
};
#endif

5
qucs/extsimkernels/xyce.cpp
View File

@ -54,6 +54,7 @@ void Xyce::determineUsedSimulations(QStringList *sim_lst)
if (sim_typ==".NOISE") simulationsQueue.append("noise");
if (sim_typ==".TR") simulationsQueue.append("tran");
if (sim_typ==".HB") simulationsQueue.append("hb");
if (sim_typ==".SP") simulationsQueue.append("sp");
if (sim_typ==".SENS_XYCE") simulationsQueue.append("sens");
if (sim_typ==".SENS_TR_XYCE") simulationsQueue.append("sens_tr");
if (sim_typ==".XYCESCR") simulationsQueue.append(pc->Name); // May be >= XYCE scripts
@ -163,6 +164,7 @@ void Xyce::createNetlist(QTextStream &stream, int , QStringList &simulations,
if ((sim_typ==".NOISE")&&(sim=="noise")) stream<<s;
if ((sim_typ==".SENS_XYCE")&&(sim=="sens")) stream<<s;
if ((sim_typ==".SENS_TR_XYCE")&&(sim=="sens_tr")) stream<<s;
if ((sim_typ==".SP")&&(sim=="sp")) stream<<s;
if (sim==pc->Name) stream<<s; // Xyce scripts
if ((sim_typ==".TR")&&(sim=="tran")){
stream<<s;
@ -239,6 +241,9 @@ void Xyce::createNetlist(QTextStream &stream, int , QStringList &simulations,
write_str.clear();
outputs.append("spice4qucs.sens_tr.cir.SENS.prn");
outputs.append("spice4qucs.sens_tr.cir.TRADJ.prn");
} else if (sim=="sp") {
write_str.clear();
outputs.append("spice4qucs_sparam.prn");
} else {
write_str = QString(".PRINT %1 format=raw file=%2 %3\n").arg(sim).arg(filename).arg(nods);
outputs.append(filename);

8
qucs/module.cpp
View File

@ -294,8 +294,8 @@ void Module::registerModules (void) {
REGISTER_SOURCE_1 (Ampere_dc);
REGISTER_SOURCE_1 (Volt_ac);
REGISTER_SOURCE_1 (Ampere_ac);
REGISTER_SOURCE_1 (Source_ac);
if (QucsSettings.DefaultSimulator == spicecompat::simQucsator) {
REGISTER_SOURCE_1 (Source_ac);
REGISTER_SOURCE_1 (Volt_noise);
REGISTER_SOURCE_1 (Ampere_noise);
}
@ -506,12 +506,12 @@ void Module::registerModules (void) {
REGISTER_SIMULATION_1 (DC_Sim);
REGISTER_SIMULATION_1 (TR_Sim);
REGISTER_SIMULATION_1 (AC_Sim);
if (QucsSettings.DefaultSimulator == spicecompat::simQucsator)
REGISTER_SIMULATION_1 (SP_Sim);
if (QucsSettings.DefaultSimulator == spicecompat::simQucsator||
QucsSettings.DefaultSimulator == spicecompat::simXyceSer||
QucsSettings.DefaultSimulator == spicecompat::simXycePar)
QucsSettings.DefaultSimulator == spicecompat::simXycePar) {
REGISTER_SIMULATION_1 (SP_Sim);
REGISTER_SIMULATION_1 (HB_Sim);
}
REGISTER_SIMULATION_1 (Param_Sweep);
REGISTER_SIMULATION_1 (Digi_Sim);
if (QucsSettings.DefaultSimulator == spicecompat::simQucsator)