%% 给定 Pierce 参量
%% 运行行波管一维非线性计算(无切断)
clc;clear;
vc=2.998e08;
n=32; % macro particle number
% temp=1100; % Temperature of cathode(K)
%============================================
freq=31e09;
U0=13500;
I0=0.204;
rb=0.17e-3;
CC=0.05227;
bb=0.1832;
dd=0.01504;
qq=0.1892;
Pin=0.05; % ipput power
Fwin=sqrt(2*(1+CC*bb)^2/(CC*I0*U0)*Pin);
tspan=[0:0.1:10]; % theta
%=============================================
[vec_classic,vec_correct]=U0tovec(U0);
vec=vec_correct; % relativistic velocity
%==========================================
u0=linspace(-0*pi/2,4*pi/2-4*pi/2/n,n); % particle distribution
% v_p=pvdistri(I0,freq,n,u0,temp);
% u1=v_p(:,2)';
% du1=(v_p(:,1)'-mean(v_p(:,1)))/(vc);
% vec=vec+mean(v_p(:,1))/vc;
%==========================================
y0(1:2:2*n-1)=u0;
y0(2:2:2*n)=0;
y0(2*n+1)=Fwin;
y0(2*n+2)=-1i*y0(2*n+1)*(bb-1i*dd);
%%%%%%%%%%%%%%%%%%%%%%
% u0=linspace(-pi/2,3*pi/2,n);
% du0=0;
% y0(1:2:2*n-1)=u0;
% y0(2:2:2*n)=du0;
% y0(2*n+1)=Fwin;
% y0(2*n+2)=-i*y0(2*n+1)*(bb-i*dd);
%y0(2*n+2)=0
%=====================================================================
[t,y]=ode45(@(t,y)TWT1_d(t,y,CC,bb,dd,qq,freq,vec,rb,u0,n),tspan,y0);
%=====================================================================
beta_e=2*pi*freq/(vec*vc);
vec_ac=vec./(1+CC*y(:,2:2:2*n)); %电子的实际速度
vec_ac_mean=mean(vec_ac,2);%电子的平均速度
%[obsBeam,obsBeam2]=EBeam(vec,vec_ac,vec_ac_mean,I0,u0);
%obsE=Ecircuit(t,y(:,1:2:2*n-1),y(:,2*n+1),U0,I0,CC,beta_e,u0);
Pout=CC*I0*U0/(2)*(abs(y(:,end-1))).^2;
%===================================================================
figure(1),hold on;
for j=1:n
plot(t,y(:,2*j-1)/pi,'r');
end
figure(2),hold on;
plot(t,abs(y(:,end-1)))
figure(3),hold on;
plot(t,angle(y(:,end-1)))
figure(4),hold on;
plot(t,vec_ac,'.r')
plot(t, vec_ac_mean,'k')
figure(5),hold on;
%plot(t,10*log10(obsBeam),'r')
%plot(t,obsBeam2,'k')
%plot(t,10*log10(Pout), 'g')
figure(6),hold on;
[val,id]=max(abs(y(:,end-1)));
plot(y(id,1:2:2*n-1)/pi,vec_ac(id,:),'or')
figure(7),hold on;
plot(t,20*log10(abs(y(:,end-1))/Fwin),'or')
麻烦各位看看这个能否在mathcad里面实现以及怎么实现

没看明白，火星文么？

%这个是注释掉的在matlab里面可以跑起来上面就是出来的图片，现在的问题是这个程序如何在mathcad里面跑起来

还是没明白，先说整体思路，再逐条讲清细节要点。关键是你要告诉我这是在干什么，说得越清楚越好，说得越多越好。必须科普，不许讨价还价。

描述少于5000字，此问题不受理。

嗯 好的

刚进来还以为走错贴吧了。看来贴吧影响力越来越大了。

我弄到现在凑不齐5000字 这个难度可真大呀

我问你答吧。
（1）什么是 Pierce 参量？
（2） 什么是行波管？
（3） 什么是一维非线性计算？
（4） 什么无切断？
先回答这四个。

freq=31e09;
U0=13500;
I0=0.204;
rb=0.17e-3;
CC=0.05227;
bb=0.1832;
dd=0.01504;
qq=0.1892;
Pin=0.05; % ipput power
Fwin=sqrt(2*(1+CC*bb)^2/(CC*I0*U0)*Pin);
tspan=[0:0.1:10]; % theta
以上这些参数都是啥？一个一个说。他们之间有啥关系？列出计算式，手写列也可以，用mathtype打出来截图也可以。

还有，上面的那一堆参数的单位是什么？

后面的那一堆M78星文一样处理，每一行都要说，否则我的问题没完没了。