<SYSTEM_TASK:>
Given the following text description, write Python code to implement the functionality described below step by step
<END_TASK>
<USER_TASK:>
Description:
 Balanced Network De-embedding
Demonstration of balanced, i.e. 2N-port, network  de-embedding.
Setup
Step2:  Build fixture network
short length of mismatched line with connector-like input shunt capacitance
some crosstalk added with nudge
Step4:  Build DUT network
some length of mismatched lines
some crosstalk added with nudge
Step6:  Build the measurement
cascade Left, DUT and Right
add some noise
Step7:  Perform de-embedding
Step8:  Display results
<ASSISTANT_TASK:>
Python Code:
import numpy as np
import skrf as rf
from skrf.media import CPW
rf.stylely()
import matplotlib.pyplot as plt
# base parameters
freq = rf.Frequency(1e-3,10,1001,'ghz')
cpw  = CPW(freq, w=0.6e-3, s=0.25e-3, ep_r=10.6)
Explanation: Balanced Network De-embedding
Demonstration of balanced, i.e. 2N-port, network  de-embedding.
Setup
End of explanation
             l1
    0----+-=======-2
         |
         = c1
         |
        GND
        
             l1
    1----+-=======-3
         |
         = c1
         |
        GND
l1 = cpw.line(20, 'mm', z0=50, embed=True)
c1 = cpw.shunt_capacitor(C=0.15e-12, z0=50)
l1 = rf.connect(c1, 1, l1, 0)
li = rf.concat_ports([l1, l1], port_order='second')
Fix = li
Fix.name = 'Fix'
Fix.nudge(1e-4)
Left = Fix
# flip fixture for right side
Right = Fix.flipped()
Explanation: Build fixture network
short length of mismatched line with connector-like input shunt capacitance
some crosstalk added with nudge
End of explanation
        l2
    0-=======-2
        l2
    1-=======-3
l2 = cpw.line(50, 'mm', z0=50, embed=True)
DUT = rf.concat_ports([l2, l2], port_order='second')
DUT.name = 'DUT'
DUT.nudge(1e-5)
Explanation: Build DUT network
some length of mismatched lines
some crosstalk added with nudge
End of explanation
            Left        Meas         Right
             l1          l2           l1
    0----+-=======-2 0-=======-2 0-=======-+----2
         |                                 |
         = c1                              = c1
         |                                 |
        GND                               GND
        
             l1          l2           l1
    1----+-=======-3 1-=======-3 1-=======-+----3
         |                                 |
         = c1                              = c1
         |                                 |
        GND                               GND
Meas = Left ** DUT ** Right
Meas.name = 'Meas'
Meas.add_noise_polar(1e-5, 2)
Explanation: Build the measurement
cascade Left, DUT and Right
add some noise
End of explanation
DUTd = Left.inv ** Meas ** Right.inv
DUTd.name = 'DUTd'
Explanation: Perform de-embedding
End of explanation
fig, axarr = plt.subplots(2,2, sharex=True, figsize=(10,6))
ax = axarr[0,0]
Meas.plot_s_db(m=0, n=0, ax=ax)
DUTd.plot_s_db(m=0, n=0, ax=ax)
DUT.plot_s_db(m=0, n=0, ax=ax, ls=':', color='0.0')
ax.set_title('Return Loss')
ax.legend(loc='lower center', ncol=3)
ax.grid(True)
ax = axarr[0,1]
Meas.plot_s_db(m=2, n=0, ax=ax)
DUTd.plot_s_db(m=2, n=0, ax=ax)
DUT.plot_s_db(m=2, n=0, ax=ax, ls=':', color='0.0')
ax.set_title('Insertion Loss')
ax.legend(loc='lower center', ncol=3)
ax.grid(True)
ax = axarr[1,0]
Meas.plot_s_db(m=1, n=0, ax=ax)
DUTd.plot_s_db(m=1, n=0, ax=ax)
DUT.plot_s_db(m=1, n=0, ax=ax, ls=':', color='0.0')
ax.set_title('Isolation')
ax.legend(loc='lower center', ncol=3)
ax.grid(True)
ax = axarr[1,1]
Meas.plot_s_deg(m=2, n=0, ax=ax)
DUTd.plot_s_deg(m=2, n=0, ax=ax, marker='o', markevery=25)
DUT.plot_s_deg(m=2, n=0, ax=ax, ls=':', color='0.0')
ax.set_title('Insertion Loss')
ax.legend(loc='lower center', ncol=3)
ax.grid(True)
fig.tight_layout()
Explanation: Display results
End of explanation
<END_TASK>