stucture_if.m

基本的源程序

% stucture_if.m
% if判断结构示例
% 求解一元二次方程a*x^2+b*x+c=0的根
% 输入方程的系数
disp('求解一元二次方程的根：');
a=input('输入系数a：');
b=input('输入系数b：');
c=input('输入系数c：');
disp('输出计算结果：')
% 计算判别式
delta=b^2-4*a*c;
% delta>0，两个不同的实数根
if(delta>0)
    x1=(-b+sqrt(delta))/(2*a);
    x2=(-b-sqrt(delta))/(2*a);
    fprintf('两个不同的实根：\n\tx1=%f\n\tx2=%f',x1,x2);
% delta=0，两个相等的实数根
elseif delta==0
    x1=(-b+sqrt(delta))/(2*a);
    fprintf('两个相同的实根：\n\tx1=x2=%f',x1);
% delta<0，两个虚根
else
    real_a=-b/(2*a);
    imag_b=sqrt(abs(delta))/(2*a);
    x1=real_a-imag_b*i;
    x2=real_a+imag_b*i;
    fprintf('两个不同的复根：\n\tx1=%s\n\tx2=%s',num2str(x1),num2str(x2));    
end