You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
24 lines
411 B
24 lines
411 B
11 months ago
|
function SensiticityCompare()
|
||
|
|
|
||
|
|
load('TPMS_CHsample.mat', 'CH','cc','num_c')
|
||
|
|
i=1;j=1;
|
||
|
|
CHij = zeros(num_c,1);
|
||
|
|
for ic=1:num_c
|
||
|
|
CHij(ic)=CH{ic}(i,j);
|
||
|
|
end
|
||
|
|
P=polyfit(cc.',CHij,9);
|
||
|
|
% plot(cc,CHij)
|
||
|
|
Sfit = polyval(P(1:9),cc);
|
||
|
|
|
||
|
|
|
||
|
|
|
||
|
|
load('TPMS_CHdc_sample.mat', 'dDHdc')
|
||
|
|
CHdcij = zeros(num_c,1);
|
||
|
|
for ic=1:num_c
|
||
|
|
CHdcij(ic)=dDHdc{ic}(i,j);
|
||
|
|
end
|
||
|
|
|
||
|
|
plot(cc.',[CHij,Sfit.',CHdcij])
|
||
|
|
legend('CHij','fitting','Analytical')
|
||
|
|
|
||
|
|
end
|