Fix docstring examples for better mkdocstrings rendering

Remove doctest format (>>> and ...) from Example sections
to improve rendering in generated API documentation.

Co-Authored-By: Claude Opus 4.5 <noreply@anthropic.com>

Author: jman4162

src/phased_array_systems/models/digital/bandwidth.py

@@ -35,8 +35,7 @@ def beam_bandwidth_product(
         Total beam-bandwidth product in Hz
 
     Example:
-        >>> beam_bandwidth_product(4, 100e6)
-        400000000.0
+        beam_bandwidth_product(4, 100e6)  # Returns 400e6
     """
     return n_beams * bandwidth_per_beam_hz
 
@@ -60,8 +59,7 @@ def max_simultaneous_beams(
         Maximum number of simultaneous beams (integer)
 
     Example:
-        >>> max_simultaneous_beams(1e9, 100e6)
-        9
+        max_simultaneous_beams(1e9, 100e6)  # Returns 9
     """
     effective_bw = processing_bandwidth_hz / overhead_factor
     return int(effective_bw / bandwidth_per_beam_hz)
@@ -93,8 +91,8 @@ def digital_beamformer_data_rate(
             - per_element_gbps: Rate per element
 
     Example:
-        >>> result = digital_beamformer_data_rate(256, 1e9, 14)
-        >>> print(f"Total data rate: {result['with_overhead_gbps']:.1f} Gbps")
+        result = digital_beamformer_data_rate(256, 1e9, 14)
+        print(f"Total data rate: {result['with_overhead_gbps']:.1f} Gbps")
     """
     raw_rate_bps = n_elements * sample_rate_hz * bits_per_sample * n_channels
     raw_rate_gbps = raw_rate_bps / 1e9
@@ -137,8 +135,8 @@ def channelizer_output_rate(
             - samples_per_channel_per_sec: Output samples per channel
 
     Example:
-        >>> result = channelizer_output_rate(1e9, 16)
-        >>> print(f"Channel BW: {result['channel_bandwidth_hz']/1e6:.1f} MHz")
+        result = channelizer_output_rate(1e9, 16)
+        print(f"Channel BW: {result['channel_bandwidth_hz']/1e6:.1f} MHz")
     """
     channel_bandwidth_hz = input_bandwidth_hz / n_channels
     channel_sample_rate_hz = channel_bandwidth_hz * overlap_factor
@@ -177,8 +175,8 @@ def processing_margin(
             - headroom_percent: Remaining capacity
 
     Example:
-        >>> result = processing_margin(1000, 600)
-        >>> print(f"Utilization: {result['utilization_percent']:.1f}%")
+        result = processing_margin(1000, 600)
+        print(f"Utilization: {result['utilization_percent']:.1f}%")
     """
     margin_ratio = available_throughput_gops / required_throughput_gops
     margin_db = 10 * math.log10(margin_ratio) if margin_ratio > 0 else float('-inf')
@@ -220,8 +218,8 @@ def beamformer_operations(
             - method: 'time_domain' or 'frequency_domain'
 
     Example:
-        >>> result = beamformer_operations(256, 4, 100e6)
-        >>> print(f"Required: {result['total_gops']:.1f} GOPS")
+        result = beamformer_operations(256, 4, 100e6)
+        print(f"Required: {result['total_gops']:.1f} GOPS")
     """
     if fft_size > 0:
         # Frequency-domain: FFT + multiply + IFFT

src/phased_array_systems/models/digital/converters.py

@@ -33,10 +33,8 @@ def enob_to_snr(enob: float) -> float:
         SNR in dB
 
     Example:
-        >>> enob_to_snr(12)
-        74.0
-        >>> enob_to_snr(14)
-        86.04
+        enob_to_snr(12)   # Returns 74.0
+        enob_to_snr(14)   # Returns 86.04
     """
     return 6.02 * enob + 1.76
 
@@ -56,8 +54,7 @@ def snr_to_enob(snr_db: float) -> float:
         Effective number of bits
 
     Example:
-        >>> snr_to_enob(74.0)
-        12.0
+        snr_to_enob(74.0)  # Returns 12.0
     """
     return (snr_db - 1.76) / 6.02
 
@@ -76,10 +73,8 @@ def enob_to_sfdr(enob: float, margin_db: float = 0.0) -> float:
         Estimated SFDR in dB
 
     Example:
-        >>> enob_to_sfdr(12)
-        74.0
-        >>> enob_to_sfdr(12, margin_db=6)  # Conservative estimate
-        68.0
+        enob_to_sfdr(12)                   # Returns 74.0
+        enob_to_sfdr(12, margin_db=6)      # Returns 68.0 (conservative)
     """
     return enob_to_snr(enob) - margin_db
 
@@ -119,8 +114,7 @@ def quantization_noise_floor(
         Noise floor spectral density in dBm/Hz
 
     Example:
-        >>> quantization_noise_floor(12, 0, 100e6, 250e6)
-        -154.0  # Approximate
+        quantization_noise_floor(12, 0, 100e6, 250e6)  # Returns ~-154.0
     """
     snr_db = enob_to_snr(enob)
     nyquist_bw = sample_rate_hz / 2
@@ -148,8 +142,7 @@ def sample_rate_for_bandwidth(
         Required sample rate in Hz
 
     Example:
-        >>> sample_rate_for_bandwidth(100e6)
-        250000000.0
+        sample_rate_for_bandwidth(100e6)  # Returns 250e6
     """
     return signal_bandwidth_hz * oversampling_ratio
 
@@ -170,8 +163,7 @@ def max_signal_bandwidth(
         Maximum signal bandwidth in Hz
 
     Example:
-        >>> max_signal_bandwidth(1e9)
-        400000000.0
+        max_signal_bandwidth(1e9)  # Returns 400e6
     """
     return sample_rate_hz / oversampling_ratio
 
@@ -201,8 +193,8 @@ def adc_dynamic_range(
             - dynamic_range_db: Usable dynamic range
 
     Example:
-        >>> result = adc_dynamic_range(14, noise_figure_db=3, bandwidth_hz=100e6)
-        >>> print(f"Dynamic range: {result['dynamic_range_db']:.1f} dB")
+        result = adc_dynamic_range(14, noise_figure_db=3, bandwidth_hz=100e6)
+        print(f"Dynamic range: {result['dynamic_range_db']:.1f} dB")
     """
     snr_db = enob_to_snr(enob)
 
@@ -255,8 +247,8 @@ def dac_output_power(
             - noise_floor_dbm: Quantization noise floor
 
     Example:
-        >>> result = dac_output_power(14, full_scale_dbm=10, backoff_db=6)
-        >>> print(f"Operating power: {result['operating_power_dbm']:.1f} dBm")
+        result = dac_output_power(14, full_scale_dbm=10, backoff_db=6)
+        print(f"Operating power: {result['operating_power_dbm']:.1f} dBm")
     """
     snr_db = enob_to_snr(enob)
     sfdr_db = enob_to_sfdr(enob)

src/phased_array_systems/models/digital/scheduling.py

@@ -117,9 +117,9 @@ def timeline_utilization(timeline: Timeline) -> dict[str, float]:
             - by_function_percent: Dict of function -> percentage
 
     Example:
-        >>> tl = Timeline(dwells=[...], frame_time_ms=100)
-        >>> util = timeline_utilization(tl)
-        >>> print(f"Utilization: {util['total_utilization']*100:.1f}%")
+        tl = Timeline(dwells=[...], frame_time_ms=100)
+        util = timeline_utilization(tl)
+        print(f"Utilization: {util['total_utilization']*100:.1f}%")
     """
     total_dwell_ms = timeline.total_dwell_time_ms
     idle_time_ms = max(0, timeline.frame_time_ms - total_dwell_ms)
@@ -169,13 +169,13 @@ def max_update_rate(
             - scan_time_s: Time for one complete scan
 
     Example:
-        >>> # Search ±60° az, ±30° el with 3° beam
-        >>> result = max_update_rate(
-        ...     scan_volume_sr=2.0,  # ~hemisphere
-        ...     beam_solid_angle_sr=0.003,  # ~3° beam
-        ...     dwell_time_us=100
-        ... )
-        >>> print(f"Update rate: {result['update_rate_hz']:.2f} Hz")
+        # Search ±60° az, ±30° el with 3° beam
+        result = max_update_rate(
+            scan_volume_sr=2.0,       # ~hemisphere
+            beam_solid_angle_sr=0.003,  # ~3° beam
+            dwell_time_us=100,
+        )
+        print(f"Update rate: {result['update_rate_hz']:.2f} Hz")
     """
     n_beam_positions = math.ceil(scan_volume_sr / beam_solid_angle_sr)
     time_per_position_us = dwell_time_us + overhead_us
@@ -218,14 +218,14 @@ def search_timeline(
         Timeline with search dwells
 
     Example:
-        >>> tl = search_timeline(
-        ...     azimuth_range_deg=(-60, 60),
-        ...     elevation_range_deg=(0, 30),
-        ...     azimuth_step_deg=3.0,
-        ...     elevation_step_deg=3.0,
-        ...     dwell_time_us=100
-        ... )
-        >>> print(f"Generated {tl.n_dwells} dwells")
+        tl = search_timeline(
+            azimuth_range_deg=(-60, 60),
+            elevation_range_deg=(0, 30),
+            azimuth_step_deg=3.0,
+            elevation_step_deg=3.0,
+            dwell_time_us=100,
+        )
+        print(f"Generated {tl.n_dwells} dwells")
     """
     dwells = []
 
@@ -283,17 +283,17 @@ def interleaved_timeline(
         Timeline with interleaved function dwells
 
     Example:
-        >>> tl = interleaved_timeline(
-        ...     functions=[
-        ...         {"function": Function.RADAR_SEARCH, "time_percent": 60,
-        ...          "dwell_time_us": 100, "dwells_per_burst": 10},
-        ...         {"function": Function.RADAR_TRACK, "time_percent": 30,
-        ...          "dwell_time_us": 50, "dwells_per_burst": 5},
-        ...         {"function": Function.ESM, "time_percent": 10,
-        ...          "dwell_time_us": 200, "dwells_per_burst": 2},
-        ...     ],
-        ...     frame_time_ms=100
-        ... )
+        tl = interleaved_timeline(
+            functions=[
+                {"function": Function.RADAR_SEARCH, "time_percent": 60,
+                 "dwell_time_us": 100, "dwells_per_burst": 10},
+                {"function": Function.RADAR_TRACK, "time_percent": 30,
+                 "dwell_time_us": 50, "dwells_per_burst": 5},
+                {"function": Function.ESM, "time_percent": 10,
+                 "dwell_time_us": 200, "dwells_per_burst": 2},
+            ],
+            frame_time_ms=100,
+        )
     """
     dwells = []
 

src/phased_array_systems/models/rf/cascade.py

@@ -43,8 +43,7 @@ def noise_figure_to_temp(nf_db: float, t0: float = T0) -> float:
         Equivalent noise temperature in Kelvin
 
     Example:
-        >>> noise_figure_to_temp(3.0)  # 3 dB NF
-        288.6  # Approximately
+        noise_figure_to_temp(3.0)  # 3 dB NF -> ~288.6 K
     """
     f_linear = 10 ** (nf_db / 10)
     return t0 * (f_linear - 1)
@@ -63,8 +62,7 @@ def noise_temp_to_figure(te: float, t0: float = T0) -> float:
         Noise figure in dB
 
     Example:
-        >>> noise_temp_to_figure(290)
-        3.01  # Approximately
+        noise_temp_to_figure(290)  # 290 K -> ~3.01 dB
     """
     f_linear = 1 + te / t0
     return 10 * math.log10(f_linear)
@@ -93,12 +91,12 @@ def friis_noise_figure(
             - stage_contributions_db: NF contribution from each stage
 
     Example:
-        >>> # LNA (15 dB gain, 2 dB NF) -> Mixer (10 dB loss, 10 dB NF)
-        >>> result = friis_noise_figure([
-        ...     (15, 2),    # LNA
-        ...     (-10, 10),  # Mixer (loss = negative gain)
-        ... ])
-        >>> print(f"System NF: {result['total_nf_db']:.2f} dB")
+        # LNA (15 dB gain, 2 dB NF) -> Mixer (10 dB loss, 10 dB NF)
+        result = friis_noise_figure([
+            (15, 2),    # LNA
+            (-10, 10),  # Mixer (loss = negative gain)
+        ])
+        print(f"System NF: {result['total_nf_db']:.2f} dB")
     """
     if not stages:
         return {
@@ -166,11 +164,11 @@ def system_noise_temperature(
             - system_nf_db: Effective system noise figure
 
     Example:
-        >>> result = system_noise_temperature(
-        ...     antenna_temp_k=50,  # Cold sky
-        ...     receiver_nf_db=2.0,
-        ...     line_loss_db=0.5
-        ... )
+        result = system_noise_temperature(
+            antenna_temp_k=50,      # Cold sky
+            receiver_nf_db=2.0,
+            line_loss_db=0.5,
+        )
     """
     # Receiver noise temperature
     t_rx = noise_figure_to_temp(receiver_nf_db)
@@ -210,8 +208,7 @@ def cascade_gain(gains_db: list[float]) -> float:
         Total gain in dB
 
     Example:
-        >>> cascade_gain([20, -3, 15, -6])  # LNA, filter, amp, cable
-        26
+        cascade_gain([20, -3, 15, -6])  # LNA, filter, amp, cable -> 26 dB
     """
     return sum(gains_db)
 
@@ -248,11 +245,11 @@ def cascade_iip3(
             - total_gain_db: Cascaded gain
 
     Example:
-        >>> # LNA (15dB, +5dBm IIP3) -> Mixer (-10dB, +10dBm IIP3)
-        >>> result = cascade_iip3([
-        ...     (15, 5),
-        ...     (-10, 10),
-        ... ])
+        # LNA (15dB, +5dBm IIP3) -> Mixer (-10dB, +10dBm IIP3)
+        result = cascade_iip3([
+            (15, 5),
+            (-10, 10),
+        ])
     """
     if not stages:
         return {"iip3_dbm": float('inf'), "oip3_dbm": float('inf'), "total_gain_db": 0}
@@ -325,11 +322,11 @@ def sfdr_from_iip3(
             - max_signal_dbm: Maximum signal before spurs exceed noise
 
     Example:
-        >>> result = sfdr_from_iip3(
-        ...     iip3_dbm=5,
-        ...     noise_floor_dbm_hz=-170,
-        ...     bandwidth_hz=1e6
-        ... )
+        result = sfdr_from_iip3(
+            iip3_dbm=5,
+            noise_floor_dbm_hz=-170,
+            bandwidth_hz=1e6,
+        )
     """
     noise_floor_dbm = noise_floor_dbm_hz + 10 * math.log10(bandwidth_hz)
     sfdr_db = (2 / 3) * (iip3_dbm - noise_floor_dbm)
@@ -387,11 +384,11 @@ def mds_from_noise_figure(
             - ktb_dbm: Thermal noise power
 
     Example:
-        >>> result = mds_from_noise_figure(
-        ...     noise_figure_db=3,
-        ...     bandwidth_hz=1e6,
-        ...     snr_required_db=10
-        ... )
+        result = mds_from_noise_figure(
+            noise_figure_db=3,
+            bandwidth_hz=1e6,
+            snr_required_db=10,
+        )
     """
     # kT in dBm/Hz at T0
     kt_dbm_hz = 10 * math.log10(K_B * t0 * 1000)  # *1000 for mW
@@ -471,15 +468,15 @@ def cascade_analysis(
             - stage_names: Names of each stage
 
     Example:
-        >>> stages = [
-        ...     RFStage("LNA", gain_db=20, noise_figure_db=1.5, iip3_dbm=-5),
-        ...     RFStage("Filter", gain_db=-2, noise_figure_db=2, iip3_dbm=30),
-        ...     RFStage("Mixer", gain_db=-8, noise_figure_db=8, iip3_dbm=15),
-        ...     RFStage("IF Amp", gain_db=30, noise_figure_db=4, iip3_dbm=10),
-        ... ]
-        >>> result = cascade_analysis(stages, bandwidth_hz=10e6)
-        >>> print(f"System NF: {result['total_nf_db']:.2f} dB")
-        >>> print(f"SFDR: {result['sfdr_db']:.1f} dB")
+        stages = [
+            RFStage("LNA", gain_db=20, noise_figure_db=1.5, iip3_dbm=-5),
+            RFStage("Filter", gain_db=-2, noise_figure_db=2, iip3_dbm=30),
+            RFStage("Mixer", gain_db=-8, noise_figure_db=8, iip3_dbm=15),
+            RFStage("IF Amp", gain_db=30, noise_figure_db=4, iip3_dbm=10),
+        ]
+        result = cascade_analysis(stages, bandwidth_hz=10e6)
+        print(f"System NF: {result['total_nf_db']:.2f} dB")
+        print(f"SFDR: {result['sfdr_db']:.1f} dB")
     """
     if not stages:
         return {}