先用外罚函数法转换为无约束的matlab程序

先用外罚函数法转换为无约束优化，再用共轭梯度法求解的matlab程序