% 读取数据


% 假设有参数和温度的数据data
data = [	
29.2,0.439548;
30,0.512072;
31,0.6235;
32.2,0.776336;
33.1,0.82584;
34.1,0.837972;
35.1,0.917072;
36,0.972324;
37.2,1.006392;
38.9,1.0119;
39.8,1.04762;
40.8,1.11426;
42.7,1.142724;
44.1,1.17678;
45,1.197368;
46.2,1.404088;
47.4,1.510588;
48.1,1.592216;
49.9,1.665832;
52,1.68084;
53.7,1.814832;
54.3,1.932868;
54.6,1.88762;
55,1.924188;
55.2,1.968284;
55.7,1.96602;
56.6,2.016636;
57,2.06644;
57,2.069168;
57.3,2.09716;





]

% 读取数据
% 进行多项式拟合，假设拟合次数为7
degree = 7;
p = polyfit(data(:,2), data(:,1), degree);

% 输出拟合图和多项式表达式
figure;
scatter(data(:,2), data(:,1), 'r');
hold on;
x = linspace(min(data(:,2)), max(data(:,2)), 100);
y = polyval(p, x);
plot(x, y, 'b');
xlabel('ARCN(db)');
ylabel('Temperature(℃)');
title('Waveguide Antenna Pork Experiment');

% 输出多项式表达式为y=x形式
coeff_str = ['y = ' num2str(p(1)) 'x^7 + ' num2str(p(2)) 'x^6 + ' num2str(p(3)) 'x^5 + '...
             num2str(p(4)) 'x^4 + ' num2str(p(5)) 'x^3 + ' num2str(p(6)) 'x^2 + '...
             num2str(p(7)) 'x + ' num2str(p(8))];
fprintf('多项式表达式为：\n%s\n', coeff_str);
%