dualband_coupling_matrix_synthesis v 1.0

Synthesizes coupling matrix for dual band bandpass filter

function K = cost_function（xOM_topIndepDsRzPsTzPepsilon1）
    
%% cost function （pass the x values from optimization function fminsearch from main）
%% M matrix from x vector

unique_M_top=unique（M_top‘stable‘）;
unique_M_top=unique_M_top‘;

count_indep=0;
for count=1:1:Indep+1 
    [rc] = find（strcmp（[M_top] unique_M_top（1count）））;
    if strcmp（unique_M_top（1count）‘O‘）
        for i=1:1:numel（r）  
            M（r（i1）c（i1））=0;
        end
    else   
        count_indep=count_indep+1; 
        for i=1:1:numel（r）
            M（r（i1）c（i1））= x（1count_indep）;
        end

    end
end

%% q matrix

sum_mii=0;
for i=1:1:O
    sum_mii=sum_mii+M（ii）;
end
qe1=2/（D（12）+1i*sum_mii）;     %reference Wang Nicholson Lancaster
qeN=qe1;
q=zeros（OO）;
q（11）=1/qe1;
q（OO）=1/qeN;

%% A matrix

syms w;
As=q+1i*w*eye（O）-1i*M;

%% reflection and transmission zeros from coupling matrix

M_1=1i*M-q;             % M_1=M‘ calculation
M_2=M_1;
M_2（1:）=[];            
M_2（:O）=[];            % M_2=M‘‘ calculation
sTzM=2*epsilon1*eig（M_2）‘/sqrt（qe1*qeN）;        % transmission zeros

temp=As;
temp（1:）=[];            
temp（:1）=[];
cof11_As=det（temp）;
Fs=vpa（det（As）-2*cof11_As/qe1）;
sRzM= （solve（Fs==0w））‘;                        % reflection zeros

%% cost function

term1=0;
for i=1:1:numel（sRzP）
    term1=term1+（abs（sRzM（1i）-sRzP（1i）））^2;
end

term2=0;
for i=1:1:numel（sTzP）
    term2=term2+（abs（sTzM（1i）-sTzP（1i）））^2;
end

K=double（term1+term2）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1547  2017-11-12 20:10  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\cost_function.m

     文件        717  2017-06-21 21:45  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\coupling_matrix_variables.m

     文件       3805  2017-10-21 11:13  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\gplotdc.m

     文件       3567  2017-10-11 19:55  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\gradient.m

     文件      11978  2017-10-21 01:25  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\graph_gui.m

     文件        481  2017-10-20 13:10  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\group_delay.m

     文件       2730  2017-11-11 16:21  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\input_variables.m

     文件       2398  2017-11-13 04:11  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\iterated_polynomial.m

     文件        334  2017-11-05 13:20  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\iteration.m

     文件      13835  2017-11-13 04:04  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\main.m

     文件        712  2017-11-11 23:05  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\sparameter_coupling_matrix.m

     文件        485  2017-11-11 23:52  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\sparameter_filter_function.m

     目录          0  2017-11-13 18:32  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm

     目录          0  2017-11-11 13:55  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0

     目录          0  2017-11-11 13:55  dualband_coupling_matrix_synthesis v 1.0

     文件       1094  2017-11-13 03:59  dualband_coupling_matrix_synthesis v 1.0\dualband_coupling_matrix_synthesis v 1.0\cm\Untitled.m

----------- ---------  ---------- -----  ----

                43683                    16