SaintYZ 发表于 2008-7-18 17:41

请高手指点（关于连续变量的表达式）

以下是程序的一部分：
functionyangzh=model
t=0.25;
......
......
for i=1:4
    p(i)=16*pi/w(i);
    y(i)=-f0(i).*p(i).*sin(w(i).*t)/k0(i,i)/(1-p(i).^2)+f0(i).*sin(16*pi.*t)/k0(i,i)/(1-p(i).^2);
end
disp(y);
我的问题是我不想t 赋初值，我想输出y 是关于t 的函数，我试了几种方法都不可以，请高手指点一下。
不甚感激！！！

[ 本帖最后由 sigma665 于 2008-7-21 20:54 编辑 ]

songzy41 发表于 2008-7-18 20:23

在disp()函数中，变量是自符串，而不是数值。可以改为
disp(num2str(y));

SaintYZ 发表于 2008-7-18 20:33

问题是y 也只是中间的一环，我所要的结果是另外一个矩阵u 也y 的乘积

sogooda 发表于 2008-7-19 09:49

回复 3楼 的帖子

请试着清楚一点儿描述你的问题，还有就是请使用可以准确反映帖子内容的标题。

SaintYZ 发表于 2008-7-19 11:38

t 相当于一个变量，我要得到的结果是关于t 的表达式，像程序中的y 矩阵，是一个其元素是关于t 的表达式。
我最后的结果x 是y 与另一个矩阵相乘得到的。
如果我最开始就定义t 等于一个具体数值的话，那就可以求出y 和x ，但我想要的是连续解，也就是解是关于t 的表达式。

friendchj 发表于 2008-7-20 10:40

自己修改了一下，看是否是你的需要
functionyangzh=myfun1
% t=0.25;
syms t
w=ones(1,4);
f0=w;
k0=2*ones(4);
for i=1:4
    p(i)=16*pi/w(i);
    yangzh(i)=-f0(i).*p(i).*sin(w(i).*t)/k0(i,i)/(1-p(i).^2)+f0(i).*sin(16*pi.*t)/k0(i,i)/(1-p(i).^2);
end
disp(yangzh);
运行结果：
>> y=myfun1
[ 8796093022208/2776947157312239*pi*sin(t)-549755813888/2776947157312239*sin(16*pi*t), 8796093022208/2776947157312239*pi*sin(t)-549755813888/2776947157312239*sin(16*pi*t), 8796093022208/2776947157312239*pi*sin(t)-549755813888/2776947157312239*sin(16*pi*t), 8796093022208/2776947157312239*pi*sin(t)-549755813888/2776947157312239*sin(16*pi*t)]


y =

[ 8796093022208/2776947157312239*pi*sin(t)-549755813888/2776947157312239*sin(16*pi*t), 8796093022208/2776947157312239*pi*sin(t)-549755813888/2776947157312239*sin(16*pi*t), 8796093022208/2776947157312239*pi*sin(t)-549755813888/2776947157312239*sin(16*pi*t), 8796093022208/2776947157312239*pi*sin(t)-549755813888/2776947157312239*sin(16*pi*t)]


>> vpa(y,3)

ans =

[ .995e-2*sin(t)-.198e-3*sin(50.2*t), .995e-2*sin(t)-.198e-3*sin(50.2*t), .995e-2*sin(t)-.198e-3*sin(50.2*t), .995e-2*sin(t)-.198e-3*sin(50.2*t)]


>> subs(y(1),0.25)

ans =

    0.0025
页: [1]
查看完整版本: 请高手指点（关于连续变量的表达式）