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