您现在的位置:首页 >> ★免费资源 >> 源码下载 >> 内容

matlab代做highspeedlogic★基于matlab编程微波系统

时间:2015-10-4 18:53:42 点击:

  核心提示:matlab代做,FPGA代做,淘宝,专业代做MATLAB、FPGA博士/硕士/本科毕业设计、项目仿真、Coursework、Assignment ...

matlab代做FPGA代做淘宝专业代做MATLABFPGA博士/硕士/本科毕业设计项目仿真CourseworkAssignment

QQ: 1224848052

function scArc(Point, LinCol)
 %scArc : Draws an arc (a circle) on the smith chart passing through a given point
%
%  SYNOPSIS:
%     Draws an arc (a circle) on the smith chart passing through a given point.
%     the point in question is defined as an (r,x) pair. This function should only
%     be called when a smith chart figure is already present.
%
%     See also scDraw, scInv, scMove, scMatchCirc
%     
%  SYNTAX:
%     scArc(Point, LinCol)
%
%  INPUT ARGUMENTS:
%     Point  : values of r and x defining a point on the smith chart as a vector [r x]
%     LinCol : color of the arc (or circle)
%
%  OUTPUT ARGUMENT:
%     none
%
%  EXAMPLE:
%         The Command sequence 
%         scDraw;
%         scArc([2 3]);
%         will draw a blank smith chart and draw a circle having center at the
%         smith chart center [1 0] and passing through the point [2 3]
%         in the smith chart corresponding to an impedance (2+j3)*Z_L. This
%         may be thought of as the locus of impedance along the lossless
%         line of line impedance Z_L.
 
%
%     Mohammad Ashfaq - (31-05-2000)
%     Mohammad Ashfaq - (13-04-2006) Modified (example included)
%

 if nargin == 1
     LinCol = 'm';
 end
     
 r1 = Point(1);
 x1 = Point(2);

 %x2 = to(2);
 %r2 = to(1);

 [u1, v1] = scPOI(r1, x1);
 %[u2,v2]=point_of_intersection(r2,x2);
 x = linspace(-sqrt(u1^2+v1^2),sqrt(u1^2+v1^2),500);
 
 plot(x,sqrt(u1^2+v1^2-x.^2),LinCol);
 plot(x,-sqrt(u1^2+v1^2-x.^2),LinCol);

 % MARK POINT
 plot(u1,v1,'r*')

 function scConCirc(r,LinCol)
 %scConCirc: draws concentric circles about the origin
%
%  SYNOPSIS:
%     Draws circle about origin in desired color.
%
%     See also scDraw, scInv, scMove, scMatchCirc  
%     
%  SYNTAX:
%     scConCirc(r,color)
%
%  INPUT ARGUMENTS:
%     r      : radius of the circle
%     LinCol : desired color of the arc, optional, default = 'm', magenta
%
%  OUTPUT ARGUMENT:
%     none
%
%
%  EXAMPLE:
%         The Command sequence 
%         scDraw;
%         scConCirc(0.4);
%         will draw a blank smith chart and draw a circle having center at the
%         smith chart center [1 0] and a radius 0.4. This may be thought of as
%         the locus of all the points having a reflection factor 0.4.
%
%     Mohammad Ashfaq - (31-05-2000)
%     Mohammad Ashfaq - (13-04-2006) Modified (example included)
%
 
 if nargin == 1
    LinCol = 'm';
 end

 x = linspace(-r,r,200);
 plot(x, sqrt(r^2-x.^2),LinCol);
 plot(x,-sqrt(r^2-x.^2),LinCol);
 

%scDraw : Draws a blank smith chart
%
%  SYNOPSIS:
%     When called without arguments, the function draws a standard smith chart.
%     Otherwise, one may want to draw only specific curves for r and/or x values.
%     For different calling syntaxes see under SYNTAX.
%
%     See also scRay, scMove, scConCirc, scMatchCirc : available for different tasks
%     
%     
%  SYNTAX:
%     [h] = scDraw(r, x, ChColor)
%           draws only the r and x lines required
%     [h] = scDraw
%           draws a standard smith chart
%
%     And one can suppress the secondary details like scales and labels by calling:
%     [h] = scDraw(0)
%           draws 'only' a standard smith chart without r and m scales etc.
%
%  INPUT ARGUMENTS:
%     r       : A vector consisting of the desired values of r 
%     x       : A vector consisting of the desired values of x
%     ChColor : Color of the smith chart
%     xL      : a vector of x values containing left termination points for r circles
%     xR      : a vector of x values containing right termination points for r circles
%     rL      : a vector of r values containing left termination points for x arcs
%     rR      : a vector of r values containing right termination points for r arcs
%
%  OUTPUT ARGUMENT:
%     h : figure handle. If no output argument is given, the handle is returned in
%         the workspace variable ans. 
%
%  EXAMPLES:
%     (1) scDraw;
%         will draw a blank smith chart
%
%     (2) h = scDraw;
%         will draw a blank smith chart and return the gigure handle h.
%
%     (3) scDraw(0);
%         will draw a blank smith chart without scales etc.
%
%     (4) scDraw([1 2 3 4 5],[.1 2 3 4]);
%         will draw only the r lines specified by [1 2 3 4 5] and x lines [.1 2 3 4]
%
%     Mohammad Ashfaq - (31-05-2000)
%     Mohammad Ashfaq - (13-04-2006) Modified (examples included)

function h = scDraw(r, x, ChColor, xR, xL, rR, rL)

FullMapWithLabels = 1;
if nargin==1
    FullMapWithLabels   = 0;
end
if nargin<=1
    % DEFINE STANDARD DEFAULT FOR r
    r= [0.00, 0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.40, 0.45, 0.50, 0.50, 0.55, 0.60, 0.65,...
            0.70, 0.75, 0.80, 0.85, 0.90, 0.95, 1.00, 1.10, 1.20, 1.30, 1.40, 1.50, 1.50, 1.60, 1.70,...
            1.80, 1.90, 2.00, 2.20, 2.40, 2.50, 2.60, 2.80, 3.00, 3.50, 4.00, 4.50, 5.00, 6.00, 7.00,...
            8.00, 9.00, 10.0, 15.0, 20.0, 50.0];
    xR=[   0,    1,    2,    1,    4,    1,    2,    1,    4,    1,    2,    5,     1,   4,    1,...
            2,    1,    4,    1,    2,    1,   10,    2,    4,    2,    4,    2,     5,   4,    2,...
            4,    2,   20,    2,    2,    5,    2,    2,   10,    5,   10,    5,    10,  10,   10,...
            10,   10,   20,   20,   50,   0];
    
    xL=[   0,    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,     4,     0,   0,    0,...
            0,    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,     0,     4,   0,    0,...
            0,    0,    0,    0,    0,    2,    0,    0,    0,    0,    0,     0,     0,   0,    0,...
            0,    0,    0,    0,    0,    0];
    rPrint = [0.1:0.1:1.0,1.5 2.0 3.0 4.0 5.0 10 20 50];
    
    % DEFINE STANDARD DEFAULT FOR x
    x=[0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.40, 0.45, 0.50,  0.50, 0.55, 0.60, 0.65, 0.70,...
            0.75, 0.80, 0.85, 0.90, 0.95, 1.00, 1.10, 1.20, 1.30, 1.40,  1.50, 1.50, 1.60, 1.70, 1.80,...
            1.90, 2.00, 2.20, 2.40, 2.50, 2.60, 2.80, 3.00, 3.20, 3.40,  3.50, 3.60, 3.80, 4.00, 4.50,...
            5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 15.0, 20.0, 50.0];
    rR=  [1,    2,    1,    3,    1,    2,    1,    3,    1,    2,     5,    1,    3,    1,    2,...
            1,    3,    1,    2,    1,   10,    2,    3,    2,    3,     2,    5,    3,    2,    3,...
            2,   10,    2,    2,    5,    2,    2,   10,    2,    2,     5,    2,    2,   10,    5,...
            10,   10,   10,   10,   10,   20,   20,   50,   0];
    
    rL=  [0,    0,    0,    0,    0,    0,    0,    0,    0,    0,     3,    0,     0,   0,    0,...
            0,    0,    0,    0,    0,    0,    0,    0,    0,    0,     0,    3,     0,   0,    0,...
            0,    0,    0,    0,    2,    0,    0,    0,    0,    0,     2,    0,     0,   0,    0,...
            0,    0,    0,    0,    0,    0,    0     0,    0];
    xPrint  = [0.1:0.1:1.0,1.5 2.0 3.0 4.0 5.0 10 20 50];
    ChColor = 'b';
end

if (nargin == 2)||(nargin ==3)
    rR = zeros(size(x));
    rL = zeros(size(x));
    xR = zeros(size(r));
    xL = zeros(size(r));
    rPrint = r;
    xPrint = x;
end
if nargin ==2 
    ChColor = 'b';
end

if isempty(get(gcf, 'name')) || nargin ==0
    h = gcf;
    set(gcf,'MenuBar','none','numbertitle','off','name','Smith Chart');
    set(gca,'position',[0.01 0.01 0.98 0.98]);
    scPresent = 0;
else
    scPresent = 1;
end


% DRAW r CIRCLES
scResCirc(r, ChColor, xR, xL);
axis equal;
axis off;


% DRAW x ARCS
scReacArc(x, 1, ChColor, rR, rL);
scReacArc(x, -1, ChColor, rR, rL);

if ~scPresent
    
    % DRAW THE OUTER CIRCLE AND THE X-AXIS
    plot([-1 1],[0 0], ChColor);
    plot(-1:.001:1, sqrt(1-(-1:.001:1).^2), ChColor);
    plot(-1:.001:1,-sqrt(1-(-1:.001:1).^2), ChColor);
    
    % PRINT r VALUES
    scPrnVal(rPrint, 'r');
    
    % PRINT x VALUES
    scPrnVal(xPrint, 'x');
    
    % DRAW THE ANGLE BOUNDARY
    scConCirc(1.01,'r');
    scConCirc(1.08,'r');
    scAngles(1.01, 1.08, 'r');
    
    % DRAW THE LENGTH BOUNDARY
    scConCirc(1.085,'m');
    scConCirc(1.16,'m');
    scLength(1.085, 1.16, 'm');
    
    if FullMapWithLabels
        scScales;
        scLabels;
    end
end

function scResCirc(r, LinCol, xR, xL)
% scResCirc: draws r circles in the smith chart
%
%  SYNOPSIS:
%     This fuction draws r circles which have their centers on the real axis and
%     the radius is always such that they pass through the point (1,0).
%     The function is called by scDraw to do so. 
%     
%  SYNTAX:
%     scResCirc(r, LinCol, xR, xL)
%
%  INPUT ARGUMENTS:
%     r      : a vector consisting of the desired values of r 
%     LinCol : desired color of the circle(s)
%     xR     : a vector containing right termination points for r circles
%     xL     : a vector containing left termination points for r circles
%
%  OUTPUT ARGUMENT:
%     none

% PLOT CIRCLES
for ii = 1:length(r)
    rc = r(ii);
    if xR(ii)~=0
        xc1 = xR(ii);
        [xco1, yco1] = scPOI(rc, xc1);
    else
        xco1 = 1;
    end
    
    if xL(ii)~=0
        [xco2, yco2] = scPOI(rc, xL(ii));
    else
        xco2 = 20;
    end
    if xco1 < xco2
        u      = linspace((rc-1)/(rc+1), xco1, 500);
        vplus  = sqrt((1/(1+rc))^2-(u-rc/(1+rc)).^2);
        vminus = (-sqrt((1/(1+rc))^2-(u-rc/(1+rc)).^2));
        plot(real(u),real(vplus),LinCol);
        hold on;
        plot(real(u),real(vminus),LinCol);
    end
    if xco2 ~= 20
        u      = linspace(xco2, xco1, 200);
        vplus  = sqrt((1/(1+rc))^2-(u-rc/(1+rc)).^2);
        vminus = (-sqrt((1/(1+rc))^2-(u-rc/(1+rc)).^2));
        plot(real(u),real(vplus),LinCol);
        hold on;
        plot(real(u),real(vminus),LinCol);
    end
end

function scReacArc(x,updown,LinCol,rR, rL)
%scReacArc: draws x arcs in the smith chart
%
%  SYNOPSIS:
%     This fuction draws x arcs which have their centers on the the line parallel
%     to the imaginary axis and passing through the point (1,0).
%
%     The function is called by scDraw to do so. 
%     
%  SYNTAX:
%     scReacArc(x, updown, LinCol, rR, rL)
%
%  INPUT ARGUMENTS:
%     x     : a vector consisting of the desired values of x
%     updown: draw arc in the upper(lower) half plane if positive(negative)
%     LinCol: desired color of the arc
%     rR    : a vector containing left termination points for x arcs
%     rL    : a vector containing right termination points for x arcs
%
%  OUTPUT ARGUMENT:
%     none

for jj=1:length(x)
    xc     = x(jj);
    if rR(jj)~=0
        rc1 = rR(jj);
        [xco1, yco1] = scPOI(rR(jj), xc);
    else
        xco1 = 1;
        yco1 = 0;
    end
    
    if rL(jj)~=0
        [xco2, yco2] = scPOI(rL(jj), xc);
    else
        xco2 = 20;
        yco2 = 0;
    end
    if xco1 < xco2 && xc~=0
        % FOR THE ARCS STARTING FROM THE OUTERMOST CIRCLE
        if xc <= 1
            % IF THE 9'O CLOCK POS. OF THE X ARC DOES NOT LIE WITHIN THE SMITH CHART
            uup = [];
            vup = [];
            udn = linspace((xc^2-1)/(xc^2+1), xco1, 500);
            vdn = sign(updown)*(1/xc-sqrt((1/xc)^2-(udn-1).^2));
            
        elseif (1/xc-sqrt((1/xc)^2-(1-1/xc-1).^2)) < yco1
            % IF THE 9'O CLOCK POS. OF THE X ARC  IS LOWER THAN POINT OF INTERSECTION
            uup = linspace(xco1,(xc^2-1)/(xc^2+1),100);
            vup = sign(updown)*(1/xc+sqrt(1/xc^2-(uup-1).^2));
            udn = [];
            vdn = [];
        else
            % IF THE 9'O CLOCK POS. OF THE X ARC  IS HIGHER THAN POINT OF INTERSECTION
            uup = linspace(1-1/xc,(xc^2-1)/(xc^2+1),100);
            vup = sign(updown)*(1/xc+sqrt(1/xc^2-(uup-1).^2));
            udn = linspace(1-1/xc,xco1,500);
            vdn = sign(updown)*(1/xc-sqrt((1/xc)^2-(udn-1).^2));
        end
    end
    if xco2~=20 && xc~=0
        if (1/xc-sqrt((1/xc)^2-(1-1/xc-1).^2)) > yco2
            uup = [];
            vup = [];
            udn = linspace(xco2, xco1, 500);
            vdn = sign(updown)*(1/xc-sqrt((1/xc)^2-(udn-1).^2));
        else
            % IF THE 9'O CLOCK POS. OF THE X ARC  IS LOWER THAN POINT OF INTERSECTION
            uup = linspace(1-1/xc,xco2,100);
            vup = sign(updown)*(1/xc+sqrt(1/xc^2-(uup-1).^2));
            udn = linspace(1-1/xc,xco1,500);
            vdn = sign(updown)*(1/xc-sqrt((1/xc)^2-(udn-1).^2));
        end
    end
    if xc~=0
        plot(real(uup),real(vup),LinCol);
        plot(real(udn),real(vdn),LinCol);
    end
end

function scPrnVal(vPrn, rORx)
%scPrnVal: Prints x and r values on the smith chart
%
%  SYNOPSIS:
%     This fuction prints x and r values on the smith chart.
%
%     THIS function is called internally from the function scDraw.
%     
%  SYNTAX:
%     scPrnVal(vPrn, rORx)
%
%  INPUT ARGUMENTS:
%     vPrn  : value to be printed
%     rORx  : whether an 'r' or an 'x' value is to be printed.
%
%  OUTPUT ARGUMENT:
%     none

if rORx == 'x'
    % PRINT x VALUES
    for jj = 1 : length(vPrn)
        xc   = vPrn(jj);
        [xco1, yco1] = scPOI(0, xc);
        [xco2, yco2] = scPOI(1, xc);
        
        if floor(xc)==xc && length(num2str(floor(xc)))>=2
            PrStr='%0.0f';
        else
            PrStr='%0.1f';
        end
        
        h = text (min(xco1), max(yco1), num2str(xc,PrStr));
        set(h,'color','r', 'rotation', 180*atan2(1-min(xco1), max(yco1)-abs(1/xc))/pi, 'fontsize', 4.5, 'HorizontalAlignment', 'right', 'VerticalAlignment', 'bottom');
        
        h = text (min(xco1), -max(yco1), num2str(xc, PrStr));
        set(h,'color','r', 'rotation', 180+180*atan2(-1+min(xco1), max(yco1)-abs(1/xc))/pi, 'fontsize', 4.5, 'HorizontalAlignment', 'left', 'VerticalAlignment', 'bottom');
        if xc <=1
            h = text (min(xco2), -max(yco2), num2str(xc, PrStr));
            set(h,'color','r', 'rotation', 180+180*atan2(-1+min(xco2), max(yco2-abs(1/xc)))/pi, 'fontsize', 4.5, 'HorizontalAlignment', 'left', 'VerticalAlignment', 'bottom');
            
            h = text (min(xco2), max(yco2), num2str(xc, PrStr));
            set(h,'color','r', 'rotation', 180*atan2(1-min(xco2), max(yco2)-abs(1/xc))/pi, 'fontsize', 4.5, 'HorizontalAlignment', 'right', 'VerticalAlignment', 'bottom');
        end
    end
elseif rORx == 'r'
    % PRINT r VALUES
    for jj = 1 : length(vPrn)
        rc   = vPrn(jj);
        xco1 = (rc-1)/(rc+1);
        yco1 = 0.01;
        [xco2, yco2] = scPOI(rc, 1);
        if floor(rc)==rc && length(num2str(floor(rc)))>=2
            PrStr='%0.0f';
        else
            PrStr='%0.1f';
        end
        h = text (xco1, yco1, num2str(rc,PrStr));
        set(h,'color','r','rotation',90,'fontsize',4.5,'HorizontalAlignment','left', 'VerticalAlignment', 'bottom');
        if rc <= 1
            h = text (xco2, yco2, num2str(rc,PrStr));
            set(h,'color','r','rotation',180*atan2(rc/(1+rc)-xco2, yco2)/pi, 'fontsize',4.5,'HorizontalAlignment','left', 'VerticalAlignment', 'bottom');
            h = text (xco2, -yco2, num2str(rc,PrStr));
            set(h,'color','r','rotation', 180+180*atan2(rc/(1+rc)-xco2, -yco2)/pi, 'fontsize',4.5,'HorizontalAlignment','left', 'VerticalAlignment', 'bottom');
        end
    end
else
    error('scPrnVal.m: The second argument must either be r or x');
end

function scAngles(inr, outr, color)
%scAngles: subdivides and shows angles on smith chart boundary
%
%  SYNOPSIS:
%     Marks angles on the outer periphery of the smith chart. The function is called
%     by scDraw to do so. May not produce desired results if called from the matlab
%     command prompt.
%
%     see also scDraw, scLength, sc
%     
%  SYNTAX:
%     scAngles(inr, outr, color)
%
%  INPUT ARGUMENTS:
%     inr   : radiur of the circle
%     outr  : radiur of the circle
%     color : desired color of the circle
%
%  OUTPUT ARGUMENT:
%     none

for ii=0:2:360
    if floor(ii/10) ~= ii/10
        plot([inr*cos(pi*ii/180), inr*1.015*cos(pi*ii/180)], [inr*sin(pi*ii/180), inr*1.015*sin(pi*ii/180)], 'k');
    else
        plot([inr*cos(pi*ii/180), inr*1.025*cos(pi*ii/180)], [inr*sin(pi*ii/180), inr*1.025*sin(pi*ii/180)], color);
        if ii~=360
            degstr = [num2str(ii),'?];
            h = text(inr*1.04*cos(pi*(ii)/180), inr*1.04*sin(pi*(ii)/180), degstr);
            set(h,'HorizontalAlignment','center','color',color,'rotation',270+ii+2.5,'fontsize',6);
        end
    end
end

function scLength(inr, outr, color)
%scLength: draw length marks on the outer boundary of the smith chart
%
%  SYNOPSIS:
%     Marks angles on the outer periphery of the smith chart
%     
%  SYNTAX:
%     scLength(inr, outr, color)
%
%  INPUT ARGUMENTS:
%     inr   : radiur of the circle
%     outr  : radiur of the circle
%     color : desired color of the circle
%
%  OUTPUT ARGUMENT:
%     none

VecLen = 0:.002:0.5;
for ii= 1:length(VecLen)
    if floor((ii-1)/5) ~= (ii-1)/5,
        plot([inr*cos(pi*(180-100*VecLen(ii)*7.2)/180), 1.015*inr*cos(pi*(180 - 100*VecLen(ii)*7.2)/180)], [inr*sin(pi*(180 - 100*VecLen(ii)*7.2)/180), 1.015*inr*sin(pi*(180 - 100*VecLen(ii)*7.2)/180)], color);
    else
        plot([inr*cos(pi*(180 - 100*VecLen(ii)*7.2)/180), 1.025*inr*cos(pi*(180 - 100*VecLen(ii)*7.2)/180)], [inr*sin(pi*(180 - 100*VecLen(ii)*7.2)/180), 1.025*inr*sin(pi*(180 - 100*VecLen(ii)*7.2)/180)], 'k');
        if VecLen(ii)~=0.5
            h = text(inr*1.04*cos(pi*((180 - 100*VecLen(ii)*7.2))/180), inr*1.04*sin(pi*(180 - 100*VecLen(ii)*7.2)/180), num2str(VecLen(ii),'%0.2f'));
            set(h,'color',color,'rotation',90+3.6-100*VecLen(ii)*7.2,'fontsize',6,'HorizontalAlignment','center');
        end
    end
end

function scScales
%scScales: draws the reflection factor and m scales on the smith chart
%
%  SYNOPSIS:
%     This function draws the reflection factor and m scales on the smith chart
%     The function is called by scDraw to do so. 
%     
%  SYNTAX:
%     scScales
%
%  INPUT ARGUMENTS:
%     none
%
%  OUTPUT ARGUMENT:
%     none


% DRAW BASE LINES
hold on;
yB = 1.50;
xB = 0.05;

plot(xB+[0 1], (yB+0.00)*[1 1], 'r');
plot(xB+[0 1], (yB+0.01)*[1 1], 'r');

% DRAW r DIVISIONS  AND PLACE r TEXT LABELS
for jj = 0:10
    plot(xB+[jj,jj]/10, [yB+0.01, yB+.04],'r');
    if jj ~=10
        plot(xB+[jj+.5,jj+.5]/10,[yB+0.01, yB+0.03],'r');
    end
    h = text(xB+jj/10, yB+0.05,num2str((10-jj)/10,'%1.1f'));
    set(h,'rotation',90,'fontsize',5,'fontname','Times New Roman');
end

% PUT CAPTION TEXT ETC. FOR REFLECTION FACTOR
h = text(xB+0.5, yB+0.18, 'Reflection Factor |r|');
set (h, 'fontname','Times New Roman','fontsize', 8,'fontname','Times New Roman','HorizontalAlignment','center');
h = text(xB+0.15, yB+0.18, '<');
set (h, 'color','b');
plot(xB+[0.15 0.25], (yB+0.18)*[1 1]);


% DRAW REFLECTED POWER DIVISIONS  AND PLACE TEXT LABELS
PVal = [1:-0.1:0.2, 0.1:-0.02:0.02, 0.01 0.0];
for jj = 1: length(PVal)
    plot(xB+[1-sqrt(PVal(jj)), 1-sqrt(PVal(jj))], [yB, yB-0.03], 'r');
    if jj ~=length(PVal)
        plot(xB+[1-sqrt((PVal(jj)+PVal(jj+1))/2), 1-sqrt((PVal(jj)+PVal(jj+1))/2)], [yB, yB-.02], 'r');
    end
    if PVal(jj)>0 && PVal(jj)<0.1
        PrStr = '%1.2f';
    else
        PrStr = '%1.1f';
    end
    h = text(xB+1-sqrt(PVal(jj)), yB-0.04,num2str(PVal(jj),'%1.2f'));
    set(h,'rotation',90,'HorizontalAlignment','right','fontname','Times New Roman','fontsize',5);
end


% PUT CAPTION TEXT ETC. FOR REFLECTED POWER
h = text(xB+0.5, yB-0.19, 'Reflected Power');
set (h, 'fontname','Times New Roman','fontsize', 8,'fontname','Times New Roman','HorizontalAlignment','center');
h = text(xB+0.15, yB-0.19, '<');
set (h, 'color','b');
plot(xB+[0.15 0.25], (yB-0.19)*[1 1]);


% ESTABLISH AND DRAW SCALES FOR m
mValMain = 0:0.1:1;
rUpMain  = abs((mValMain-1)./(mValMain+1));
mValSub  = 0.05:0.1:0.95;
rUpSub   = abs((mValSub-1)./(mValSub+1));

mdBMain  = [0:2:10 15 20 30 40];
rDoMain  = abs((10.^(-mdBMain/20)-1)./(1+10.^(-mdBMain/20)));
mdBSub   = [1 3 5 7 9 11 12 13 14 16 17 18 19 22 24 26 28];
rDoSub   = abs((10.^(-mdBSub/20)-1)./(1+10.^(-mdBSub/20)));

plot(-xB-[0 1], (yB+0.00)*[1 1], 'r');
plot(-xB-[0 1], (yB+0.01)*[1 1], 'r');


% DRAW m DIVISIONS  AND PLACE m TEXT LABELS
for jj = 1:length(mValMain)
    plot(-xB-[rUpMain(jj),rUpMain(jj)], [yB+0.01, yB+.04],'r');
    h = text(-xB-rUpMain(jj), yB+0.05,num2str(mValMain(jj),'%1.1f'));
    set(h,'rotation',90,'fontsize',5,'fontname','Times New Roman');
end
for jj = 1:length(mValSub)
    plot(-xB-[rUpSub(jj),rUpSub(jj)], [yB+0.01, yB+.03],'r');
end

% PUT CAPTION TEXT ETC. FOR REFLECTION FACTOR
h = text(-xB-0.5, yB+0.18, 'm = |u_{min}/u_{max}|');
set (h, 'fontname','Times New Roman','fontsize', 8,'fontname','Times New Roman','HorizontalAlignment','center');
h = text(-xB-0.145, yB+0.18, '<');
set (h, 'color','b','rotation',180);
plot(-xB-[0.15 0.25], (yB+0.18)*[1 1]);

% DRAW mdB DIVISIONS  AND PLACE mdB TEXT LABELS
for jj = 1:length(mdBMain)
    plot(-xB-[rDoMain(jj),rDoMain(jj)], [yB, yB-.03],'r');
    h = text(-xB-rDoMain(jj), yB-0.04, num2str(mdBMain(jj), '%.0f'));
    set(h,'rotation',90,'fontsize',5,'fontname','Times New Roman','HorizontalAlignment','right');
end
plot(-xB-[1,1], [yB, yB-.03],'r');
h = text(-xB-1, yB-0.04, '\infty');
set(h,'rotation',90,'fontsize',5,'HorizontalAlignment','right');


for jj = 1:length(mdBSub)
    plot(-xB-[rDoSub(jj),rDoSub(jj)], [yB, yB-.02],'r');
end

% PUT CAPTION TEXT ETC. FOR mdB
h = text(-xB-0.5, yB-0.19, 'm in db');
set (h, 'fontname','Times New Roman','fontsize', 8,'fontname','Times New Roman','HorizontalAlignment','center');
h = text(-xB-0.75, yB-0.19, '<');
set (h, 'color','b');
plot(-xB-[0.65 0.75], (yB-0.19)*[1 1]);

function scLabels
%scLabels: Draws different standard labels on smith chart
%
%  SYNOPSIS:
%     This function draws different standard labels on smith chart. This function is
%     called by scDraw to to so. If called from the matlab commant prompt, it may not
%     produce the desired results.
%
%     
%  SYNTAX:
%     scLabels
%
%  INPUT ARGUMENTS:
%     none
%
%  OUTPUT ARGUMENT:
%     none

% LABEL FOR INNER SCALE
r  = 1.25;
th = (145:-.1:140)*pi/180;
plot (r*cos(th), r*sin(th),'k');
h = text(r*cos(th(1)), r*sin(th(1)), '<');
set(h,'rotation',(th(1)-pi/2)*180/pi,'VerticalAlignment','middle');

PStr = 'arg(r) Inner Scale';
th = (linspace(140,120,length(PStr)))*pi/180;
for jj = 1: length(PStr)
    h = text(r*cos(th(jj)), r*sin(th(jj)), PStr(jj));
    set(h,'rotation',(th(jj)-pi/2)*180/pi, 'fontsize',9, 'fontname','Times New Roman');
end

% LABEL FOR OUTER SCALE
r  = 1.35;
th = (125:-.1:120)*pi/180;
plot (r*cos(th), r*sin(th),'k');
h = text(r*cos(th(length(th))), r*sin(th(length(th))), '<');
set(h,'rotation',(th(length(th))+pi/2)*180/pi,'VerticalAlignment','middle');

th = [146 145 144]*pi/180;
h = text(r*cos(th(1)), r*sin(th(1)), 'l');
set(h,'rotation',(th(1)-pi/2)*180/pi, 'fontsize',9, 'fontname','mt extra');
h = text(r*cos(th(2)), r*sin(th(2)), '/');
set(h,'rotation',(th(2)-pi/2)*180/pi, 'fontsize',9, 'fontname','Times New Roman');
h = text(r*cos(th(3)), r*sin(th(3)), '\lambda');
set(h,'rotation',(th(3)-pi/2)*180/pi, 'fontsize',9);


PStr = 'Outer scale';
th = (linspace(142,126,length(PStr)))*pi/180;
for jj = 1: length(PStr)
    h = text(r*cos(th(jj)), r*sin(th(jj)), PStr(jj));
    set(h,'rotation',(th(jj)-pi/2)*180/pi, 'fontsize',9, 'fontname','Times New Roman');
end


% PRINT ZL =   AT THE BOTTOM RIGHT CORNER

h = text (1.0, -1, 'Z_L =      \Omega');
set(h, 'fontname', 'Times New Roman', 'fontsize', 14,'HorizontalAlignment', 'right');


% DRAW THE TRANSMISSION LINE SYMBOL AT THE BOTTOM LEFT CORNER

% TL OUTLINE
Lx = -1.00;
Rx = -0.75;
Ty = -0.95;
By = -1.05;
plot([Lx Rx],[Ty Ty], 'r');
plot([Lx Rx],[By By], 'r');

% CONNECTION POINTS
Cx = 0.03;
h=plot([Rx-Cx Rx-Cx Lx Lx],[Ty By Ty By], '.');
set(h,'linewidth',1.5);


% CONNECTIVE LINES
Cy = 0.02;
plot([Rx Rx],[Ty Ty-Cy], 'r');
plot([Rx Rx],[By By+Cy], 'r');

% IMPEDANCE HORIZONTAL LINES
Ix = 0.01;
plot([Rx-Ix Rx+Ix], [By+Cy By+Cy] , 'b');
plot([Rx-Ix Rx+Ix], [Ty-Cy Ty-Cy] , 'b');

% IMPEDANCE VERTICAL LINES
plot([Rx-Ix Rx-Ix], [Ty-Cy By+Cy] , 'b');
plot([Rx+Ix Rx+Ix], [Ty-Cy By+Cy] , 'b');

% LABEL R+jX
h = text (Lx,(Ty+By)/2,'R + jX');
set (h, 'fontname', 'Times New Roman', 'fontsize', 5);

% LABEL l and lines etc.
h = text ((Rx+Lx-Cx)/2, By -0.05,'l');
set (h, 'fontname', 'mt extra', 'fontsize', 5, 'HorizontalAlignment', 'center');
plot([Lx    (Rx+Lx-Cx)/2-0.01],[By-0.05 By-0.05], 'k');
plot([Rx-Cx (Rx+Lx-Cx)/2+0.01],[By-0.05 By-0.05], 'k');
plot([Lx    Lx   ], [By-0.04 By-0.06],'k');
plot([Rx-Cx Rx-Cx], [By-0.04 By-0.06],'k');

Tags:微波 
作者:highspeedlogic 来源:highspeedlogic
  • 您是如何找到本站的?
  • 百度搜索
  • Google搜索
  • 查阅资料过程中
  • 论坛发现
  • 百度贴吧发现
  • 朋友介绍
本站最新成功开发工程项目案例
相关文章
  • 没有相关文章
相关评论
发表我的评论
  • 大名:
  • 内容:
  • MATLAB代做,MATLAB专业代做,MATLAB淘宝代做(www.hslogic.com) © 2018 版权所有 All Rights Reserved.
  • Email:highspeed_logic@163.com 站长QQ: 1224848052

    专业代做/代写/承接、MATLAB、SIMULINK、FPGA项目、博士/硕士/本科毕业设计、课题设计、论文,毕业论文,Coursework、Eassy、Assignment