26 views (last 30 days)
Show older comments
AYYAPPA RUDRASIMHA on 28 Aug 2024 at 13:06
Answered: Divyajyoti Nayak on 30 Aug 2024 at 7:27
- SOGI_S_Function.zip
When I implement simple SOGI block to generate quadrature signal for given sinusoidla input, output magnitude comming exactly half. SOGI will give two outputs one is alpha and another is beta component. Beta component magnitude and phase is correct but alpha component magnitude is comming exactly half. Same discrete equations when I implemented in matlab function or tf block in simulink its giving correct. Same equations and everything same but when I implemented in s function its giving wrong.
Correct output with MATLAB Function: (Valpha and vin both having same phase and maginitude)
Wrong output with S Function (C code): (Valpha and vin both having same phase and maginitude)
I attached the simulation files of s function, This error is only because of s function related issues. Please check and let me know.
Thanks in advance.
0 Comments Show -2 older commentsHide -2 older comments
Show -2 older commentsHide -2 older comments
Sign in to comment.
Sign in to answer this question.
Answers (1)
Divyajyoti Nayak on 30 Aug 2024 at 7:27
Open in MATLAB Online
Hi @AYYAPPA RUDRASIMHA,
The simulation results using the MATLAB function block and S-function block are different because the C code you have written for the S-function block is not equivalent to the code in MATLAB function block due to a small difference.
In the MATLAB function block, the ‘err’ variable is given the value of the input variable ‘ref’ after the assignment of variables ‘err_n2’and ‘err_n1’.
err_n2 = err_n1 ; %Initially Zero
err_n1 = err ; %Initially Zero because err = 0
err = ref; %The input ref is assigned to err
While in the S-function block, the equivalent variables ‘SOGI_v.vin’, ‘SOGI_v.vin_n1’, and ‘SOGI_v.vin_n2’ are handled differently. SOGI_v.vin is assigned the input variable ‘vin’ before ‘SOGI_v.vin_n1’ and ‘SOGI_v.vin_n2’ have been assigned. This makes the initial state of the variables different and hence leading to the wrong result.
SOGI_v.vin = vin; %Assigned before SOGI_v.vin_n1 and SOGI_v.vin_n2
SOGI_v.d0 = (4 + 2 * K * Ts * w + (w * w * Ts * Ts));
d0 = (4 + (2 * K * Ts * w) + (w * w * Ts * Ts));
SOGI_v.b1 = (-8 + (2 * w * w * Ts * Ts)) / SOGI_v.d0;
SOGI_v.b2 = ((4 - (2 * K * Ts * w) + (w * w * Ts * Ts))) / SOGI_v.d0;
SOGI_v.a_alpha = (2 * K * Ts * w) / SOGI_v.d0;
SOGI_v.a_beta = (K * w * w * Ts * Ts) / SOGI_v.d0;
b1 = (-8 + (2 * w * w * Ts * Ts)) / (4 + (2 * K * Ts * w) + (w * w * Ts * Ts));
b2 = ((4 - (2 * K * Ts * w) + (w * w * Ts * Ts))) / (4 + (2 * K * Ts * w) + (w * w * Ts * Ts));
a1 = (2 * K * Ts * w) / (4 + (2 * K * Ts * w) + (w * w * Ts * Ts));
SOGI_v.vin_n2 = SOGI_v.vin_n1; %Initially zero
SOGI_v.vin_n1 = SOGI_v.vin; %Not initially zero as vin is already assigned
To fix the issue, just write the assignment of ‘SOGI_v.vin’ after assigning ‘SOGI_v.vin_n2’ and ‘SOGI_v.vin_n1’ and it should give the correct result.
0 Comments Show -2 older commentsHide -2 older comments
Show -2 older commentsHide -2 older comments
Sign in to comment.
Sign in to answer this question.
See Also
Categories
Code GenerationEmbedded CoderDeployment, Integration, and Supported HardwareEmbedded Coder Supported HardwareTexas Instruments C2000 Processors
Find more on Texas Instruments C2000 Processors in Help Center and File Exchange
Tags
- s-function
- simulink
- code
- matlab coder
- matlab function
- embedded coder
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
An Error Occurred
Unable to complete the action because of changes made to the page. Reload the page to see its updated state.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
You can also select a web site from the following list
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- Deutsch
- English
- Français
- United Kingdom(English)
Asia Pacific
- Australia (English)
- India (English)
- New Zealand (English)
- 中国
- 日本Japanese (日本語)
- 한국Korean (한국어)
Contact your local office
