随便问问

sys1983 · 发表于 2006-5-19 15:38

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

老师要求用动画表示（最好）界面也可以关于数值分析的 好难下面是代码有高手帮忙看下知道不 function [x,det,index] = Gauss(A,b) % 求线性方程组的列主元 Gauss 消去法，其中， % A 为方程组的系数矩阵； % b 为方程组的右端项； % x 为方程组的解； % det 为系数矩阵A的行列式的值； % index 为指标变量，index=0表示计算失败，index=1表示计算成功。 [n,m] = size(A);nb = length(b); % 当方程组行与列的维数不相等时，停止计算，并输出出错信息。 if n~=m error('The rows and columns of matrix A must be equal!'); return; end %当方程组与右端项的维数不匹配时，停止计算，并输出出错信息 if m~=nb error('The columns of A must be equal the length of b!'); return; end %开始计算，先赋初值 index = 1; det = 1;x=zeros(n,1); for k=1:n-1 % 选主元 a_max = 0; for i = k:n if abs(A(i,k))>a_max a_max=abs(A(i,k));r=1; end end if a_max<1e-10 index = 0;return; end % 交换两行 if r>k for j=k:n z=A(k,j);A(k,j)=A(r,j);A(r,j)=z; end z=b(k);b(k)=b(r);b(r)=z;det=-det; end %消元过程 for i=k+1:n m=A(i,k)/A(k,k); for j = k+1:n A(i,j)=A(i,j)-m*A(k,j); end b(i)=b(i)-m*b(k); end det=det*A(k,k); end det=det*A(n,n); % 回代过程。 if abs(A(n,n))<1e-10 index=0;return; end for k=n:-1:1 for j=k+1:n b(k)=b(k)-A(k,j)*x(j); end x(k)=b(k)/A(k,k); end

happy · 发表于 2006-5-19 18:42

什么意思？你想问什么问题？

sys1983 · 发表于 2006-5-20 09:09

教授你好对上面的线性方程组求解的过程加点代码 我们老师要求用动画表示（最好），如或者做个界面也可以输入矩阵A， b 后面得到结果x的解

happy · 发表于 2006-5-20 14:04

是不是上述函数永动画的形式画图？ 这个前面有一个摆线的例子，你搜索一下，把里边的函数换成你的函数就基本差不多了，可能个别参数需要调整一下

sys1983 · 发表于 2006-5-20 14:19

教授再问一下也是数值分析的 引用了MATLAB的函数ode23 这个列向量为yprime。同样，y1和y2合并写成列向量y。所得函数M文件是： 
function yprime=vdpol(t , y);
% VDPOL(t , y) returns derivatives of the Van der Pol equation:
% 
% x ‘‘-mu *(1-x ^2)*x ‘+x=0 (‘ = d/dx , ‘‘ = d^2/dx^2)
% 
% let y(1)=x and y(2)=x‘
% 
% then y(1) ‘ = y(2)
% y(2) ‘ = MU*(1-y(1)^2)*y(2)-y(1)
 
global MU % choose 0<MU<10 in Command workspace
 
yprime=[y(2) MU*(1-y(1)^2)*y(2)-y(1)]; % output must be a column
 
计算结果如下：
 
>>global MU % define MU as a global variable in the Command Workspace
>>MU=2; % set global parameter to desired value
>>[t , y]=ode23(‘ vdpol ’ , 0 , 30 , [1 ; 0]); % to=0 , tf=30 , yo=[1 ; 0]
>>y1=y( : , 1); % first column is y(1) versus time points in t
>>y2=y( : , 2); % second column is y(2)
>>plot(t , y1 ， t , y2 , ‘ -- ‘)
>>xlabel(‘ Time , Second ‘) , ylabel(‘ Y(1) and Y(2) ‘)
>>title(‘ Van der Pol Solution for mu=2 ‘) 前面是对的按例应该可以出图的（后面有问题） 
[t , y]=ode23(‘ vdpol ‘ , 0 , 30 , [1 ; 0]); % to=0 , tf=30 , yo=[1 ; 0] 我替换成了 [t , y]=ode23(‘vdpol ’, [0 , 30] , [1 ， 0]); 但是后面还是有错 能帮我看下指导下吗？

happy · 发表于 2006-5-20 14:27

函数改为 function yprime=vdpol(t , y); global MU % choose 0<MU<10 in Command workspace yprime=[y(2); MU*(1-y(1)^2)*y(2)-y(1)]; % output must be a column 调用改为 [t , y]=ode23('vdpol' , [0,30] , [1,0]);

sys1983 · 发表于 2006-5-21 11:04

好的谢谢教授

账号		自动登录	找回密码
密码			我要加入

随便问问

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（sys1983）随便问问

回复：（sys1983）随便问问

回复：（sys1983）随便问问