% 定义方程
syms wn t
G = 0.8 * (0.36/wn - (1.04)^t * 0.12) * (wn/(0.6 * 0.9^2 * 0.1^0.1))^(-1/0.46) - 0.2;
% 初始化 t 的值范围
t_values = linspace(-10, 10, 20);
% 计算对应的 wn 值
wn_values = zeros(size(t_values));
for i = 1:length(t_values)
t_val = t_values(i);
wn_val = double(solve(subs(G, t, t_val), wn));
wn_values(i) = wn_val(1);
end
% 绘制图形
figure;
plot(t_values, wn_values, 'o-');
xlabel('t');
ylabel('wn');
title('wn vs t');
grid on;