Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

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

随便问问

老师要求用动画表示（最好）界面也可以关于数值分析的 好难下面是代码有高手帮忙看下知道不 function = Gauss(A,b) % 求线性方程组的列主元 Gauss 消去法，其中， % A 为方程组的系数矩阵； % b 为方程组的右端项； % x 为方程组的解； % det 为系数矩阵A的行列式的值； % index 为指标变量，index=0表示计算失败，index=1表示计算成功。 = 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）随便问问

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

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

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

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

回复：（sys1983）随便问问

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

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=;%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
>>=ode23(‘ vdpol ’ , 0 , 30 , );%to=0 , tf=30 , yo=
>>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 ‘) 前面是对的按例应该可以出图的（后面有问题） 
=ode23(‘ vdpol ‘ , 0 , 30 , );%to=0 , tf=30 , yo= 我替换成了 =ode23(‘vdpol ’, , ); 但是后面还是有错 能帮我看下指导下吗？

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

回复：（sys1983）随便问问

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

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

好的谢谢教授

页: [1]

声振论坛's Archiver

随便问问

回复：（sys1983）随便问问

回复：（sys1983）随便问问

回复：（sys1983）随便问问