jennychung 发表于 2007-4-8 14:49

sym变量如何显示的问题

运行以下程序
syms h
h1=-1/(4*3/(pi*h)*sqrt(1-(1/h)^2))-7.07;
=solve(h1,h)

结果是:
h =

1.0391
3.6756

我要将h的结果显示在Edit Text中
但是我用
set(findobj('tag','edit12'),'string',h);

显示错误
??? Error using ==> set
Conversion to double from sym is not possible.

Error in ==> feixianxing>pushbutton3_Callback at 607
set(findobj('tag','edit12'),'string',z);

Error in ==> gui_mainfcn at 75
      feval(varargin{:});

Error in ==> feixianxing at 42
    gui_mainfcn(gui_State, varargin{:});

我把h转为h=num2str(h);
也不对
怎么弄啊!

[ 本帖最后由 eight 于 2007-4-27 18:55 编辑 ]

eight 发表于 2007-4-8 16:12

原帖由 jennychung 于 2007-4-8 14:49 发表
运行以下程序
syms h
h1=-1/(4*3/(pi*h)*sqrt(1-(1/h)^2))-7.07;
=solve(h1,h)

结果是:
h =

1.0391
3.6756

我要将h的结果显示在Edit Text中
但是我用
set(findobj('tag','edit12'),'string ...


help subs, the answer is: h = subs(h)

jennychung 发表于 2007-4-8 16:49

subs()不是替换吗
不对啊

eight 发表于 2007-4-8 16:52

原帖由 jennychung 于 2007-4-8 16:49 发表
subs()不是替换吗
不对啊

以下是我运行的结果:


>> syms h
h1=-1/(4*3/(pi*h)*sqrt(1-(1/h)^2))-7.07;
=solve(h1,h)

h =

-1/50*2121^(1/2)/pi*(2121+50*pi)^(1/2)-1/50/pi*(4498641-106050*pi)^(1/2)
-1/50*2121^(1/2)/pi*(2121+50*pi)^(1/2)+1/50/pi*(4498641-106050*pi)^(1/2)


>> vpa(h)

ans =

-26.986864109735927904304945785532
-1.000687246729425731400525817088


>> subs(h)
ans =
-26.9869
   -1.0007


最后的ans是double型

jennychung 发表于 2007-4-8 16:58

我的程序的一部分是
syms h
      w=1/sqrt(T1*T2);
      Gn=K/(j*w*(T1*j*w+1)*(j*w*T2+1));
      U=real(Gn);
      U=vpa(U,4);
      h1=-pi*h/(4*k)*sqrt(1-(e0/h)^2)-U;
      =solve(h1,h);
end
ping1=num2str(w);
set(findobj('tag','edit11'),'string',ping1);
set(findobj('tag','edit12'),'string',h);

我要将h的结果显示在Edit Text中
但是显示错误
不知道怎么弄了

eight 发表于 2007-4-8 17:05

原帖由 jennychung 于 2007-4-8 16:58 发表
我的程序的一部分是
syms h
      w=1/sqrt(T1*T2);
      Gn=K/(j*w*(T1*j*w+1)*(j*w*T2+1));
      U=real(Gn);
      U=vpa(U,4);
      h1=-pi*h/(4*k)*sqrt(1-(e0/h)^2)-U;
      =sol ...

按照我的方法应该没有问题的,你贴的代码不完整,无法运行,因此无法知道出错信息

jennychung 发表于 2007-4-8 17:11

代码还是很长啊
简单一下:

function pushbutton3_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton3 (see GCBO)
% eventdatareserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
hk1=findobj('tag','edit4');K=str2num(get(hk1,'string'));
ht1=findobj('tag','edit5');T1=str2num(get(ht1,'string'));
ht2=findobj('tag','edit6');T2=str2num(get(ht2,'string'));
hlei1=findobj('tag','popupmenu2');lei1=get(hlei1,'value');
axes(handles.axes1);
switch lei1
    case 1
      ks=2;e0=1;%gai
      for A=1.1:0.1:20;
            x=2*ks/pi*(asin(e0/A)+e0/A*sqrt(1-(e0/A)^2));
            y=0;
            plot(-1/x,y,'m');
            hold on
      end
      num=;den=conv(conv(,),);
      G=tf(num,den);
      w1=1.8:0.1:20;
      nyquist(G,w1);hold off;      
      %A w
      syms h
      w=1/sqrt(T1*T2);
      Gn=K/(j*w*(T1*j*w+1)*(j*w*T2+1));
      U=real(Gn);
      U=vpa(U,4);
      h1=-1/(2*ks/pi*(asin(e0/h)+e0/h*sqrt(1-(e0/h)^2)))-U;
      =solve(h1,h);
    case 2      
      .....
    case 3      
      ....
    case 4      
      ....
    case 5      
      ....
end
ping1=num2str(w);z=subs(h);
set(findobj('tag','edit11'),'string',ping1);
set(findobj('tag','edit12'),'string',z);

eight 发表于 2007-4-8 17:15

原帖由 jennychung 于 2007-4-8 17:11 发表
代码还是很长啊
简单一下:

function pushbutton3_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton3 (see GCBO)
% eventdatareserved - to be defined in a future version...


你这个是GUI中某个函数的代码,贴出来也没用,我一样无法运行,你不如自己调试一下,在z=subs(h)那里设置一个断点,然后看出现什么问题了。如果能够自己解决就最好,不能解决的话再把运行到哪一步出错、还有其出错信息贴上来吧

jennychung 发表于 2007-4-8 17:18

设段点在z=subs(h)

显示错误
??? Error using ==> set
Conversion to double from sym is not possible.

Error in ==> feixianxing>pushbutton3_Callback at 607
set(findobj('tag','edit12'),'string',z);

Error in ==> gui_mainfcn at 75
      feval(varargin{:});

Error in ==> feixianxing at 42
    gui_mainfcn(gui_State, varargin{:});

[ 本帖最后由 jennychung 于 2007-4-8 17:23 编辑 ]

eight 发表于 2007-4-8 17:30

不可能吧,你用以下代码试试,看得到的z是否double型的:

syms h
h1=-1/(4*3/(pi*h)*sqrt(1-(1/h)^2))-7.07;
=solve(h1,h);
a=subs(h);
whos a

我的结果:
Name      Size                  BytesClass
a         2x1                        16double array
Grand total is 2 elements using 16 bytes

jennychung 发表于 2007-4-8 17:38

就运
syms h
h1=-1/(4*3/(pi*h)*sqrt(1-(1/h)^2))-7.07;
=solve(h1,h);
a=subs(h);
whos a
的话,是double型的

但是在我的程序里运出来就是不对啊

double型的是不是也是用num2str转成 string型的呢?
edit text 的string属性是要输入string类型的

eight 发表于 2007-4-8 17:41

原帖由 jennychung 于 2007-4-8 17:38 发表
就运
syms h
h1=-1/(4*3/(pi*h)*sqrt(1-(1/h)^2))-7.07;
=solve(h1,h);
a=subs(h);
whos a
的话,是double型的

但是在我的程序里运出来就是不对啊

double型的是不是也是用num2str转成 string型的呢 ...

一般情况下不需要人工转换,double型的变量在text、edit控件上显示的时候默认会转换为string型的

如果上面的测试得到的是double型,那证明你的问题不在那里,你先 clear all 删除所有变量的值,然后 clc 清屏,再设置断点,再运行一次看看吧

jennychung 发表于 2007-4-8 18:01

我不用GUI运行
就在matlab里运行如下代码的话是对的

syms w w1 x y A U Gn h1 h;
K=input('K=')
T1=input('T1=')
T2=input('T2=')
k=2;e0=1;%gai
for A=1.1:0.1:20;
    x=2*k/pi*(asin(e0/A)+e0/A*sqrt(1-(e0/A)^2));
    y=0;plot(-1/x,y,'m')
    hold on
end
num=;den=conv(conv(,),);
G=tf(num,den);
w1=1.8:0.1:20;
nyquist(G,w1)
%A w
Kc=sym(num2str(K));Tc1=sym(num2str(T1));Tc2=sym(num2str(T2));
w=1/sqrt(Tc1*Tc2)
Gn=Kc/(j*w*(Tc1*j*w+1)*(j*w*Tc2+1));
U=real(Gn);
U=vpa(U,4);
h1=-1/(2*k/pi*(asin(e0/h)+e0/h*sqrt(1-(e0/h)^2)))-U;
=solve(h1,h)

但是在GUI里运行
function pushbutton3_Callback(hObject, eventdata, handles)
% hObject    handle to pushbutton3 (see GCBO)
% eventdatareserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
hk1=findobj('tag','edit4');K=str2num(get(hk1,'string'));
ht1=findobj('tag','edit5');T1=str2num(get(ht1,'string'));
ht2=findobj('tag','edit6');T2=str2num(get(ht2,'string'));
hlei1=findobj('tag','popupmenu2');lei1=get(hlei1,'value');
axes(handles.axes1);
switch lei1
    case 1
      ks=2;e0=1;%gai
      for A=1.1:0.1:20;
            x=2*ks/pi*(asin(e0/A)+e0/A*sqrt(1-(e0/A)^2));
            y=0;
            plot(-1/x,y,'m');
            hold on
      end
      num=;den=conv(conv(,),);
      G=tf(num,den);
      w1=1.8:0.1:20;
      nyquist(G,w1);hold off;      
      %A w
      syms h
      w=1/sqrt(T1*T2);
      Gn=K/(j*w*(T1*j*w+1)*(j*w*T2+1));
      U=real(Gn);
      U=vpa(U,4);
      h1=-1/(2*ks/pi*(asin(e0/h)+e0/h*sqrt(1-(e0/h)^2)))-U;
      =solve(h1,h);
end
ping1=num2str(w);z=subs(h);
set(findobj('tag','edit11'),'string',ping1);
set(findobj('tag','edit12'),'string',z);
就是不对啊

eight 发表于 2007-4-8 18:05

原帖由 jennychung 于 2007-4-8 18:01 发表
我不用GUI运行
就在matlab里运行如下代码的话是对的

syms w w1 x y A U Gn h1 h;
K=input('K=')
T1=input('T1=')
T2=input('T2=')
k=2;e0=1;%gai
for A=1.1:0.1:20;
    x=2*k/pi*(asin(e0/A)+e0/A ...


你还是老老实实地在GUI中该函数的第一个语句设置断点,看看出什么问题了吧。我猜你运行的时候没有运行case 1里面的代码

jennychung 发表于 2007-4-8 18:10

应该运行了case 1啊
因为正确显示了w的值啊.
就是不显示h的值,老是报错
??? Error using ==> set
Conversion to double from sym is not possible.

Error in ==> feixianxing>pushbutton3_Callback at 608
set(findobj('tag','edit12'),'string',z);

Error in ==> gui_mainfcn at 75
      feval(varargin{:});

Error in ==> feixianxing at 42
    gui_mainfcn(gui_State, varargin{:});

??? Error while evaluating uicontrol Callback.
页: [1] 2
查看完整版本: sym变量如何显示的问题