里兹法的MATLAB程序实现

李兹法的思想首先由Lord Rayleigh提出，后又由W. Ritz独立地从更一般的意义上阐明其思想，因此此法又称为Rayleigh-Ritz法。它无需导出微分方程（即Euler-Lagrange方程）来进行求解，而是直接从泛函出发，通过假设所求驻值函数为已知函数的组合，将变分问题化为函数求极值问题。

function [DC]=RitzTorsionD（abn）
% This function calculates the torsional rigidity of rectangular section
%  based on Ritz method in which the trial function is
%    （x^2-a^2）*（y^2-b^2）*（ C1 + C2*x^2*y^2 + C3*x^4*y^4 + C4*x^6*y^6 + ...）
% Usage:
%       [CD] = RitzTorsion（abn）
% Input:
%       a ------ width of the section
%       b ------ height of the section
%       n ------ the number of the trial functions
% Output:
%       D ------ the torsional rigidity
%       C ------ the coefficients of the trial functions
%
%---------------------------------------------------------------
%  Any comments can be sent to xurongqia